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Abstract 


To ensure a smooth and safe operation in railway transport the knowledge of 
the velocity and a track-selective position of the railway vehicles is essential, 
especially due to the increase of digitalisation and automation in railway trans- 
port. For a completely autonomous operating railway vehicle, it is important to 
avoid costintensive modernisation of the railway track infrastructure. Instead, 
methods which can be used directly on the railway vehicle are necessary. The 
accuracy and availability of the existing on-board sensors (e. g. wheel encoders, 
inertial measuring unit, satellite-based positioning systems or optical methods) 
are not sufficient for track-selective positioning of the railway vehicle. Uncer- 
tainties due to multipath propagation, masking, availability or poor lighting 
conditions lead to ambiguities or even a complete loss of position. Slip between 
wheel and rail or wear of the wheels lead to incorrect velocitys and constantly 
growing uncertainties in the travelled distance. So a track selective position 
determination in areas with several parallel tracks is no longer possible. In 
order to avoid the limitations of present sensors, in this work a new sensor con- 
cept is proposed. The concept is based on the ferromagnetic inhomogeneities 
along the railway track, which are mainly caused by the already existing track- 
side infrastructure elements. For example for the rails, rail fastenings, check 
rails, frogs, the Difference-Inductance-Sensor (DIS) was developed. Due to the 
arrangement of the components along the track and the system-related binding 
of the railway vehicle, each rail segment has a unique ferromagnetic signa- 
ture. To determine the velocity of the railway vehicle, the DIS is equipped 
with two detectors covering the ferromagnetic inhomogeneities of the rails by 
means of an electromagnetic measuring method. The velocity of the railway 
vehicle results from the time delay of the signatures and the distance between 
the detectors. At the same time, the ferromagnetic inhomogeneities of the rails 
serve as landmarks. Each track section in the railway network can be uniquely 
identified, thus providing a track-selective method for position determination. 
The procedure was evaluated under controlled laboratory conditions as well as 
on a georeferenced test track. Here, the determination of the velocity and the 
location procedure was tested using only the DIS. The velocity and position 
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deviation of the DIS to the reference sensors has the same order of magnitude 
as the deviation of the two reference sensors from each other. The velocity 
deviation is less than |0.5| m/s, that of the position less than |1.0| m to the 
reference sensors. In addition, the DIS and the ferromagnetic inhomogeneities 
of the railway provide a track-selective position of the rail vehicle. A model 
is used to determine the maximum length of ferromagnetic inhomogeneities 
of 200 m, which is required for a unique comparison with the ferromagnetic 
track map. An optimization procedure is presented, which improves the search 
on extended ferromagnetic track maps and simultaneously reduces the amount 
of data to be stored. With the completely self-sufficient operating DIS, the 
velocity as well as the track-selective position with a ferromagnetic track map 
is provided. Based on the ferromagnetic inhomogeneities, the DIS is able to 
clearly distinguish between the left and right rail of a track segment. In contrast 
to satellite-based localization methods, parallel or closely spaced rail segments 
can be clearly distinguished from each other by means of the ferromagnetic 
signature. This means that track selectivity is no longer limited by the spa- 
tial proximity of the tracks. Compared to the previous methods, the DIS and 
the ferromagnetic inhomogeneities of the railway track ensure a continuous 
and track-selective position determination of the railway vehicle, which are not 
restricted to individual sections in the vicinity of switches or balises. With an 
absolute position accuracy of less than 1.0 m and a velocity deviation of less 
than 0.5 m/s, the DIS meets the requirements of modern traffic management 
and control systems. The DIS is thus able to replace a satellite-based localiza- 
tion method or to complement it by the complementary measurement method. 
Ambiguities in the track occupancy of the railway vehicle, due to uncertainties 
in the position, are omitted. 
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Kurzfassung 


Mit der zunehmenden Digitalisierung und Automatisierung des Schienen- 
verkehrs ist die Kenntnis der Geschwindigkeit und der gleisselektiven Posi- 
tion der Schienenfahrzeuge grundlegend, um einen reibungsfreien und sicheren 
Betriebsablauf zu gewährleisten. Für einen vollständig autonomen Betrieb- 
sablauf sind insbesondere Verfahren erforderlich, die auf dem Schienen- 
fahrzeug zum Einsatz kommen und dadurch kostenintensive Modernisierun- 
gen der streckenseitigen Infrastruktur umgehen. Für eine gleisselektive Posi- 
tionsbestimmung des Schienenfahrzeugs reichen Genauigkeit und Verfüg- 
barkeit der bisherigen fahrzeugseitigen Sensoren (Radencoder, Inertiale Mes- 
seinheit, satellitengestützte Positionierungssysteme oder optische Verfahren) 
nicht aus. Unsicherheiten aufgrund von Mehrwegausbreitung, Verdeckung, 
geringer Verfügbarkeit, Ausfällen oder schlechte Lichtverhältnissen führen 
zu Mehrdeutigkeiten oder einem kompletten Verlust der Position. Schlupf 
zwischen Rad und Schiene oder Verschleiß des Rades führen zu fehler- 
haften Geschwindigkeiten und stetig wachsenden Unsicherheiten in der zurück- 
gelegten Wegstrecke. Insbesondere in Bereichen mit mehreren parallel zueinan- 
der verlaufenden Gleisen ist eine gleisselektive Positionsbestimmung nicht mehr 
möglich. 

Um die Einschränkungen aktueller Sensoren zu umgehen, wird in dieser 
Arbeit ein neues Sensorkonzept vorgeschlagen. Basierend auf den ferromag- 
netischen Inhomogenitäten entlang der Gleisstrecke, die maßgeblich durch 
die bereits vorhandenen streckenseitigen Infrastrukturelemente, darunter die 
Schienen, Schienenbefestigung, Radlenker und Herzstücke, verursacht werden, 
entstand der Differenz-Induktivitäts-Sensor (DIS). Durch die systembedingte 
Bindung des Schienenfahrzeugs und die Anordnung der Bauelemente entlang 
der Gleisstrecke besitzt jedes Schienenelement eine einzigartige ferromagne- 
tische Signatur. Um daraus die Geschwindigkeit des Schienenfahrzeugs zu 
ermitteln, ist der DIS mit zwei Detektoren ausgestattet, die über ein elektro- 
magnetisches Messverfahren berührungslos an zwei unterschiedlichen Posi- 
tionen die ferromagnetischen Inhomogenitäten der Schienen erfassen. Aus 
dem zeitliche Versatz der Signaturen und dem Abstand der Detektoren ergibt 
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sich die Geschwindigkeit des Schienenfahrzeugs. Gleichzeitig dienen die fer- 
romagnetischen Inhomogenitäten der Schienen als Landmarken, mit denen 
jeder Streckenabschnitt im Streckennetz eindeutig identifizierbar ist. Dadurch 
wird ein gleisselektives Verfahren bereitgestellt, um die Position des Schienen- 
fahrzeugs zu bestimmen. Unter kontrollierten Laborbedingungen und auf einer 
georeferenzierten Teststrecke wird die Erstellung der ferromagnetischen Streck- 
enkarten, die Bestimmung der Geschwindigkeit und das Lokalisierungsver- 
fahren unter alleiniger Verwendung des DIS evaluiert. Die Abweichung in der 
Geschwindigkeit und Position des DIS zu den Referenzsensoren liegt dabei 
in derselben Größenordnung wie die Abweichungen der beiden Referenzsen- 
soren zueinander. Die Abweichung der Geschwindigkeit beträgt weniger als 
|0,5| m/s, die der Position weniger als |1, 0| m, zu der Referenzsensorik. Zusät- 
zlich liefert der DIS über die ferromagnetischen Inhomogenitäten der Schiene 
eine gleisselektive Position des Schienenfahrzeugs. Über ein Modell wird eine 
maximale Länge der ferromagnetischen Inhomogenitäten von 200 m ermittelt, 
die für einen eindeutigen Abgleich mit der ferromagnetischen Streckenkarte 
erforderlich ist. Zusätzlich wird ein Optimierungsverfahren vorgestellt, mit dem 
der Abgleich auf ausgedehnten Streckenkarten beschleunigt und gleichzeitig die 
zu speichernde Datenmenge der ferromagnetischen Streckenkarten reduziert 
wird. 

Mit dem vollständig autark operierenden DIS ist es nun möglich, berührungs- 
los die Geschwindigkeit und die gleisselektive Position mit einer ferromag- 
netischen Streckenkarte zu ermitteln. Anhand der ferromagnetischen Inho- 
mogenitäten ist der DIS in der Lage, die linke und rechte Schiene eines 
Gleissegmentes eindeutig voneinander zu trennen. Im Gegensatz zu satel- 
litengestützten Lokalisierungsverfahren sind parallel verlaufende oder nah 
beieinander liegende Schienensegmente eindeutig unterscheidbar. Die Gleis- 
selektivität ist dadurch nicht mehr durch die räumliche Nähe der Gleise 
beschränkt. Im Vergleich zu den bisherigen Verfahren ist mit dem DIS 
und den ferromagnetischen Inhomogenitäten der Schiene eine kontinuier- 
liche und gleisselektive Position des Schienenfahrzeugs gegeben und nicht 
auf einzelne Abschnitte in der Nähe von Weichen oder Balisen beschränkt. 
Mit einer absoluten Positionsgenauigkeit von unter einem Meter und einer 
Geschwindigkeitsabweichung von weniger als 0,5 m/s erfüllt der DIS die 
Anforderungen moderner Verkehrsmanagement- und Kontrollsysteme. Der 
DIS ist dadurch in der Lage, ein satellitenbasiertes Lokalisierungsverfahren 
zu ersetzen oder dieses über das komplementäre Messverfahren sinnvoll zu 
ergänzen. Mehrdeutigkeiten in der Gleisbelegung des Schienenfahrzeugs auf- 
grund von Unsicherheiten in der Position entfallen. 
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1 Einleitung 


Das markanteste Merkmal spurgeführter Fahrzeuge besteht in der Bindung 
des Fahrzeugs an die Fahrbahn. Für Schienenfahrzeuge ist dies durch das 
Gleis — zwei parallel zueinander verlaufende Stahlschienen — gegeben. Wer- 
den Gleise über Weichen miteinander verbunden, entsteht ein ausgedehntes 
Streckennetz, auf dem sich die Schienenfahrzeuge fortbewegen. Aufgrund der 
begrenzten Anzahl an Gleisen und der bisherigen technischen Umsetzung der 
Leit- und Sicherheitstechnik (LST), die zur Absicherung der Gleisstrecke ein- 
gesetzt wird, ist ein erhöhtes Verkehrsaufkommen mit einer Modernisierung 
der Gleisinfrastruktur und somit erheblichen Kosten verbunden. Um auch in 
Zukunft ein erhöhtes Verkehrsaufkommen zu bewältigen, besteht eine Mög- 
lichkeit darin, den Durchsatz auf der Gleisstrecke durch Steigerung der Taktrate 
der Schienenfahrzeuge zu erhöhen. Zur Gewährleistung eines sicheren Betrieb- 
sablaufes werden neue Ansätze, unter anderem für die Geschwindigkeits- 
und gleisselektive Positionsbestimmung der Schienenfahrzeuge, benötigt. Die 
Konzeptionierung und Bereitstellung eines bordautonomen Verfahrens für die 
Geschwindigkeits- und gleisselektive Positionsbestimmung von Schienenfahr- 
zeugen wird in dieser Arbeit realisiert. 


1.1 Motivation 


Der Verkehrsträger Schiene ist ein essenzieller Bestandteil für die Beförderung 
von Personen und den Transport von Gütern. Dies gilt insbesondere für große 
und schwere Nutzlasten, die über weite Strecken hinweg befördert werden. 
Durch die vorgegebene Wegstrecke (Trasse) ist die Haft- bzw. Gleitreibung! 
zwischen den Stahlrädern und der Stahlschiene um den Faktor 10 geringer als 
zwischen Gummireifen und Asphalt. Auf der einen Seite sinkt dadurch der 


l Stahl auf Stahl: Haftreibungszahl uy r = 0,15; Gleitreibungszahl ug pr = 0,03 bis 0,09 [1]. 
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benötigte Energieaufwand für den Transport, erhöht aber gleichzeitig den 
Bremsweg des Schienenfahrzeugs. Fahren auf Sicht ist dadurch im Regelbe- 
trieb nicht möglich. Trotz dieser Einschränkung zählt der Schienenverkehr mit 
2,7 Verletzten und 0,04 Toten pro eine Milliarde Personenkilometer — im Ver- 
gleich zu dem Personenkraftwagen (Pkw) mit 275,8 Verletzten und 2,93 Toten 
pro eine Milliarde Personenkilometer (Stand 2010) in Deutschland — zu 
den sichersten Verkehrsmitteln [2]. Dies wird erreicht, indem streckenseitige 
Sicherungssysteme zum Einsatz kommen. Signale, Balisen, Achszähler, Gleis- 
freimeldeanlagen entlang der Gleisstrecke und das Fahren im Blockabstand? 
gewährleisten einen sicheren Betriebsablauf, indem die Gleisstrecke für jedes 
Schienenfahrzeug abgesichert wird [3]. Über die Länge des Blockabstandes 
und die zulässige Höchstgeschwindigkeit ist der maximale Durchsatz des Stre- 
ckenabschnittes gegeben. Aufgrund der stetig steigenden Verkehrsnachfrage, 
des Ausbaus des multimodalen Mobilitätsangebots, der Reduktion der CO2- 
Emissionen und der Entlastung der Straßen ist eine Zunahme des Schienen- 
verkehrs und eine Verlagerung des Verkehrsaufkommen auf die Schiene erfor- 
derlich. Um nicht den Anschluss zu verlieren, ist es notwendig, die bestehende 
Infrastruktur effizienter zu nutzen, Kosten zu senken und die Zuverlässigkeit 
sowie die Verfügbarkeit über moderne Verkehrsmanagement- und Kontrollsys- 
teme zu regeln [4]. 

Zur Bewältigung des immer stärker anwachsenden Verkehrsaufkommens ist 
es erforderlich, die bestehende Streckeninfrastruktur effizienter auszulasten. 
Eines der erfolgversprechendsten Ansätze stellt dabei der Wechsel vom fes- 
ten Blockabstand auf den wandernden Raumabstand dar. Dabei wandert der 
zu sperrende Streckenabschnitt mit dem Schienenfahrzeug [5] und ermöglicht 
dadurch eine dichtere Taktung der Schienenfahrzeuge. 

Unter Berücksichtigung der zunehmenden Digitalisierung und Automatisie- 
rung des Schienenverkehrs ist die Kenntnis der exakten Geschwindigkeit und 
Position des Schienenfahrzeugs ein essenzieller Bestandteil, um einen siche- 
ren und reibungsfreien Betriebsablauf zu gewährleisten [6]. Auf Basis dieser 
Informationen werden Entscheidungen gefällt und Situationen bewertet, die das 
automatische Schienenfahrzeug anschließend umsetzt. Durch die Verlagerung 
der Sicherungssysteme auf das Schienenfahrzeug ist ein Abbau der strecken- 
seitigen Infrastruktur möglich, wodurch die Kosten bei der Installation und 


? Der Blockabstand oder Streckenblock sperrt den dazugehörigen Streckenabschnitt für die anderen 
Schienenfahrzeuge. 


1.2 Problemstellung 


Instandhaltung gesenkt werden. Um auch in Zukunft einen effizienten und siche- 
ren Betriebsablauf zu gewährleisten, sind neue Verfahren und Ansätze in den 
Bereichen der Geschwindigkeits- und Positionsbestimmung erforderlich. In die- 
ser Arbeit wird ein berührungsloses und bordautonomes Sensorsystem für die 
Geschwindigkeits- und gleisselektive Positionsbestimmung von Schienenfahr- 
zeugen vorgestellt, welches den Anforderungen moderner Zugbeeinflussungs- 
und Zugsicherungssystemen entspricht. 


1.2 Problemstellung 


Eine exakte Geschwindigkeits- und gleisselektive Positionsbestimmung, die auf 
dem Schienenfahrzeug zum Einsatz kommt, ist mit den bisherigen Sensoren 
nicht gegeben und wird aus Mangel entsprechender Alternativen über die stre- 
ckenseitige Infrastruktur (Balisen, LZB, PZB etc.) realisiert. Momentan wird 
über Drehgeber (Radencoder), die an der Motorwelle oder an einer mitlaufenden 
Achse befestigt sind, die Geschwindigkeit und die zurückgelegte Wegstrecke 
auf dem Schienenfahrzeug ermittelt. Aufgrund von Schlupf zwischen Rad und 
Schiene und der Verringerung des Radumfangs durch Verschleiß der Räder ist 
die gemessene Relativgeschwindigkeit zur Trasse und die ermittelte zurück- 
gelegte Wegstrecke mit erheblichen Unsicherheiten verbunden. Infolgedessen 
steigt die Unsicherheit der Position des Schienenfahrzeugs trotz einer vorab 
bekannten Wegstrecke und der Kenntnis der Weicheinstellungen im Gleisnetz 
kontinuierlich an. Eine gleisselektive Position des Schienenfahrzeugs ist bereits 
nach kürzester Zeit nicht mehr gegeben und erfüllt dadurch nicht die Anforde- 
rung moderner Zugbeeinflussungs- und Zugsicherungssysteme [7]. Erst durch 
neue Fertigungstechniken im Bereich der Mikroelektronik und der Massen- 
fertigung entstanden kostengünstige Alternativen, die im Bahnbereich ein- 
gesetzt werden konnten. Satellitengestützte Positionierungssysteme? (GNSS) 
in Verbindung mit Trägheitsnavigationssystemen (IMUs) sind ein Beispiel 
für die fahrzeugseitige Bestimmung der absoluten Position. Unabhängig von 
dem Betriebszustand des Schienenfahrzeugs oder der Weichenstellung werden 


3 Diese umfassen das US-amerikanische GPS, das europäische Galileo, das chinesische BeiDou 
und das russische GLONASS. 
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Abbildung 1.1: Schematische Darstellung der Positionsunsicherheit von satellitengestützten Loka- 
lisierungsverfahren. Eine gleisselektive Position des Schienenfahrzeugs ist nicht gegeben. Die 
Positionshypothesen verteilen sich primär auf die Streckenabschnitte B, C, D und E. 


Abweichungen zu dem vorgegebenen Fahrweg erkannt und folglich die Zuver- 
lässigkeit und Sicherheit im Schienennetz erhöht. Durch die kostengünstige 
Sensorik* stellen GNSS-basierte Verfahren häufig die Grundlage aktueller For- 
schung [8], [9] im Bereich der Selbstlokalisierung von Schienenfahrzeugen dar. 
Mehrwegeausbreitung, Abschattung der Satellitensignale in urbanen Gebieten, 
Wäldern, Bahnhöfen oder in Tunneln beeinträchtigen die Verfügbarkeit und 
erhöhen die Unsicherheiten in der Positionsgenauigkeit und Geschwindigkeit 
des Schienenfahrzeugs. Bereits kleinste Abweichungen, die in der Größenord- 
nung der Messunsicherheit der Sensoren liegen, können zu schweren Unfäl- 
len führen. Selbst die Kombination aus Radencoder, IMU und GNSS sind 
nicht ausreichend, um eine gleisselektive Position und dadurch einen siche- 
ren Betriebsablauf auf Streckenabschnitten mit mehreren parallel zueinander 
verlaufenden Gleisen zu gewährleisten (vgl. Abbildung 1.1). Trotz der vielfälti- 
gen Möglichkeiten und der einfachen Handhabung konnten sich GNSS-basierte 
Lokalisierungsverfahren in der Bahnindustrie bisher nicht durchsetzen [9], [10]. 


1.3 Stand der Forschung 


Um die Schwächen der satellitengestützten Positionsbestimmung zu kompensie- 
ren, werden ortsfeste Strukturen, bspw. Stahlbrücken, Bahnübergänge, Gebäude 
oder benachbarte Gleise, als Landmarken verwendet [11], [12], [13], [14], 


4 Die Entwicklung und Bereitstellung der gesamten Infrastruktur sowie die laufenden Kosten 
werden dabei nicht berücksichtigt. 


1.3 Stand der Forschung 


Tabelle 1.1: Übersicht der gebräuchlichsten Sensoren für die inkrementelle (I) und absolute (A) 
Bestimmung der Position und Geschwindigkeit von Schienenfahrzeugen. Dargestellt wird eine 
qualitative Gegenüberstellung der Positionsgenauigkeit, des Aufwands bei der Inbetriebnahme und 
der Erstellung der Karten, der Verfügbarkeit und der Witterungsbeständigkeit. 


Sensor Beitrag zur Lokalisierung Ink./Abs. Pos.Genaui gkeitAufwandVerfügbarkeitWterungs- 


einfluss 
Radencoder misst Geschwindigkeit/Wegstrecke I gering gering hoch hoch 
IMU misst Beschleunigungen/Drehraten I gering gering hoch hoch 
Doppler-Radarmisst Geschwindigkeit/Wegstrecke T gering gering mittel mittel 
ECS misst Schienenbefestigungen IA) mittel mittel mittel hoch 
GNSS ermittelt geographische Koordinaten A mittel hoch mittel hoch 
LIDAR erkennt Strukturen (Schienen, Weichen) A hoch mittel mittel gering 
Kamera erkennt Strukturen/Objekte A hoch mittel mittel mittel 
Balise stellt Informationspunkte bereit A hoch gering gering hoch 
DIS misst ferr.-mag. Inhomogenitaten A&I hoch gering hoch hoch 


[15]. Die immanenten Eigenschaften der Gleisstrecke, darunter die Krüm- 
mungen, verschiedene Weichenbauteile, Überhöhungen- und Unebenheiten der 
Schienen besitzen eine ausreichende Struktur, um die Position des Schienen- 
fahrzeugs bestimmen zu können [12], [16], [17]. Die wichtigsten fahrzeug- 
seitigen Sensoren, die zur Bestimmung der Geschwindigkeit und der Posi- 
tion für Schienenfahrzeuge verwendet werden, sind u. a. Doppler-Radar [18], 
[19], Light Detection and Ranging (LIDAR) [11], [20], Kamera [15], [21], 
[22], Wirbelstromsensor (ECS) [16], [17], [23], [24], Anisotropic Magneto- 
Resistance (AMR) Sensoren [12], [25], Balisen, Punktförmige/Linienförmige 
Zugbeeinflussung (PZB/LZB) [3], Radio Frequency Identification (RFID)- 
Tags [26], inertiale Messeinheit (IMU) [27], [28], [29], [30], [31]. Eine kontinu- 
ierliche Geschwindigkeits- und gleisselektive Positionsbestimmung auf Basis 
einer flächendeckenden Installation streckenseitiger Infrastrukturelemente ent- 
fällt aufgrund der hohen Installations- und Instandhaltungskosten der Hard- 
ware. PZB/LZB, Balisen oder RFID-Tags entlang der Gleisstrecke stellen des- 
halb keine realistische Alternative für Zukünftige Lokalisierungsverfahren dar. 
Alternativen sind Sensoren und Messverfahren, die direkt auf dem Schienen- 
fahrzeug zum Einsatz kommen und nicht auf streckenseitige Infrastrukturele- 
mente angewiesen sind. Über das Doppler-Radar wird bspw. die Relativge- 
schwindigkeit zwischen Trasse und Schienenfahrzeug bestimmt. Dieser benö- 
tigt jedoch eine Mindestgeschwindigkeit von 16 km/h und ist zudem anfällig 
gegenüber Umwelteinflüssen, wie Vereisung der Trasse oder Spritzwasser [18]. 
Über die Beschleunigungsraten und die Drehlage der IMU oder Achslagerbe- 
schleunigungsaufnehmer (ABA) sind einzelne Streckenabschnitte anhand der 
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jeweiligen Krümmung und ihrer individuellen Vibrationssignaturen voneinan- 
der differenzierbar. Die Ergebnisse hängen dabei stark von dem jeweils einge- 
setzten Fahrzeugtyp ab. Durch die Möglichkeit, den nächsten Streckenabschnitt 
vorausschauend zu untersuchen und dadurch frühzeitig Objekte im Gleis oder 
in der unmittelbaren Umgebung zu erfassen, gehören optische Messverfah- 
ren, darunter Kamera und Light Detection and Ranging (LIDAR), mit zu den 
vielversprechendsten Ansätzen. Die Messverfahren sind jedoch auf gute Sicht- 
verhältnisse angewiesen. Verunreinigungen der optischen Bauelemente durch 
Schmutz, Staub oder Öl, schlechte Lichtverhältnisse oder Umwelteinflüsse, 
bspw. Schnee oder Regen, beeinträchtigen die Funktionsweise der Sensoren. 
AMR und ECS verwenden ein elektromagnetisches Messverfahren, um Gleis- 
und Infrastrukturelemente entlang der Gleisstrecke bzw. Veränderungen des 
Erdmagnetfelds zu erfassen. Zur Ermittlung der Position aus den bereitgestellten 
Sensorinformationen sind bei allen vorgestellten Verfahren Karten erforderlich. 
Abhängig vom Sensor werden die erfassten Merkmale durch ein Referenzel- 
lipsoid (GNSS), eine 3D-Punktwolke (LIDAR), Bildmerkmale (Kamera), eine 
zurückgelegte Wegstrecke (Radencoder und Doppler-Radar), Signaturen oder 
Marker (AMR, ECS, Balisen und RFID-Tags) einer georeferenzierten Position 
in einer Karte zugeordnet. Die Tabelle 1.1 zeigt eine Übersicht der gebräuch- 
lichsten Sensoren für die Bestimmung der Geschwindigkeit und Position für 
Schienenfahrzeuge. 

Die Erfassung der streckenseitigen Infrastruktur über ein elektromagnetisches 
Messverfahren, welches unempfindlich gegenüber Verschmutzung, Witterungs- 
verhältnissen und der Umgebungsbeleuchtung ist, besitzt ein großes Potenzial, 
um neben der Geschwindigkeit auch eine gleisselektive Position des Schie- 
nenfahrzeugs zu bestimmen. Aktuell fehlen hierfür jedoch die entsprechenden 
Sensoren, um dieses Vorhaben in die Tat umzusetzen. 


1.4 Relevante Arbeiten 


Die Erfassung der streckenseitigen Infrastruktur über ein elektromagnetisches 
Messverfahren besitzt ein großes Potenzial, um neben der Geschwindigkeit 
auch eine Position des Schienenfahrzeugs zu ermitteln. Wie bereits erwähnt, 
verwenden die bisherigen Arbeiten u. a. einen Wirbelstromsensor (ECS), um 
grobe Inhomogenitäten, darunter die Schienenbefestigungen entlang der Gleis- 
strecke zu erfassen [16], [17], [32], [33]. Das auffälligste Merkmal des ECS 
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sind die periodischen Sensorsignale, welche von der Anordnung der Schienen- 
befestigungen entlang der Gleisstrecke stammen. Über die Kreuzkorrelations- 
funktion (KKF) wird der Versatz der beiden Sensorsignale bestimmt und mit 
dem Abstand der Sensoren die Geschwindigkeit ermittelt. Die Kreuzkorrelati- 
onsfunktion weist dadurch ebenfalls einen periodischen Verlauf auf, wodurch 
sich das Hauptmaximum und die Nebenmaxima stark ähneln [16]. Treten Stö- 
rungen in den Sensorsignalen auf, ist die Bestimmung des wahren Maximum 
nur bedingt möglich. Nur an Weichen wird die Periodizität der Sensorsignale 
durch die zusätzlichen Bauteile der Weiche unterbrochen. Aus der Klassifikation 
der Bauelemente der Weiche sind diese anhand der Sensorsignale identifizier- 
bar [32]. Neben der Geschwindigkeit ergibt sich aus der spezifischen Anord- 
nung der Schienenelemente (Radlenker, Herzstück etc.) die Befahrrichtung der 
Weiche. Aus der Abfolge der Bauelemente und der zurückgelegten Wegstre- 
cke werden über ein stochastisches Lokalisierungsverfahren Rückschlüsse auf 
die Position des Schienenfahrzeugs innerhalb des Gleisnetzes gezogen [34]. 
Aufgrund der teils großen Unsicherheiten bei der Klassifikation der Bauteile, 
Unsicherheiten in der Geschwindigkeit und der damit verbundenen zurückge- 
legten Wegstrecke (Koppelnavigation) sind mehrere Positionshypothesen des 
Schienenfahrzeugs gegeben. Insbesondere auf sehr großen Streckenkarten und 
ohne die Kenntnis einer initialen Position des Schienenfahrzeugs ist eine ein- 
deutige Position des Schienenfahrzeugs erst durch das überfahren mehrere Wei- 
chen möglich. Darüber hinaus ist durch die geringe Variabilität in den Signalen 
des ECS eine kontinuierliche Positionsbestimmung des Schienenfahrzeugs auf 
freier Strecke nicht gegeben. Zusammenfassend erfüllt der ECS nicht die Anfor- 
derungen, die für eine kontinuierliche und gleisselektive Positionsbestimmung 
des Schienenfahrzeugs erforderlich sind. 


1.5 Aufbau der Arbeit 


Die vorliegende Arbeit behandelt die Konzeptionierung eines neuartigen Sen- 
sors für die bordautonome Bestimmung der Geschwindigkeit und der gleisselek- 
tiven Position von Schienenfahrzeugen. Über ein elektromagnetisches Messver- 
fahren werden dabei die ferromagnetischen Strukturen der Schienen und die der 
unmittelbaren Umgebung erfasst. Durch die Bindung des Schienenfahrzeugs an 
die Gleisstrecke wird diese zu einem essenziellen Bestandteil für die Bestim- 
mung der Geschwindigkeit und Position des Schienenfahrzeugs. Der Aufbau 


7 


1 Einleitung 


des Sensors, das zugrunde liegende Messverfahren, die Funktionsweise der 
Sensorik und die Aufbereitung der Messdaten wird in dem Kapitel 2 beschrie- 
ben. Die Sensorik wird anhand von Simulationen untersucht und ausgehend 
von den Ergebnissen in einem finalen Sensoraufbau umgesetzt. Kapitel 3 geht 
explizit auf die Bestimmung der Fahrtrichtung und der Geschwindigkeit ein. 
Dabei wird ebenfalls die Montageposition des Sensor und die Auswirkung auf 
die Geschwindigkeit untersucht. Anschließend wird über die ferromagnetischen 
Inhomogenitäten der Schiene der Sensor und das Verfahren unter kontrollierten 
Laborbedingungen an verschiedene Szenarien getestet und unter realen Bedin- 
gungen auf einer georeferenzierten Teststrecke validiert. Kapitel 4 behandelt die 
Aufbereitung der ferromagnetischen Inhomogenitäten und die Erstellung der 
dazugehörigen ferromagnetischen Streckenkarten. Mit Hilfe der daraus gewon- 
nen Landmarken wird das Lokalisierungsverfahren über die ferromagnetischen 
Inhomogenitäten der Schienen untersucht und mit den hochgenauen Referenz- 
sensoren verglichen. In der Abbildung 1.2 ist der gesamte Informationsfluss von 
der Signalfilterung, Erkennung der Fahrtrichtung, Erstellung der ferromagneti- 
schen Streckenkarten, Bestimmung der Geschwindigkeit und Position grafisch 
dargestellt. Aus den gemessenen ferromagnetischen Inhomogenitäten wird in 
Kapitel 5 ein Modell der ferromagnetischen Inhomogenitaten erstellt, welches 
die Struktur entlang der Gleisstrecke wiedergibt. Neben der maximalen ferro- 
magnetischen Signallänge, die erforderlich ist, um einen eindeutigen Abgleich 
mit der Streckenkarte zu bekommen, wird ein Verfahren zur Beschleunigung des 
Abgleiches mit der ferromagnetischen Streckenkarte untersucht. Zusammenfas- 
send wird in Kapitel 6 die Leistungsfahigkeit des neuen Sensors in Verbindung 
mit den ferromagnetischen Inhomogenitäten bewertet und ein Ausblick über 
zukünftige Arbeiten gegeben. 


1.5 Aufbau der Arbeit 


Kapitel 3 Kapitel 2 
Geschwindigkeitsbestimmung „u „u 


Sensorkonzept 


Richtung 
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Abbildung 1.2: Grafische Darstellung der Informationsverarbeitungskette der Kapitel 2, 3 und 4. 
Der hellblaue Kasten kennzeichnet die optionalen Sensoren und die gestrichelten Linien kenn- 
zeichnen den optionalen Informationsfluss. Grüne Kästen repräsentieren die Ergebnisse der der 
einzelnen Kapitel. 


2 Sensorkonzept 


In diesem Kapitel wird das Sensorkonzept, das physikalische Messprin- 
zip, der Aufbau und die Funktionsweise des Differenz-Induktivitäts-Sensors 
(engl. Difference-Inductance-Sensor)!, kurz DIS [36], [37], [38], beschrieben. 
Bei dem DIS handelt es sich um einen neuartigen Sensor, der speziell für die 
Anforderungen im Schienenverkehr im Rahmen dieser Arbeit entwickelt wurde. 
Das bordautonome Sensorkonzept des DIS verwendet dabei die immanenten 
ferromagnetischen Inhomogenitäten, die sich entlang und in der unmittelba- 
ren Umgebung der Schiene befinden. Anhand dieser wird die Geschwindigkeit 
und die absolute Position von spurgeführten Fahrzeugen (Schienenfahrzeugen) 
berührungslos bestimmt. 


2.1 Anforderungen an die Sensorik 


Sensoren, die in sicherheitsrelevanten Bereichen, darunter die Geschwindig- 
keits- und Positionsbestimmung, zum Einsatz kommen, sind besonderen 
sicherheitstechnischen Anforderungen unterworfen. Insbesondere in der 
Bahnindustrie müssen alle elektrischen Baugruppen, vom Mikrotonkol- 
ler über die Widerstände und Kondensatoren, für einen erweiterten Tem- 
peraturbereich (-40 °C bis +85 °C?) konzipiert sein. Gleichzeitig müs- 
sen mechanische Anforderungen und die Prüfung der Betriebsmittel für 
Hard- und Software nach (DIN EN 50155) erfüllt werden. Bis ein 


l Patent DE102016105413A1 [35] 
? Übertemperatur von +15 K für die Dauer von 10 Minuten. 
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neuer Sensor für den operativen Betrieb freigegeben wird, muss die- 
ser einen langwierigen Entwicklungs- und Zertifizierungsprozess durchlau- 
fen. Um diese Hürden zu umgehen, wurde bereits bei der Konzepter- 
stellung des DIS darauf geachtet, dass dieser als vollständig autark ope- 
rierendes Sensorsystem verwendet werden kann. Über das berührungslose 
Messverfahren wird der DIS galvanisch von dem Schienenfahrzeug, der rest- 
lichen Bordelektronik und der Infrastruktur getrennt. Dadurch wird verhin- 
dert, dass der DIS in einem bereits zertifizierten Betriebsprozess eingreift und 
bestehende Sicherheitskonzepte auf dem Schienenfahrzeug oder der Streckenin- 
frastruktur gefährdet. Eine erneute Zertifizierung des Gesamtsystems Schie- 
nenfahrzeug unter Berücksichtigung aller Teilsysteme entfällt. Die Zertifizie- 
rung der funktionalen Sicherheit nach DIN EN 50126 (generischer RAMS- 
Prozess) [39] ist zunächst nur für den DIS erforderlich. Um diesen Anfor- 
derungen gerecht zu werden, muss der DIS alle Informationen, die für eine 
Geschwindigkeits- und absolute Positionsbestimmung benötigt werden, aus den 
eigens erfassten Messwerten ableiten können. Informationen, die das Schienen- 
fahrzeug bereitstellt, darunter die Fahrtrichtung, Drehzahl des Motors und die 
Stillstandserkennung, sind für den DIS nicht zugänglich. Als Basis für den DIS 
wird ein energieeffizienter Mikrocontroller STM32F405RG [40] (MCU) ver- 
wendet, der für einen Temperaturbereich von -40 °C bis +85 °C bzw. +105 °C 
ausgelegt ist (entspricht DIN EN 50155). Über das Sensor-Framework wird eine 
Virtuelle Maschine (VM), die ein Bestandteil des European Space Software 
Repository (ESSR) [41] ist, bereitgestellt. Mit der bereits zertifizierten Hard- 
wareabstraktionsschicht (HAL) müssen nur die firmwarespezifischen Kompo- 
nenten des DIS und die Anwendungssoftware zertifiziert werden [42]. 


2.2 Messung ferromagnetischer Merkmalspunkte 
und Landmarken 


Das markanteste Merkmal von spurgeführten Fahrzeugen, insbesondere von 
Schienenfahrzeugen, ist die vorgegebene Trasse (Verkehrsweg). Diese besteht 
aus zwei parallel zueinander verlaufenden Stahlschienen, die in einem festen 
Abstand von 600 bis 1676 mm [43] auf dem Gleisbett liegen. Die Bahnkurve 
des Schienenfahrzeugs wird dadurch, bis auf wenige Millimeter in der Quer- 
richtung, auf eine reine längsseitige Bewegung entlang des Gleises beschränkt. 
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Abbildung 2.1: Während einer Messfahrt aufgenommene und Min-Max-normierte (Inter- 
vall [-1, 1]) Signalverläufe der ferromagnetischen Inhomogenitäten der beiden Detektoren A und B 
für eine Dauer von 20 s. 


Werden mehrere Gleisabschnitte über Weichen miteinander verbunden, entsteht 
ein komplexes Gleisnetz. Durch die starke Bindung des Schienenfahrzeugs an 
das Gleis ist gewährleistet, dass die zurückgelegte Wegstrecke stets dieselben 
ferromagnetischen Merkmalspunkte passiert. Durch die Variation der Bau- 
teile (Schienen, Schienenbefestigungen, Radlenker, Herzstücke etc.) entlang 
der Gleisstrecke und die Superposition aus den ferromagnetischen Struktu- 
ren entsteht eine kontinuierliche Veränderung der ferromagnetischen Zusam- 
mensetzung entlang der Schiene und damit eine einzigartige und eindeutige 
Signatur (Fingerabdruck) der Schiene für jeden Streckenabschnitt. Werden die 
ferromagnetischen Inhomogenitäten, siehe Abbildung 2.1, mit den korrespon- 
dierenden ortsfesten Schienenelementen verbunden, ergeben sich für jedes Stre- 
ckensegment spezifische Landmarken. Die dabei entstehenden Streckenkarten 
beinhalten genügend Informationen, um das Schienenfahrzeug bei einer spä- 
teren Befahrung im Gleisnetz zu lokalisieren. In der Abbildung 2.2 sind die 
ferromagnetischen Strukturen eines Streckenabschnittes farbcodiert zu sehen. 
Dabei ist zu beachten, dass jede Schiene eines Streckenabschnittes eine eigene 
individuelle Signatur besitzt. 


2.3 Messprinzip 


Um den Herausforderungen einer genauen Geschwindigkeits- und absoluten 
Positionsbestimmung auf Basis von ferromagnetischen Inhomogenitäten ent- 
lang der Trasse zu begegnen, ist eine hohe Abtastrate und eine zuverlässige 
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Bestimmung der ferromagnetischen Inhomogenitäten erforderlich. Mit dem 
DIS werden die ferromagnetischen Signaturen an zwei räumlich voneinander 
getrennten Positionen gemessen. Mit der MCU und einer separaten Messschal- 
tung (siehe Abbildung 2.3 und Bestückungsplan Abbildung A.5) für jeden der 
Detektoren werden über ein elektromagnetisches Messverfahren auch verdeckte 
ferromagnetische Strukturen im Gleisbett erfasst. Der DIS ist dadurch in der 
Lage, auch bei schlechten Witterungsverhältnissen, bspw. Schnee oder Regen, 
oder bei starker Verschmutzung, bspw. durch Staub und Öl, die ferromagneti- 
schen Strukturen zu erfassen. Dazu wird jeweils die Induktivität: 


Ni. A 
ls 


L(x) = Ho ` Mr(x) - (2.1) 
an den beiden Detektoren, die aus einer elektrischen Spule bestehen, gemessen. 
Die Induktivität des Detektors hängt dabei von den intrinsischen Parametern 
der Spulengeometrie ab. Dazu gehört die Querschnittfläche A der Spule, die 
Länge der Spule /,, die Anzahl der Windungen N und die magnetische Feldkon- 
stante uo. Der extrinsische Faktor und damit die zu bestimmende Messgröße 
ist die Permeabilitätszahl u,(x). Diese ist von der Position x der ferromagne- 
tischen Materialien entlang der Schiene abhängig [44]. Die Permeabilitätszahl 


Abbildung 2.2: Orthographische 3D-Darstellung der aufgenommenen ferromagnetischen Inhomo- 
genitäten des T4-Rings auf dem Prüf- und Validationscenter Wegberg-Wildenrath (Map data ©2016 
Google Z:32U; E:306076; N:5666477). Die Signalverläufe der ferromagnetischen Inhomogenitäten 
(Laufzeitdifferenz AT") sind farbcodiert in der Karte aufgetragen. 


14 


2.4 Vorzüge des DIS 


wird dabei messtechnisch über die Phase (Phasenverschiebung) %,,; bzw. die 
Laufzeitdifferenz: 


= Apı,i 
~ 2n fo 


bestimmt. Diese beinhaltet die konstante Arbeitsfrequenz fo, die Spannung U 
und Stromstärke / im Stromkreis der Detektorspule. Das theoretische Wirkprin- 
zip und die Funktionsweise ist im Anhang A.1 und A 2 detailliert beschrieben. 
Für eine einheitliche Darstellung der gemessenen ferromagnetischen Inho- 
mogenitäten werden die Laufzeitdifferenzen Min-Max-normiert im Intervall 
von [-1,1] (AT*), vgl. Abbildung 2.1, dargestellt. 


Ar 


(2.2) 


2.4 Vorzüge des DIS 


Für die messtechnische Umsetzung wird eine Schaltung entworfen, deren 
Signalausgabe der Laufzeitdifferenz aus Gleichung 2.2 entspricht (siehe Abbil- 
dung 2.3). Im Vergleich zu den bisherigen Messverfahren verzichtet der DIS bei 
der Signalerfassung auf analoge Baugruppen. Über ein integriertes Gate-Array 
in der MCU ist eine zeitsynchrone Erfassung der Laufzeitdifferenz für mehrere 
Detektoren gewährleistet. Durch die hardwarebasierte und taktsynchrone Mess- 
werterfassung und durch den Verzicht auf Leistungselektronik ist der DIS im 
Gegensatz zum ECS als Low-Power Sensor konzipiert. Durch die Verwendung 
der ferromagnetischen Inhomogenitäten als Landmarken entfällt die Instal- 
lation von kostenintensiven Infrastrukturelementen entlang der Strecke. Mit 
dem kostengünstigen Aufbau des DIS und der energieeffizienten MCU sind 
Anwendungen auf Güterwagons, die keine Stromversorgung besitzen, denkbar. 
Im Nachfolgenden wird die Funktionsweise der wichtigsten Baugruppen der 
Messschaltung für einen der beiden Detektoren erläutert. Der zweite Detektor 
verwendet eine identische Messschaltung. 
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Abbildung 2.3: Elektrischer Schaltplan einer der beiden Messschaltungen (A/B) des DIS. Grüner 
Block: Einkopplung der Messfrequenz. Blauer Block: Bestimmung der Laufzeitdifferenz. Gelber 
Block: Auswertung der Laufzeitdifferenz. 
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Die Funktionsweise der Messschaltung ist dabei in die drei Bereiche unterteilt: 


e Einkopplung der Messfrequenz in den Detektorstromkreis (grün) 
in Abbildung 2.3 


e Bestimmung der Laufzeitdifferenz zwischen der Spanung und 
Stromstärke (blau) in Abbildung 2.3 


e Auswertung der Laufzeitdifferenz (gelb) in Abbildung 2.3 


2.4.1 Einkopplung der Messfrequenz 


In dem Eingangsbereich der Messschaltung (Abbildung 2.3 grüner Block) wird 
eine feste Erregerfrequenz? (an der Leitung des ARM_TAKT) kapazitiv (über 
den Kondensator C1) in den Schwingkreis eingekoppelt. Dadurch gelangt ledig- 
lich der Anteil der Wechselstromkomponente in die Messschaltung. Die Schal- 
tung des DIS wird mit einer Pulsweitenmodulation (PWM?) direkt über die 
MCU betrieben. Durch den gemeinsamen Basistakt ist die Messfrequenz der 
MCU synchron zu der Erregerfrequenz der Messschaltung. Kommt es zu einer 
Verschiebung der Laufzeitdifferenz zu der Erregerfrequenz wird dies direkt von 
der MCU erfasst. 


2.4.2 Bestimmung der Laufzeitdifferenz 


Die Erregerfrequenz entspricht der in Kapitel A.2 beschriebenen Resonanzfre- 
quenz des Serienschwingkreises und wird über den Transformator T1 in den 
Sekundärkreis mit der Detektorspule L1 (siehe Abbildung 2.3) eingespeist. Der 
Transformator trennt den Detektorstromkreis galvanisch von der Messschal- 
tung, wodurch diese vor Überspannung geschützt ist. Das phasenverschobene 


3 Die Erregerfrequenz wird von der MCU generiert. 
* Tastgrad 50:50 zwischen den Werten von 0,0 V und 3,3 V mit einer konstanten Erregerfrequenz. 
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Abbildung 2.4: Links: Signalverläufe der Messschaltung aus der Schaltungsssimulation für ver- 
schiedene Induktivitätsänderungen. Rechts: Detailansicht des Ausganges des Logikgatters. 


Signal an der Primärseite des Transformators wird vom hochfrequenten Rau- 
schen befreit und die negative Halbwelle über eine Schottky-Diode (D1) termi- 
niert. Die Nulldurchgänge des phasenverschobenen Signalverlaufs der verblei- 
benden positiven Halbwelle werden über den Komparator UlA (LM2901) [45] 
gegeniiber GND bestimmt. Am Ausgang des Komparators liegt anschlieBend 
ein PWM-Signal (der positiven Halbwelle) des phasenverschobenen Eingangs- 
signals an (siehe Abbildung 2.3 blauer Block). Eine detaillierte Darstellung 
der Signalverläufe ist in der Abbildung 2.4 (links) zu sehen. Das Erreger- 
signal (ARM_TAKT) und der Ausgang des Komparators (Detektorsignal; 
LM2901_Ausgang) werden über ein Logikgatter (NL7SZ97 [46]; AND-Gatter 
mit invertiertem LM2901_Ausgang) miteinander verglichen. Die Zeitspanne 
des aktiven Zustandes (Pulsbreite) am Ausgang des Logikgatters entspricht der 
Laufzeitdifferenz aus der Gleichung 2.2. Besitzen beide Signalverläufe dieselbe 
Phasenlage, ist der Ausgang des Logikgatters Null. Je größer die Laufzeitdif- 
ferenz an den Detektorspulen, desto größer ist die Pulsbreite am Ausgang des 
Logikgatters. 
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Abbildung 2.5: Oben: Zählerbasierte Bestimmung der Pulsbreite (Laufzeitdifferenz Ar, ) für einen 
Detektor am Ausgang des Logikgatters (unten) bei einer internen Abtastfrequenz von 84 MHz. 


2.4.3 Auswertung der Laufzeitdifferenz 


Der Ausgang des Logikgatters wird direkt mit der MCU verbunden. Dort wird 
mit einem internen Zähler (fz = 84 MHz) die Pulsbreite ermittelt. Mit einer 
Erregerfrequenz von fo = 30 kHz ergibt sich mit der Gleichung 2.3 und 2.4 eine 
Auflösung der Phase von Ymin ~ 0,13° bzw. der Laufzeitdifferenz Atmin X 12 ns. 
In der Abbildung 2.5 ist das Verfahren zur Bestimmung der Pulsbreite Ar,, über 
einen der internen Hardwarezähler der MCU zu sehen. 


1 360° 

Ymin = Ze (2.3) 
ke 
Zy =1 

Atmin = ES (2.4) 


Für den eingeschwungenen Zustand wird die Laufzeitdifferenz einmal pro Peri- 
ode über den Zählerstand Z,, (proportional zu Ar,„) ermittelt und zusammen mit 
einem internen Zeitstempel in einen Ringpuffer in der MCU abgelegt. Die 
Messfrequenz des DIS ist identisch mit der Erregerfrequenz der Messschal- 
tung. Ist der Ringpuffer zur Hälfte gefüllt, werden die Zählerstände und die 
Zeitstempel an einen Auswerterechner übertragen. Neben den Laufzeitdiffe- 
renzen, die durch die ferromagnetischen Inhomogenitäten an der Detektorspule 
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Abbildung 2.6: Phase zwischen Spannung und Stromstärke für verschiedene Induktivitätsände- 
rungen an einer Detektospule für den Bereich der Arbeitsfrequenz der Messschaltung. 


verursacht werden, enthält die Pulsbreite am Ausgang des Logikgatters einen 
konstanten Offset. Dieser wird durch die Gatterlaufzeiten der übrigen elektri- 
schen Baugruppen in der Messschaltung verursacht und ist konstant für eine 
feste Erregerfrequenz. Vor der weiteren Verarbeitung wird der Offset von dem 
eigentlichen Nutzsignal der Detektoren entfernt. 


2.5 Ergebnisse der simulierten Messschaltung 


Die Messschaltung für einen der beiden Detektoren, welche in der Abbil- 
dung 2.3 dargestellt ist, wird mit einer SPICES (Simulation Program with 
Integrated Circuit Emphasis) basierenden Schaltungssimulation untersucht. 
Anhand der Simulation wird der gesamte Messprozess, von der Änderung 
der Detektorspuleninduktivität über den Signalvergleich bis zur Ermittlung 


5 LTspice XVII des Halbleiterherstellers Linear Technology (Seit 2017: Analog Devices). 


20 


2.5 Ergebnisse der simulierten Messschaltung 


Noo o - == un ER a — 


Abbildung 2.7: Montageposition des DIS an einem Schienenfahrzeug (gang car) 10 cm tiber dem 
Schienenkopf. Der rote und blaue Kreis markieren die Detektoren A und B. Der gelbe einge- 
färbte Bereich beinhaltet die Messschaltung des DIS. Die ferromagnetischen Inhomogenitaten sind 
schematisch unterhalb des DIS entlang der Schiene aufgetragen. 


der Laufzeitdifferenz, unter Beriicksichtigung aller elektrischer Komponenten 
untersucht. Mit einer AC-Frequenzanalyse® wird zunächst die Arbeitsfrequenz 
der Messschaltung ermittelt. Eine minimale Anderung an der Spuleninduk- 
tivität AL führt dabei zu einer möglichst großen Phasenverschiebung an der 
Detektorspule. In Abbildung 2.6 ist die Phasenlage an der Detektorspule für 
unterschiedliche Verschiebungen aus der Grundinduktivität bei einer Arbeits- 
frequenz von 30 kHz zu sehen. Eine detaillierte Darstellung der Arbeitsfre- 
quenz (ARM_TAKT), des Ausgangs des Komparators und des Logikgatters 
findet sich in der Abbildung 2.4. Die rot- und blau gestrichelten Linien kenn- 
zeichnen dabei den Zeitpunkt einer steigenden bzw. einer fallenden Signalflanke 
des ARM_TAKT. Die schwarz gestrichelten Linien markieren den Zeitpunkt 
einer steigenden Flanke am Ausgang des Komparators. Der Zeitraum zwischen 
einer rot- und einer schwarz gestrichelten Linie entspricht der Laufzeitdifferenz 


6 Untersuchung des Schaltungsverhaltens als Funktion der Frequenz. 
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Abbildung 2.8: Normierte Häufigkeit der Zeitdifferenz zwei aufeinanderfolgender Messungen des 
DIS ausgewertet über drei reale Messfahrten mit insgesamt 169 Millionen Messpunkten [38]. 


(siehe Gleichung 2.4), welche durch die Detektorspule verursacht wird. In der 
Abbildung 2.4 (rechts) ist die Laufzeitdifferenz am Ausgang des Logikgatters in 
einer Detailansicht zu sehen. Abhängig von der Induktivität der Detektorspule 
ändert sich die Laufzeitdifferenz, welche direkt von der MCU (siehe Kapi- 
tel 2.4.2) erfasst wird. Dadurch können über die Messschaltung auch geringe 
Veränderungen in den ferromagnetischen Strukturen erfasst werden. 


2.6 Ergebnisse der realen Messschaltung 


Ausgehend von den Ergebnissen aus den Simulationen und der Schaltung 
aus Abbildung 2.3 wird jeweils eine finale Messschaltung für die Detekto- 
ren (A und B) aus konkreten elektrischen Baugruppen aufgebaut und die beiden 
Ausgänge der Logikgatter (J4) mit der MCU (Pin AO und A3 [40]) verbunden. 
Beide Detektoren besitzen einen Abstand von d = 0,60 m zueinander und wer- 
den zusammen mit der restlichen Elektronik auf einer Trägerplattform mon- 
tiert. Der fertige Aufbau des DIS ist in der Abbildung 2.7, montiert an einem 
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Abbildung 2.9: Normierte Häufigkeitsverteilung für eine Messung bestehend aus 7,2 Millionen 
Messpunkten an einer festen Position. Für den Detektor A beträgt die mittlere Laufzeitdifferenz 
Ha = 51,5 ns, o, = 7,7 ns und für den Detektor B up = 52,0 ns, og = 8,1 ns [38]. 


Erprobungsfahrzeug’, zu sehen. Das Testen der Firmware und der Messschal- 
tung der Detektorspulen soll Aufschluss tiber die Messgenauigkeit des DIS 
geben. Insbesondere die zeitdiskrete und liickenlose Erfassung der Messwerte 
ist für die spätere Geschwindigkeits- und absolute Positionsbestimmung ent- 
scheidend. Durch die zeitdiskrete Messwerterfassung tiber die MCU besitzt 
jeder der Messpunkte einen internen Zeitstempel mit einer zeitlichen Auf- 
lösung von 1 us. In der Abbildung 2.8 sind die Zeitdifferenzen von insge- 
samt 169 Millionen Messpunkten zu sehen. Bedingt durch die Auflösung des 
internen Zeitstempels besitzen 2/3 der Messpunkte eine Zeitdifferenz von 33 us 
und 1/3 aller Messpunkte eine Zeitdifferenz von 34 us. Anhand der Zeitstempel 
ergibt sich eine kontinuierliche und verlustfreie Erfassung der Messpunkte, die 
der zuvor eingestellten Arbeitsfrequenz von 30 kHz entspricht. Davon profi- 
tiert insbesondere das nachfolgende Filter- und Auswerteverfahren, welches für 
die Geschwindigkeits- und absolute Positionsbestimmung verwendet wird. Für 
Geschwindigkeiten von 400 km/h wird jeder Schienenmeter mit 270 Messpunk- 
ten abgetastet. Um das Rauschverhalten der beiden Detektoren zu untersuchen, 


7 Gang car der Firma ARIC (Automotive and Rail Innovation Center). 
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wurden die Detektoren in einem Abstand von 10 cm über dem Schienenkopf 
an ein Schienenfahrzeug montiert. In der Abbildung 2.9 ist die Verteilung der 
Messwerte bei einer festen Position für ca. 7,2 Millionen Messpunkte zu sehen. 
Für den Detektor A ergibt dies eine mittlere Laufzeitdifferenz von ua = 51,5 ns 
mit einer Standardabweichung von oa = 7,7 ns. Die mittlere Laufzeitdiffe- 
renz des Detektor B beträgt ug = 52,0 ns mit einer Standardabweichung von 
op = 8,1 ns. 


2.7 Signalfilterung 


Bei der Konzeption des DIS wurde auf eine spezielle hardwareseitige Signal- 
filterung und Vorverarbeitung der Messdaten verzichtet. Die Aufbereitung und 
Interpretation der mit dem DIS erfassten Messdaten wird auf einem separaten 
Auswerterechner durchgeführt. Dadurch kann die Hardware des DIS flexibler 
für den Mess- und Auswerteprozess eingesetzt werden. Ohne den Aufbau der 
Messschaltung zu verändern, können verschiedene Verfahren für die Analyse 
der Messsignale untersucht werden. In der Abbildung 2.10 ist das normierte 
Messsignal (schwarze Kurve) eines Detektors des DIS für eine Zeitspanne von 
100 s (entspricht 3 Millionen Messpunkten) zu sehen. Der dargestellte Signal- 
verlauf setzt sich aus einem hochfrequenten Quantisierungsrauschen, verur- 
sacht durch die Messeinheit der MCU, und einem niederfrequenten Nutzanteil 
zusammen. Für die weitere Auswertung wird zunächst der Anteil des Rauschens 
aus dem Messsignal (siehe Abbildung 2.10 magenta Kurven) herausgefiltert, 
ohne dabei markante Signalstrukturen (bspw. Kanten oder abrupte Signalän- 
derungen) zu glätten. Hierfür wird ein auf Wavelets basiertes Verfahren zur 
Filterung und Analyse der Messsignale verwendet. Wavelets wurden bereits 
erfolgreich in der Infrarot (IR) Spektroskopie [47], [48] und bei der Auswer- 
tung von Elektroenzephalogrammen (EEG) [49], [50] angewandt. Mit der Zeit- 
und Frequenzdarstellung der Signalstrukturen bieten die Wavelet-Koeffizienten 
Möglichkeiten für die Signalanalyse, die bereits in den genannten Teilgebieten 
erfolgreich umgesetzt wurden. Dazu gehört die Entfernung von Artefakten, die 
Verringerung des Rauschanteils, die Datenreduktion sowie die gezielte Extrak- 
tion von Merkmalen aus dem Orts- bzw. Frequenzbereich des Nutzsignals. 
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Abbildung 2.10: Min-Max-normiertes [-1, 1] Messsignal eines Detektors dl (schwarz) und mit 
der Gleichung 2.20 gefiltertes Signal w (cyan). Rekonstruierte Signalverläufe der Detailkoeffizi- 
enten D; (magenta) und Approximationskoeffizienten A, (rot) für die Level / = {1, 4, 7} [36]. 


2.7.1 Wavelet-Transformation 


Bei der Wavelet-Transformation handelt es sich um eine lineare Transforma- 
tion, ähnlich der Fourier-Transformation, zur Uberfiihrung einer Funktion vom 
Ortsbereich (bzw. Zeitbereich) in den Frequenzbereich. Die Orts- und Fre- 
quenzdarstellung sind äquivalent, besitzen also denselben Informationsgehalt 
und können ohne Informationsverlust ineinander überführt werden. Bei der 
Fourier-Transformation führt eine lokale Änderungen eines Funktionswertes 
im Orts- oder Frequenzbereich zu einer globale Änderung im jeweiligen ande- 
ren Darstellungsbereich. Da im Ortsbereich keinerlei Frequenzinformation und 
im Frequenzbereich keine Ortsinformationen vorhanden sind, können Orts- und 
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Frequenzinformationen bei der Fourier-Transformation nur bedingt als Merk- 
male verwendet werden. Im Vergleich zur Window- bzw. Short-Time-Fourier- 
Transformation (STFT) [51], bei der örtlich- bzw. zeitlich begrenzte Segmente 
untersucht werden (dadurch ist eine lokale Frequenzanalyse des Signals mög- 
lich), bietet die Wavelet-Transformation eine Alternative Ort/Zeit-Frequenz- 
Analyse für die Untersuchung von nichtstationären Signalen. Dabei handelt es 
sich um eine Multiskalen-Analyse, bei der das Signal auf verschiedenen Skalen 
(hier Level /) unterschiedlicher Auflösungen durch die Wavelet-Koeffizienten 
repräsentiert wird [52]. Die Besonderheit der Wavelet-Transformation liegt 
dabei in der gleichzeitigen Analyse der Orts- und Frequenzanteile eines Signals. 


2.7.2 Diskrete Wavelet-Transformation 


Für die Signalanalyse werden die Messpunkte v von einem der Detektoren des 
DIS über digitale Filterbänke (Tief- und Hochpassfilter) der Diskreten Wavelet- 
Transformation (DWT) in unterschiedliche Frequenzbänder zerlegt [53]. Das 
Ergebnis der DWT hängt dabei entscheidend von der Wavelet-Filter-Basis 
(Haar-Wavelet, Daubechies-Wavelets, etc.), d. h. von den Skalierungskoeffi- 
zienten h, und den Wavelet-Koeffizienten gn, ab. Die Anzahl der Waveletko- 
effizienten für die jeweilige Wavelet-Filter-Basis ist über den Index n gegeben. 
Durch das Anwenden der Wavelet-Transformationsmatrix: 


hy ho hz h4 hs he 0 0 
81 82 83 84 85 86 0 0 
0 0 h h hi, hy 0o 0 
0 O gi g2 93 8 0 0 
W; = 0 0 0 0 hy ho 0 0 (2.5) 

000 0g e 0 0 

2 00 
0 0 0 0 0 0 h;-ı hz 
0 0 0 0 0 0 82-1 gz 
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auf die Messpunkte v (mit der Gesamtanzahl der Messpunkte j) des jeweiligen 
Levels / ergeben sich die Approximationskoeffizienten a! und Detailkoeffizien- 
ten d! für das Frequenzband des nächsten Levels: 


vj d 
+1 
V5 d 
l Vi al! L+1) ql+1 1+1)ql+1 
Wr =Wiı| u" |= Au = (a |d)= (v ld”) (2.6) 
k+1 k 
l IH 
Zi Dë 
v d. 
mit: 
alt! = (a, ie, ay (2.7) 
alt! = (as. SER sa), (2.8) 


Abhängig von der Wavelet-Filter-Basis ergibt sich die Anzahl der 
Approximations- und Detailkoeffizienten aus der Anzahl der Waveletkoeffi- 
zienten n und der Anzahl der Messpunkte j über: 
j+n-1 | 


5 (2.9) 


Aus den Wavelet-Koeffizienten a"! und d’*! wird der vorherige Signalzu- 
stand (v!) über: 


Rull = yi, 0 < L< lmax (2.10) 


mit der Rücktransformationsmatrix R; verlustfrei rekonstruiert. Für die 
Wavelet-Transformationsmatrix und die Rücktransformationsmatrix gilt: 


WR = ¥. (2.11) 


Durch die Detailkoeffizienten werden feine Strukturen (sogenannte Details), 
bspw. hochfrequente Signalanteile (u. a. Rauschen der Messeinheit in der 
MCU), dargestellt. Die Approximationskoeffizienten hingegen beinhalten 
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Abbildung 2.11: Darstellung der Multiskalen-Analyse für einen synthetischen Signalverlauf dl. 
Links: Abgebildet sind die Approximationskoeffizienten a! und Detailkoeffizienten d der DWT. 
Rechts: Darstellung der riicktransformierten (IDWT) Approximationskoeffizienten A! und Detail- 
koeffizienten D! jeweils für die Level / = {1, . . ., 4}. Die Skalierung der Achsen ist für die Darstel- 
lungen der Approximations- und Detailkoeffizienten der DWT und IDWT konstant. 
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niedrige Frequenzanteile, welche den ferromagnetischen Signalstrukturen 
im Ursprungssignal entsprechen. Durch rekursives Anwenden der Wavelet- 
Transformation auf die Approximationskoeffizienten des letzten Zustandes wer- 
den die Wavelet-Koeffizienten der tieferen Level / (vgl. Gleichungen 2.12 
bis 2.15) ermittelt. 


y? Wo (a'|d') (2.12) 
Ww 

al 1 (a?|d?) (2.13) 

a2 W2 (aè EN (2. 14) 

al Wis WEIER EC ES (2.15) 


In der Abbildung 2.11 ist die schematische Abfolge der Multiskalen-Analyse 
für die Level / = {1,...,4} zu sehen. Die linke Seite zeigt die Approximations- 
und Detailkoeffizienten (a und d) aus der DWT. Auf der rechten Seite sind 
die rücktransformierten (Inversen Diskreten Wavelet-Transformation (IDWT)) 
Approximations- und Detailkoeffizienten (A und D) zu sehen. Das Ursprungs- 
signal wird durch die Wavelet-Koeffizienten in immer feinere Frequenzbän- 
der fa,ı und fa,ı (siehe Gleichung 2.16 und 2.17) unterteilt. 


fal = 0 4 (2.16) 
fa = $ = (2.17) 
Js = h (2.18) 


Die Messfrequenz fo des DIS liegt bei 30 kHz und f, ist die Nyquist-Frequenz. 
Niedrigere Level decken den Bereich mit höheren Frequenzen ab. Durch die 
überlappenden Frequenzbereiche der Approximations- und Detailkoeflizienten 
wird eine hohe Auflösung im Zeit- und Frequenzbereich des Ursprungssignals 
erreicht. In der Abbildung 2.10 sind die mit der IDWT rücktransformierten 
Signalanteile der Detailkoeffizienten D’ (magenta) und Approximationskoeffi- 
zienten A! (rot) für die Level / = {1,4,7} zu sehen. Das unterteilen in immer 
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feinere Approximations- und Detailkoeffizienten endet sobald das maximale 
Level von: 


Imax = to, (= i) (2.19) 


erreicht ist. Für eine vollständige und verlustfreie Rekonstruktion (IDWT) 
des Ursprungssignals wird nur der Approximationskoeffizient eines Levels 
zusammen mit den Detailkoeffizienten der darunterliegenden Level benötigt, 
v? = IDWT(a'mx, dies. gms! d?, d1). Werden vor der Rücktransforma- 
tion die Detailkoeffizienten dl eines Level / unterdrückt, verschwindet das dazu- 
gehörige Frequenzband f.,ı unabhängig von der Amplitude aus dem Ursprungs- 
signal [48]. Einzelne Frequenzbänder werden dadurch unterschiedlich stark 
gewichtet. Für die Messsignale des DIS wird ein soft-thresholding (wave- 
let shrinkage) nach Donoho [54] auf die Deatilkoeffizienten (siehe Glei- 
chung 2.20 bis 2.23) angewendet [55], [56]. 


di! = sgn(d;) - fson(d!) (2.20) 


By 
= 
ô= o - V2 -Iog0) (2.22) 


0 III 6 <0 
Fi d! = l 
fson(d;) Wi ld'|-5 >0 


(2.21) 


(2.23) 


Das soft-thresholding hängt vom Betrag der Detailkoeffizienten ab und wirkt 
auf den gesamten Frequenzbereich des Signals. Der Betrag der Wavelet- 
Koeffizienten ist ein Maß für den Informationsgehalt. Kleine Beträge werden 
dabei vernachlässigt, ohne dabei die Signalcharakteristika des Ursprungssi- 
gnals zu verfälschen. Höhere Frequenzanteile, bspw. das Quantisierungsrau- 
schen, werden dadurch stärker gedämpft. Die Frequenzanteile des Nutzsignals 
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hingegeben bleiben erhalten. Durch Rekonstruktion der gefilterten Wavelet- 
Koeffizienten: 


8 Rom: 
(ads, d'De, d'De, dng? (2.24) 


wimax! Nmax-1 


wie? (2.25) 


Wow (2.26) 


verschwinden die hochfrequenten Signalanteile aus dem Messsignal. Durch 
das soft-thresholding ergeben sich für die Rekonstruktion im allgemeinen glat- 
tere Signalverläufe und markante Signalstrukturen bleiben weiterhin erhalten, 
wobei sehr feine Details im Ursprungssignal verloren gehen können. In der 
Abbildung 2.10 (cyan Kurve) ist der geglättete Signalverlauf des Ursprungssi- 
gnals zu sehen. 


2.7.3 Wavelet-Koeffizienten 


Anhand der Wavelet-Koeffizienten (a!|d!,...,d!) und der dazugehörigen 
Wavelet-Filter-Basis ist das Ursprungssignal des DIS mit allen Details voll- 
ständig rekonstruierbar. Wird für die Rücktransformation ausschließlich einer 
der Approximationskoeffizienten a’ verwendet, ergeben sich die Signalverläufe 
A! (vgl. Abbildung 2.11). Mit jedem Level reduziert sich das Signalrauschen 
und die zu speichernde Datenmenge der Approximationskoeffizienten a’. Durch 
den geringeren Detailgrad und die weiterhin gegebenen Ortsauflösung können 
auch die Approximationskoeffizienten, ohne vorherige Rücktransformation, für 
die Analyse der Signalstrukturen und insbesondere für einen Signalabgleich 
verwendet werden. In der Abbildung 2.10 sind für einen der Detektoren des 
DIS die Signalverläufe für die Approximationskoeffizienten al. af und a’ und 
die dazugehörigen Signalverläufe der Rücktransformationen A', AT und A’ 
zu sehen. Verglichen mit dem Wavelet-gefilterten Signalverlauf w° besitzt der 
Signalverlauf A’ genügend Struktur, um eine Geschwindigkeit oder absolute 
Positionsbestimmung durchzuführen. Für die Filterung der Messsignale und zur 
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Repräsentation der Approximations- und Detailkoeffizienten werden Daube- 
chies Wavelets (db Wavelets) verwendet. Daubechies Wavelets eignen sich ins- 
besondere zur Darstellung von Polynomstrukturen, wie sie innerhalb der gemes- 
senen ferromagnetischen Signalstrukturen auftreten (vgl. Abbildung 2.10) [57]. 
Um komplexe Strukturen bzw. einzelnen Signalkomponenten akkurat wieder- 
zugeben, wird die db20-Wavelet verwendet. 


2.8 Fazit Sensorkonzept 


Es wurde das theoretische Wirkprinzip zur Erfassung der ferromagnetischen 
Inhomogenitäten und die Funktionsweise der Messschaltung beschrieben. Auf 
Basis der Simulationsergebnisse wurde ein funktionsfähiger Prototyp des DIS 
aufgebaut. Die Sensorplattform besteht aus zwei Detektoren A und B, den 
dazugehörigen Messschaltungen und einer MCU. Die Ergebnisse der realen 
Messschaltung stimmen mit den Ergebnissen aus den Simulationen überein. 
Aus der vorgegebenen Arbeitsfrequenz fo = 30 KHz konnten die äquidistanten 
zeitlichen Abstände von 33,3 us für jede Messung durch die MCU bestätigt 
werden. Die Unterschiede in den Laufzeitdifferenzen an den beiden Detekto- 
ren sind auf Fertigungstoleranzen der elektronischen Bauteile und die unter- 
schiedlichen Positionen entlang der Schiene zurückzuführen. Die ermittelte 
Standardabweichung der realen Messschaltung liegt für beide Detektoren unter- 
halb der theoretischen Auflösung der Messschaltung (vgl. Gleichung 2.4). Zur 
Unterdrückung des Quantisierungsrauschens in den Messschaltungen wird ein 
Wavelet-basiertes Filterverfahren auf die Messdaten des DIS angewandt. Mit 
einem soft-thresholding der Detailkoeffizienten werden definierte Frequenz- 
bänder aus dem Messsignal herausgefiltert, ohne markante Signalstrukturen in 
den Messsignalen zu verfälschen. Gleichzeitig liefert das Wavelet-basiertes Fil- 
terverfahren über die Approximationskoeflizienten ein effektives Verfahren um 
die Anzahl der zu speichernden Datenpunkte zu reduzieren. Das vorgestellte 
Sensorkonzept und der DIS erfüllen alle technischen Voraussetzungen, um eine 
bordautonome Erfassung der Geschwindigkeit und der absoluten Position über 
die ferromagnetischen Inhomogenitäten zu tätigen. 
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Für die Untersuchung der berührungslosen und bordautonomen Geschwin- 
digkeitsmessung kommt ausschließlich der DIS auf dem Schienenfahrzeug 
zum Einsatz. Das bedeutet, es stehen keine zusätzlichen Sensorinformatio- 
nen im Schienenfahrzeug oder entlang der Gleisstrecke zu Verfügung, um 
die Geschwindigkeit oder die absolute Position des Schienenfahrzeugs zu 
ermitteln. Für eine präzise Geschwindigkeitsmessung ist deshalb ein Mess- 
verfahren erforderlich, welches robust gegenüber äußeren Umwelteinflüssen 
(Schnee, Regen, Spritzwasser, Öl und Staub), verschleißfrei und unempfind- 
lich gegenüber Schlupf zwischen Rad und Schiene ist. Um diesen Herausfor- 
derungen zu begegnen, ist der DIS mit zwei baugleichen Detektoren (siehe 
Kapitel 2.3) ausgestattet, die längs zur Fahrtrichtung des Schienenfahrzeugs 
in einem festen Abstand zueinander angebracht sind (siehe Abbildung 2.7 
und 3.1). Mit dem in Kapitel 2.3 beschriebenen elektromagnetischen Messver- 
fahren zur Bestimmung der ferromagnetischen Inhomogenitäten, wird über ein 
Template-Matching-Verfahren berührungslos die Geschwindigkeit des Schie- 
nenfahrzeugs bestimmt. 

Template-Matching- bzw. Laufzeit-Korrelations-Verfahren besitzen ein breites 
Anwendungsspektrum, u. a. zur Messung der Entfernung bei Radaren [58], 
zum Auffinden von korrespondierenden Bildelementen in der digitalen Bild- 
verarbeitung [59] oder zur berührungslosen Geschwindigkeitsbestimmung von 
Oberflächen [60]. Mit dem Eddy Current Sensor (ECS) wurde das Laufzeit- 
Korrelations-Verfahren bereits angewandt, um berührungslos die Geschwin- 
digkeit von Schienenfahrzeugen zu bestimmen [16], [17]. Allerdings treten 
beim ECS aufgrund der periodischen Sensorsignale Mehrdeutigkeiten bei der 
Signalzuordnung auf, wodurch eine eindeutige Bestimmung der Geschwindig- 
keit erschwert wird. 
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Abbildung 3.1: Seitliche Montageansicht des DIS an einem Schienenfahrzeug in der Nähe des 
Drehgestells zusammen mit einer schematischen Darstellung der ferromagnetischen Inhomogeni- 
täten entlang der Schiene. 


3.1 Prinzip der Geschwindigkeitsmessung 


Die an den beiden Detektoren des DIS gemessenen ferromagnetischen Inhomo- 
genitäten werden in ein Image- und ein Template-Signal aufgeteilt. Das Image- 
Signal dient als Referenz und beinhaltet die zuletzt gemessenen ferromagneti- 
schen Inhomogenitäten der Schiene. Das Template-Signal beinhaltet ebenfalls 
die zuletzt gemessenen ferromagnetischen Inhomogenitäten der Schiene, jedoch 
an dem zweiten Detektor. Bei dem Template-Matching werden die Signalver- 
läufe der beiden Detektoren miteinander verglichen und der zeitliche (oder 
räumliche) Versatz der Signale bestimmt. Dabei wird das kürzere Template- 
Signal als Schablone verwendet und in dem Image-Signal gesucht. Die Signal- 
längen des Image- und Template-Signals können dabei je nach Anwendungsfall 
und Zustand des Schienenfahrzeugs (Geschwindigkeit oder Rechenleistung des 
Auswerterechners) variieren. Dabei sind vor allem die letzten Sekunden der 
beiden Detektorsignale für eine zeitnahe Geschwindigkeitsbestimmung von 
Bedeutung. In der Abbildung 3.2 (oben) sind die Wavelet-gefilterten und Min- 
Max-normierten Signalverläufe (im Intervall [-1, 1]) der Detektoren A (rot, 
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Abbildung 3.2: Oben: Signalverlauf der beiden Detektoren A (Image / 4 s / rot) und B (Template 
/ 2 s / blau) des DIS. Unten links: RC-Signalverlauf der Korrelationskoeffizienten, quadratische 
Polynome fj und fy der lokalen Maxima und f; des globalen Maximums. Unten rechts: Detailansicht 
des quadratischen Polynoms f}, an dem eine Übereinstimmung von RC = 0,98 erreicht wurde. Es 
ergibt sich ein zeitlicher Versatz von tT = 0,12 s. 


Image-Signal) und B (blau, Template-Signal) der letzten vier bzw. zwei Sekun- 
den! des DIS zu sehen. Bedingt durch den festen Abstand der Detektoren 
zueinander und der Montageposition des DIS an dem Drehgestell des Schie- 
nenfahrzeugs (siehe Abbildung 3.1) unterscheiden sich die Signalverläufe der 
beiden Detektoren, bei konstanter Geschwindigkeit des Schienenfahrzeugs, nur 
durch einen zeitlichen Versatz r voneinander. Bereits beim ECS wurde ein 
Korrelationsverfahren verwendet, um den zeitlichen Signalversatz der Signale 
zu ermitteln [16], [17]. Aufgrund der periodischen Signalstruktur des ECS, 
vgl. Kapitel 1.4, besitzt die Korrelationsfunktion mehrere nahezu identische 
Maxima an unterschiedlichen Positionen [17]. Eine eindeutiger Signalversatz 


l Die zeitlichen Abstufungen können abhängig von der Geschwindigkeit und der Güte des Signales 
variieren. 
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ist aufgrund der Mehrdeutigkeiten nicht gegeben. Beim ECS wird dabei stets 
das erste Maxima aus der Korrelationsfunktion verwendet. Im Vergleich zum 
ECS besitzt die Korrelationsfunktion RC des DIS durch die Messungen der 
einzigartigen ferromagnetischen Signaturen der Schiene ein eindeutiges Maxi- 
mum, vgl. Abbildung 3.2. Bedingt durch die Signallängen, den Abstand der 
Detektoren zueinander und die Trägheit des Schienenfahrzeugs sind Streckun- 
gen und Stauchungen zwischen den Signale aufgrund von Beschleunigungen 
oder Abbremsvorgängen vernachlässigbar. Die Ähnlichkeit der Signalverläufe 
wird mit einem Zero Mean Normalized Cross Correlation (ZNCC) [59], [61] 
Template-Matching-Verfahren bestimmt. Dieses ergibt sich aus: 


STR]: IL + x’) 


Vas Te DT (Lx + x’) 


RC[x,T, I] = 


-19 ee 3.1) 


mit der diskreten Verschiebung x der Signale zueinander, der Elementlänge des 
Template Sr, den mittelwertfreien Signalkomponenten 7, T und den dazuge- 
hörigen Standardabweichungen oz, or. Der Korrelationskoeffizient RC kann 
Werte aus dem Intervall [-1,1] annehmen und gibt den Grad des linearen Zusam- 
menhangs der Detektorsignale an. Der Schwerpunkt des ZNCC liegt auf dem 
Vergleich der Signalstruktur. Das Messverfahren ist dadurch robust gegenüber 
linearen Verschiebungen in den Messsignalen, die bspw. auf Fertigungstoleran- 
zen der elektrischen Bauteile zurückzuführen sind. Dadurch entfällt die Not- 
wendigkeit, die Detektoren zu kalibrieren oder gegeneinander abzugleichen. 
Für den Fall, dass RC = 1 ist, unterscheiden sich die Detektorsignale A und B 
höchstens durch einen Skalierungsfaktor oder Offset voneinander. Kleine bzw. 
negative RC-Koeffizienten repräsentieren einen geringen bzw. antisymmetri- 
schen (negativen) Zusammenhang der Detektorsignale und damit den Bereich 
in dem keine Übereinstimmung zwischen den Signalen vorliegt. Durch das 
Messverfahren ist der Versatz x der RC-Funktion (Gleichung 3.1) auf diskrete 
Werte {0,1,2,...} beschränkt und limitiert dadurch die spätere Geschwin- 
digkeitsauflösung. Um dies zu umgehen, werden die RC-Koeffizienten in der 
Umgebung der jeweiligen Maxima i mit einem quadratischen Polynom f;(x) 
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angefittet. Die neue Position des Maxima (siehe Abbildung 3.2 unten links und 
rechts) wird über die Gleichung: 


Xmax 


arg max(f;(x)) 
arg max(RC[x, T, I]) (3.2) 


Q 


bestimmt?. Zur robusten Prädiktion des zeitlichen Versatzes werden beim DIS 
zusätzlich die Nebenmaxima, X max, ;, im RC-Signalverlauf (siehe Abbildung 3.2 
fı und f2) über mehrere Zeitschritte hinweg geclustert. Die Cluster werden so 
gewählt, dass diejenigen Merkmals-Vektoren (Xmax,i, fi(Xmax,i))r, die einander 
am ähnlichsten sind (Versatz und RC-Koeffizienten über mehrere Zeitschritte t 
hinweg), einer gemeinsamen Gruppe zugeordnet werden. Dabei variiert die 
Anzahl der Maxima in der RC-Funktion aufgrund der komplexen ferroma- 
gnetischen Signalstrukturen zwischen den Zeitschritten. Zur Gruppierung der 
unstrukturierten Merkmals-Vektoren wird das Dirichlet Process Mixture Model 
(DPMM) [62] Clusterverfahren gewählt. Das DPMM verwendet im Gegensatz 
zu den herkömmlichen Clusteranalysemethoden keine vorab definierte Clus- 
tergröße, sondern bestimmt diese aus den Merkmals-Vektoren. In der Abbil- 
dung 3.3 sind für 40 Zeitschritte (à 0,2 s) die RC-Korrelationsfunktionen mit 
den dazugehörigen Maxima (für RC-Koeffizienten > 0) zu sehen. Das DPMM- 
Clustering unterteilt die Merkmals-Vektoren in 7 Cluster, die sich aus den 
Nebenmaxima C1, C3, C4, C5, C6 und C7 sowie dem Hauptmaximum C2 
zusammensetzen und dadurch die Ergebnisse aus den vorherigen Zeitschritten 
beriicksichtigen. Fehldetektionen oder Mehrdeutigkeiten aufgrund von zeitlich 
begrenzten Signalschwankungen (siehe schwarze Kästen in der Abbildung 3.3) 
werden über eine Prädiktion der Merkmals-Vektoren unterbunden. Liegt kei- 
ner der aktuellen Merkmals-Vektoren in dem vorhergesagten Cluster bzw. ist 
die Vorhersagewahrscheinlichkeit nicht aussagekräftig genug, wird der vorhe- 
rige Wert verwendet und die aktuelle Messung verworfen. Das Cluster C2 
(roter Kasten) in der Abbildung 3.3 an der Position ux = 2808,49 besitzt einen 
durchschnittlichen Korrelationswert von urc = 0,97 und dadurch die höchste 


? Im Vergleich zum DIS besitzt die Korrelationsfunktion des ECS mehrere nahezu identische 
Maxima an unterschiedlichen Positionen in der Korrelationsfunktion [17]. Störungen im Messi- 
gnal können dabei zu einer falschen Signalversatz führen. Eine Bestimmung der Geschwindigkeit 
ist dann nicht mehr gegeben. 
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Abbildung 3.3: Zeitliches DPMM-Clustering der lokalen Maxima fj (Xmax,i) der Korrelations- 
funktion RC über eine Dauer von 8 s, mit der einfachen o Umgebung (inneren Ellipse) und der 
doppelten o Umgebung (äußeren Ellipse). Für eine bessere Darstellung wurde die Messfrequenz 
des DIS auf f) = 1,5 kHz reduziert. 


Aufenthaltswahrscheinlichkeit. Mit der Messfrequenz fo des DIS und den Start- 
zeiten tr und t; der Signalverläufe des Template- und Imagesignals, ergibt sich 
der zeitliche Versatz der Detektrosignale aus: 


1 
T = (tr — t1) Amen ` 7- (3.3) 
fo 
Mit dem Abstand d der Detektoren ist die Geschwindigkeit: 
d 
v=- (3.4) 
T 


des Schienenfahrzeugs gegeben. Das in der Abbildung 3.3 gezeigte Szenario 
zeigt über den Cluster C2 den Geschwindigkeitsverlauf der letzten acht Sekun- 
den an. Es ist anzumerken, dass im Gegensatz zum ECS? der Abstand der 
Detektoren beim DIS aufgrund der nicht periodischen Signalstruktur beliebig 
gewählt werden kann. Bei hohen Geschwindigkeiten profitiert insbesondere 
das Template-Matching-Verfahren von einem großen Abstand zwischen den 
Detektoren. 


3 Beim ECS ist der Abstand durch den mittleren Schwellenabstand (ca. 60 cm) gegeben. 
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Abbildung 3.4: Darstellung des Schienenprüfstandes zusammen mit einer Detailansicht des DIS. 
3.2 Versuche am Schienenprüfstand 


Die Funktionsweise des DIS wurde auf dem Schienenprüfstand unter kon- 
trollierten Laborbedingungen getestet. Der Schienenprüfstand, siehe Abbil- 
dung 3.4, ist mit einem beweglichen Führungsschlitten ausgestattet. An diesem 
wird der DIS befestigt und periodisch mit einer Maximalgeschwindigkeit von 
3 m/s über die darunterliegenden Schienenelemente geführt [34]. Mit einer 
effektiven Länge von 6,7 m werden über den DIS die ferromagnetischen Inho- 
mogenitäten der darunterliegenden Vignolschienen und eines Weichenelemen- 
tes (Herzstück) erfasst. Anhand des Schienenprüfstandes wird die Hardware 
des DIS, die Signalfilterung und die Auswertesoftware für die Geschwindig- 
keitsbestimmung unter kontrollierten Laborbedingungen untersucht und mit 
der Referenzgeschwindigkeit des Schienenprüfstandes verglichen. Der Aufbau 
ermöglicht es, unterschiedliche Montagepositionen über dem Schienenkopf und 
quer zur Fahrtrichtung miteinander zu vergleichen. 
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Abbildung 3.5: Signalverlauf der Detektoren A (rot) und B (blau) auf dem Schienenprüfstand 
mit den Wendepunkten des Führungsschlittens (schwarze gestrichelten Linien). Die ferromagneti- 
schen Strukturen des Detektors B (gelber Signalbereich) sind vollständig in den ferromagnetischen 
Strukturen des Detektors A (magentafarbener Signalbereich) enthalten [37]. 


3.2.1 Fahrtrichtungserkennung 


In der Abbildung 3.5 sind für die beiden Detektoren des DIS die gemesse- 
nen ferromagnetischen Inhomogenitäten auf dem Schienenprüfstand zu sehen. 
Durch die kontinuierliche Messwerterfassung des DIS und den ständigen Rich- 
tungswechsel des Führungsschlittens auf dem Schienenprüfstand tritt an den 
Umkehrpunkten ein Phasensprung zwischen den beiden Signalverläufen auf. 
Kommt es zu einem Richtungswechsel (schwarz gestrichelte Linien), muss 
die Zuordnung der Detektoren zu den Signalkomponenten T und / für das 
Template-Matching (Gleichung 3.1) getauscht werden. Anderenfalls beinhaltet 
das gesuchte Template-Signal Merkmalspunkte der Schiene, die in dem Image- 
Signal noch nicht gemessen wurden. Ein Signalmatching bzw. die Geschwin- 
digkeitsbestimmung ist in diesem Fall nur noch bedingt möglich. Die Erkennung 
der Fahrtrichtung und damit die Zuordnung der Signale zu den Detektoren wird 
bei der Initialisierung des DIS (beim Anfahren) durchgeführt. Zunächst wird 
eines der Detektorsignale willkürlich als Image-Signal gewählt. Die Länge der 
ferromagnetischen Signatur ist so gewählt, dass die Signalkomponente des zwei- 
ten Detektorsignals (Template-Signal) vollständig vom ersten eingeschlossen 
ist (vgl. Abbildung 3.5 gelber und magentafarbener Signalbereich). Unabhängig 
von der Fahrtrichtung und der Zuordnung der Detektoren sind die ferromagne- 
tischen Signalstrukturen des Template-Signals dadurch immer vollständig im 
Image-Signal enthalten. Aus dem zeitlichen Versatz (Gleichung 3.3) ergibt sich 
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Abbildung 3.6: Querschnitt durch eine Vignolschiene. Dargestellt sind die unterschiedlichen Mon- 
tagepositionen des DIS über dem Schienenkopf [37]. 


die Bewegungsrichtung des Schienenfahrzeugs in Abhängigkeit des Vorzei- 
chens: 


—1 negative Fahrtrichtung 


sgn(T) = (3.5) 


1 positive Fahrtrichtung ` 


Über die Fahrtrichtung sind die Signalkomponenten des Template- und des 
Image-Signals eindeutig den Detektoren zuzuordnen. Dies ist entscheidend, 
um bei niedrigen Geschwindigkeiten die Position des Maximums exakt mit der 
RC-Korrelationsfunktion (siehe Abbildung 3.3) bestimmen zu können. 
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Montageposition 1 Montageposition 2 Montageposition 3 
1.2 1.2 1.2 
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£ 1.0 Be ee en & 1.0 = Seele — SE 1.0 geen eee = 
=09 =09 =09 ui eem 
0.8 0.8 0.8 
2 3 4 1 2 3 1 2 
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Montageposition 4 Montageposition 5 Montageposition 6 
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1.2 T 1.2 1.2 = 
E08 j . £0. 
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0.2 
1 2 3 4 
Durchlauf Durchlauf Durchlauf 


Abbildung 3.7: Betrag der mittleren Geschwindigkeit Joore! mit der jeweiligen Standardabwei- 
chung o, des DIS für die Montagepositionen 1 bis 6 an dem Schienenprüfstand für vier Durchläufe. 
Rote gestrichelte Linie: Referenzgeschwindigkeit vp = 1 m/s an dem Schienenprüfstand. 


3.2.2 Einfluss der Montageposition auf die 
Geschwindigkeitsbestimmung 


Die ferromagnetischen Inhomogenitäten entlang der Schiene dienen als Land- 
marken für die Geschwindigkeitsbestimmung des Schienenfahrzeugs. Für eine 
optimale Erfassung der ferromagnetischen Inhomogenitäten ist eine Monta- 
geposition des DIS in der Nähe des Drehgestells an dem Schienenfahrzeug 
vorgesehen. Während einer Kurvenfahrt wird dadurch sichergestellt, dass beide 
Detektoren dieselben ferromagnetischen Strukturen erfassen. Zudem wird eine 
zu große Auslenkung der Detektoren quer zur Schiene verhindert. Um die Aus- 
wirkungen der Montagepositionen auf die Signalqualität an den Detektoren 
zu untersuchen, werden für verschiedene Montageposition (vgl. Abbildung 3.6 
und A.6) auf dem Schienenprüfstand Geschwindigkeitsmessungen mit dem DIS 
durchgeführt. In der Abbildung 3.7 sind die Beträge der mittleren Geschwin- 
digkeiten pts für die Montagepositionen (1 bis 6) zusammen mit den Standard- 
abweichungen o für jeweils vier Fahrten zu sehen. Die vom Schienenprüfstand 
vorgegebene Geschwindigkeit liegt für alle Messungen innerhalb der ermit- 
telten Standardabweichung des DIS. Für eine asymmetrische Anordnung der 
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Geschwindigkeit v / (m/s) 


Zeit t / (s) 


Abbildung 3.8: Gemessenes Geschwindigkeitsprofil des DIS an dem Schienenprüfstand 10 cm über 
dem Schienenkopf. Eine detaillierte Übersicht über die eingestellten Referenzgeschwindigkeiten 
des Schienenprüfstandes sind in der Tabelle 3.1 zu sehen. 


Detektoren über dem Schienenkopf (Sensorposition 1 und 3 bzw. Sensorpo- 
sition 4 und 6) ist zu erkennen, dass die Messunsicherheit mit zunehmender 
Entfernung zum Schienenkopf steigt. Die Montagepositionen (2 und 5) direkt 
über dem Schienenkopf liefern mit einer Standardabweichung von a < 0,2 m/s 
bzw. 0,1 m/s die beste Übereinstimmung zu der Referenzgeschwindigkeit von 
1,0 m/s des Schienenprüfstandes. Alle weiteren Versuchsfahrten mit dem DIS 
werden mit einem Abstand von ca. 10 cm über dem Schienenkopf durchgeführt. 
Dies entspricht einer Position, die zwischen der Montageposition 2 und 5 liegt. 
Dadurch wird das vorgeschriebene Lichtraumprofil und der einzuhaltenden 
Sicherheitsabstand von 10 cm über dem Schienenkopf eingehalten [63]. 


3.2.3 Geschwindigkeitsmessungen auf dem 
Schienenprüfstand 


Das Gesamtsystem DIS - inklusive der Erfassung der ferromagnetischen 
Inhomogenitäten, Signalfilterung, Richtungserkennung und Bestimmung der 
Geschwindigkeit — wurde auf dem Schienenprüfstand, 10 cm über dem Schie- 
nenkopf, evaluiert. Das gemessene Geschwindigkeitsprofil des DIS ist in der 
Abbildung 3.8 zu sehen. Dabei wurden im Niedriggeschwindigkeitsbereich acht 
Richtungswechsel (a bis h) durchlaufen. In der Tabelle 3.1 sind die mittleren 
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Tabelle 3.1: Mittlerer Geschwindigkeitsverlauf vs des DIS auf dem Schienenprüfstand für die 
Durchläufe a bis h in einem Abstand von 10 cm über dem Schienenkopf. Mit der Standard- 
abweichung œ, der eingestellten Referenzgeschwindigkeit vp des Schienenprüfstandes und der 
Geschwindigkeitsabweichung |Av| [38]. 


a b c d € f g h 


pos / (m/s) 1,09 -1,21 127 -1,28 147 -1,46 1,44 -1,60 
o /(m/s) 0,12 0,13 0,11 010 0,06 0,12 012 0,19 
vp /(m/s) 1,00 -1,20 1,20 -1,20 1,40 -1,40 -1,60 1,60 
|Av| / (m/s) 0,09 0,01 0,07 0,08 0,07 0,06 0,16 0,00 


Geschwindigkeiten des DIS vpıs zusammen mit den Referenzgeschwindig- 
keiten des Schienenprüfstandes vp für die einzelnen Zeitabschnitte hinterlegt. 
Für die Durchläufe a bis f und h liegt die Referenzgeschwindigkeit innerhalb 
der mit dem DIS ermittelten Standardabweichung bei einer durchschnittlichen 
Geschwidigkeitsabweichung von 0,05 m/s. Bei dem Durchlauf g ergibt sich 
die höchste Geschwindigkeitsabweichung des DIS mit einer Abweichung von 
Av =0,16 m/s zu der eingestellten Geschwindigkeit auf dem Schienenprüfstand. 
Die erhöhte Abweichung ist auf abrupte Richtungswechsel an den Umkehrpunk- 
ten des Schienenprüfstandes und fehlende Signalvariabilität an den Detektoren 
zurückzuführen®. 


3.2.4 Fazit Schienenprüfstand 


Den Kern der Untersuchungen am Schienenprüfstand bildete das berührungs- 
lose Messverfahren zur Erfassen der ferromagnetischen Inhomogenitäten ent- 
lang der Schiene, die Signalfilterung (vgl. Kapitel 2.7), die Richtungserkennung 
und die Bestimmung der Geschwindigkeit. Darüber hinaus wurden verschie- 
dene Montagepositionen und deren Einfluss auf die Geschwindigkeitsmessung 
untersucht. Unabhängig von der Montageposition des DIS ergaben sich für alle 
Fahrten auf dem Schienenprüfstand im anspruchsvollen Niedriggeschwindig- 
keitsbereich (< 2 m/s) Abweichungen, die unterhalb der geforderten 0,56 m/s 
liegen und damit den Anforderungen moderner Zugbeeinflussungssysteme, wie 


4 Umkehrpunkte des Schienenprüfstandes repräsentiert ein Szenario, welches im regulären Betrieb 
nicht vorgesehen ist. 
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Schiene 
mittlere Schiene 
@ linke Schiene 


N-Distanz / (m) 


im mn — EE 


u a i m 
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Abbildung 3.9: Satellitenansicht des T4-Rings (Map data ©2016 Google Z:32U; E:306076; 
N:5666477). Die Gleisstrecke ist in drei Schienenabschnitte, die rechte, mittlere und linke Schiene, 
unterteilt. Heatmap: Digitales Höhenmodell und Schienenvektorisierung zusammen mit den Kol- 
limationsmarkierungen (schwarze Kreise) [64]. 


dem ETCS, für niedrige Geschwindigkeiten (< 30 km/h) genügen [7]. Unter 
Berücksichtigung der sicherheitstechnischen Vorgaben wurde eine ideale Mon- 
tageposition des DIS zur Erfassung der ferromagnetischen Inhomogenitäten in 
der Nähe des Drehgestells über dem Schienenkopf ermittelt. Dadurch ist eine 
minimale seitliche Auslenkung der Detektoren bei Kurvenfahrten gewährleistet. 


3.3 Feldversuche 


Ausgehend von den Ergebnissen auf dem Schienenprüfstand wurden Untersu- 
chungen mit dem DIS auf einem realen Schienenfahrzeug auf der Testanlage 
des PCW? durchgeführt. Die Versuchsfahrten fanden bei starkem Regen auf 
der Teststrecke T4 statt. Dieser Streckenabschnitt besitzt ein redundantes und 
diversitär aufgebautes Messsystem zur absoluten Positionsbestimmung, über 
welches ebenfalls die Geschwindigkeit des Schienenfahrzeugs ermittelt wird. 


5 Prüf- und Validationcenter Wegberg-Wildenrath (PCW). 
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Die Auflösung und insbesondere die Genauigkeit der Referenzsensorik, die auf 
dem T4 zu Verfügung steht, liegt dabei weit über dem Standard kommerzieller 
Sensoren, die im Schienenverkehr zum Einsatz kommen. Das georeferenzierte 
Testgelände besitzt dadurch ideale Voraussetzungen, um die Leistungsfähigkeit 
des DIS gefahrlos unter realen Bedingungen erproben zu können. 


3.3.1 Versuchsaufbau und Referenzsensorik 


Für die Erfassung der ferromagnetischen Strukturen (siehe Kapitel 2.4) wird die 
ermittelte ideale Montageposition aus Kapitel 3.2.2 für den DIS (in der Nähe 
der Vorderachse des Erprobungsfahrzeugs® und einem Abstand von ca. 10 cm 
über dem Schienenkopf) verwendet (siehe Abbildung 2.7). Das Erprobungsfahr- 
zeug verfügt über einen elektrischen Antrieb und wird mit einer Fernsteuerung 
über die Gleise bewegt. Als Teststrecke wurde das Testgleis T47 (siehe Abbil- 
dung 3.9), welches sich innerhalb des railGATE [65] befindet, verwendet. Der 
T4 besitzt ein eingemessenes, hoch präzises geodätisches Netz, welches aus 
sechs georeferenzierten Messstationen (Beobachtungspfeiler) besteht. Zusam- 
men mit dem railGATE bilden diese die Referenz für die Geschwindigkeit 
und die Lokalisierung des Erprobungsfahrzeugs. Das Gesamtsystem wurde 
mit einer Genauigkeit im Submillimeterbereich [66], [64] eingemessen. Im 
Nachfolgenden werden die Referenzsensoren des Testfeldes bzw. des Erpro- 
bungsfahrzeugs vorgestellt: 


e Rundprisma / Leica Viva TS15 (TS15): 

Das Leica Viva TS15 ist ein LASER-Tachymeter der Schweizer Firma 
Leica Geosystems AG [67]. Über eine von sechs georeferenzierten Mess- 
stationen (Beobachtungspfeiler) wird mit Hilfe eines Lasers der Winkel 
und die Entfernung zu einem Rundprisma ermittelt [66]. Mit den Beob- 
achtungspfeilern und den Messpunkten des TS15 wird das Rundprisma, 
welches sich zusammen mit der GNSS-Antenne auf dem Erprobungs- 
fahrzeug befindet, in ein kartesisches Koordinatensystem (UTM32) über- 
führt. 


6 Gang car (Meterspur) der Firma ARIC / Automotive and Rail Innovation Center. 
7 Länge 553 m / Radius 50 m. 
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e railGATE / Novatel GPS 703 GGG (GNSS): 

Bei dem railGATE handelt es sich um eine Versuchsanlage zur Untersu- 
chung, Erprobung und Beurteilung der nächsten Generation von GNSS- 
Receivern und deren Anwendung im Bahnbereich. Das railGATE besteht 
aus acht, 40 bis 60 Meter hohen Türmen, den sogenannten Galileo Pseu- 
doliten (terrestrischer Sender), die als künstliche Satelliten fungieren und 
Galileo-konforme Signale bereitstellen. Die Versuchsanlage des PCW 
ist dadurch bereits in der Lage, das Satellitennavigationssystem Gali- 
leo zu nutzen. Aus den Standorten der Pseudoliten und der Laufzeit- 
messung der einzelnen Signale ergibt sich die Position der GPS-703- 
GGG Antenne [68]. Das railGATE besitzt eine Positionsgenauigkeit von 
1,3 m (Base-Mode) bzw. 0,3 m im Assisted Mode [69]. 


Das Rundprisma und die GNSS-Antenne befinden sich auf derselben verti- 
kalen Achse und sind gemeinsam an dem Erprobungsfahrzeug montiert. Das 
Gesamtsystem besitzt eine absolute Genauigkeit von +2 cm und eine relative 
Genauigkeit von +2 mm zu den Beobachtungspfeilern und damit eine höhere 
Auflösung als herkömmliche Sensoren, die für die Positionsbestimmung von 
Schienenfahrzeugen in diesem Bereich zum Einsatz kommen [66]. Aus dem 
zeitlichen Abstand von zwei aufeinanderfolgenden Koordinatenmesspunkten 
wird die Geschwindigkeit des Erprobungsfahrzeugs ermittelt. 


3.3.2 Versuchsfahrten 


Aus der unregelmäßigen Verteilung der ferromagnetischen Strukturen, die sich 
entlang der Schiene befinden, sind für das eingesetzte Mess- und Auswer- 
teverfahren zur Geschwindigkeitsbestimmung vor allem diejenigen Signalab- 
schnitte von besonderem Interesse, bei denen eine geringe Signalvariabilität an 
den Detektoren gemessen wird. Dies ist insbesondere im Niedriggeschwindig- 
keitsbereich der Fall, weshalb es ein besonders anspruchsvolles Testszenario 
darstellt. Um die Leistungsfähigkeit des DIS speziell unter diesen erschwer- 
ten Bedingungen zu untersuchen, wurden auf der Teststrecke unterschiedliche 
Anfahrvorgänge (0 m/s bis x 1 m/s) und Rangierfahrten im Niedriggeschwin- 
digkeitsbereich (~ 1 m/s bis 4 m/s) durchgeführt. Dazwischen befinden sich 
Abschnitte, in denen mit konstanter Geschwindigkeit über die Teststrecke gefah- 
ren wurde. Die maximale Geschwindigkeit betrug dabei 8 m/s. Eine detaillierte 
Ansicht der einzelnen Streckenabschnitte (ID 1 bis 33) mit den dazugehörigen 
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Geschwindigkeitsprofilen ist in der Abbildung 3.10 zu sehen. Der Parameter O 
(0 rechte Schiene und 1 linke Schiene) kennzeichnet die Montageposition des 
DIS über der jeweiligen Schiene. Durch das Rangieren auf dem T4 wechselt die 
Schiene, über der sich der DIS befindet. Beide Schienen eines Streckenabschnit- 
tes besitzen eine individuelle ferromagnetische Signatur. In der Abbildung 3.9 
ist dies durch die grüne und schwarze Kurve gekennzeichnet. 


N-Distanz d / (m) 


Geschwindigkeit v / (m/s) 


50 100 150 200 250 0 50 100 150 200 250 0 50 100 150 200 250 


0 50 100 150 200 250 


E-Distanz d / (m) 


Abbildung 3.10: Position und Geschwindigkeit des Erprobungsfahrzeugs auf der Teststrecke T4 
(Koordinatenoffset Z:32U; E:306026; N:5666467) unterteilt in zusammenhängende und nicht über- 
lappende Sequenzen mit einer eindeutigen ID. O ist die Position des DIS tiber der rechten (0) bzw. 
linken (1) Schiene eines Gleisabschnittes. 
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Abbildung 3.11: (a) Detailansicht der ferromagnetischen Inhomogenitäten für das Fahrmanöver 
ID = 3. Aus den ferromagnetischen Inhomogenitäten werden die Cluster (b) und die Geschwin- 


digkeit (c) des Erprobungsfahrzeugs ermittelt. (d) zeigt die Geschwindigkeitsabweichung zu der 
Referenzsensorik. 
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3.3.3 Geschwindigkeitsmessungen im Feldversuch 


Für die Geschwindigkeitsuntersuchungen des DIS werden die oben beschrie- 
benen Fahrmanöver auf der Teststrecke T4 des PCW durchgeführt und mit 
den Messungen der Referenzsensoren verglichen. Eine detaillierte Ansicht der 
gemessenen ferromagnetischen Strukturen des Fahrmanöver mit der ID = 3 
(siehe Abbildung 3.10) ist in der Abbildung 3.11 a) zu sehen. Die Signal- 
komponenten des Image- und Templatesignals besitzen für den Vorgang der 
Geschwindigkeitsmessung jeweils eine zeitliche Länge von 7,0 s und 4,5 s. Die 
einzelnen Maxima der RC-Funktionen werden mit einem DPMM-Clustering 
abschnittsweise in die Cluster {C1, C2, C3, ...} aufgeteilt und für die gesamte 
Dauer des Fahrmanövers miteinander verknüpft (siehe b). Dadurch werden 
Mehrdeutigkeiten im Niedriggeschwindigkeitsbereich (siehe t = 174 s), die 
durch fehlenden Signalvariabilität auftreten, kompensiert. Das Cluster C1, mit 
einem mittleren RC-Koeffizienten von 0,96 und einer Standardabweichung von 
0,07, entspricht dem wahrscheinlichsten zeitlichen Versatz für das gesamte Fahr- 
manöver. Aus dem Abstand der Detektoren (d = 0,6 m) ergibt sich für jeden Zeit- 
schritt die Geschwindigkeit des Erprobungsfahrzeugs über die Gleichung 3.4. 
In der Abbildung 3.11 (c) ist der dazugehörige zeitliche Verlauf zusammen mit 
der Referenzgeschwindigkeit zu sehen. Für Geschwindigkeiten kleiner 0,5 m/s 
(t=0, 158, 192 s) ist ein Rückgang der Signalvariabilität an den Detektoren und 
ein Anstieg der Geschwindigkeitsabweichung zur Referenzgeschwindigkeit zu 
beobachten. Für das gesamte Fahrmanöver ergibt sich eine mittlere Geschwin- 
digkeitsabweichung des DIS zur Referenzgeschwindigkeit von Av = 0 m/s 
und eine Standardabweichung von a = 0,12 m/s, siehe Abbildung 3.11 (d). 
Für das gesamte Fahrmanöver liegt die Geschwindigkeitsabweichung innerhalb 
der geforderten minimal bzw. maximal erlaubten Geschwindigkeitsabweichung 
(+0,56 m/s), die für das ETCS vorgesehen ist [7]. In der Abbildung 3.12 ist ein 
weiteres Rangier-Szenario, bestehend aus den Fahrmanövern ID = 25 bis 33 
(vgl. Abbildung 3.10), zu sehen. Bei diesen Fahrmanövern handelt es sich 
um mehrere aufeinanderfolgende Fahrtrichtungswechsel, in denen das Schie- 
nenfahrzeug beschleunigt, abbremst und anschließend in die Gegenrichtung 
erneut beschleunigt wurde. Die einzelnen Richtungswechsel gehen dabei naht- 
los ineinander über, wodurch das Erprobungsfahrzeug nur für eine sehr kurze 
Zeit zum Stehen kommt. Die Richtungswechsel sind an den Positionen t = 23, 
71, 110, 161, 197, 270, 304, 337 s, insbesondere durch die Signalverläufe an 
den Detektoren (a) und dem Cluster C1 (b), zu erkennen. Das Image- und 
Templatesignal besitzen auch hier eine zeitliche Länge von 7,0 s und 4,5 s. 
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Abbildung 3.12: (a) Detailansicht der ferromagnetischen Inhomogenitäten für die Fahrmanöver 
ID = 25 bis 33. Aus den ferromagnetischen Inhomogenitäten werden die Cluster (b) und die 
Geschwindigkeit (c) des Erprobungsfahrzeugs ermittelt. (d) zeigt die Geschwindigkeitsabweichung 
zu der Referenzsensorik. 
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Tabelle 3.2: Mittlere Geschwindigkeitsabweichung Av zu den Referenzsensoren GNSS und TS15 
sowie die Standardabweichung a. Untersucht wurden jeweils die Fahrmanöver ID = 1 bis 33 für 
unterschiedliche Templatelängen. * Bei den Übereinstimmungen handelt es sich um die prozen- 
tualen Anteile der Geschwindigkeiten des DIS, die eine Abweichung von maximal + 0,56 m/s 
aufweisen und damit die Anforderungen des ETCS erfüllen. * Ist die Templatelänge mit den besten 
Geschwindigkeitsergebnissen. 


Templatelänge / (s) 
1 2 3 EN 4,5 5 
GNSS -0,06 -0,10 zs 0,02 0,01 0,01 
Av / (m/s) 
TS15 -0,07 =x 0,03 z0 0,02 0,06 
GNSS 0,37 045 0,30 0,19 023 0,24 
a / (m/s) 
TS15 0,79 0,5 0,56 0,38 0,40 0,60 
Uberein- GNSS 89 88 95 97 97 96 


stimmungen* /(%) TS15 72 90 93 93 93 92 


Durch den häufigen Richtungswechsel liegt der mittlere RC-Koeffizient des 
C1-Clusters bei 0,90 mit einer Standardabweichung von 0,16. Im Vergleich zu 
der verwendeten Referenzsensorik liefert der DIS eine Richtung der ermittelten 
Geschwindigkeit des Erprobungsfahrzeugs und erlaubt somit, die Fahrtrichtung 
aus den ferromagnetischen Signalen zu extrahieren. In der Abbildung 3.12 (c) 
sind zur besseren Darstellung beide Geschwindigkeitsgrößen des DIS (DIS und 
|DIS|) und die Referenzgeschwindigkeit abgebildet. Die mittlere Geschwindig- 
keitsabweichung Av für das Fahrmanöver ID = 25 bis 33 beträgt dabei 0,05 m/s 
mit einer Standardabweichung von 0,23 m/s (vgl. Tabelle 3.3). In der Abbil- 
dung 3.13 links ist die Häufigkeit der Geschwindigkeitsabweichung Avpısxer. 
über die Referenzgeschwindigkeit für die Fahrmanöver ID = 1 bis 33, für 
einen Zeitraum von ca. 40 Minuten, farbcodiert dargestellt. Auf der rechten 
Seite ist das dazugehörige Histogramm der Häufigkeit über die akkumulierten 
Geschwindigkeitsabweichungen zu sehen. In beiden Darstellungen sind zusätz- 
lich die minimal und maximal zulässigen Geschwindigkeitsabweichung für das 
ETCS, blau und rot gestrichelte Linie, dargestellt. Die vereinzelt auftretenden 
Abweichungen in den Bereichen zwischen 0 bis 1 m/s und 2,5 m/s zu der 
Referenzgeschwindigkeit resultieren aus der geringen Signalvariabilität an den 
Detektoren. Durch ein längeres Templatesignal können die Ausreißer teilweise 
kompensiert werden. Diese verursachen aber einen längeren Nachlauf bei der 
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Abbildung 3.13: Links: Häufigkeit der Geschwindigkeitsabweichung aufgetragen über die Refe- 
renzgeschwindigkeit vger, des Schienenfahrzeugs für die Streckenabschnitte ID = 1 bis 33. Rechts: 
Histogramm der Häufigkeit akkumuliert über alle Geschwindigkeiten. Die blau und rot gestri- 
chelte Linie kennzeichnen die minimale und maximale zulässige Geschwindigkeitsabweichung 
nach ETCS [7]. 


Geschwindigkeitsmessung. Um auf schnell wechselnde Geschwindigkeitsän- 
derungen reagieren zu können, ist ein kurzes Templatsignal zu bevorzugen. In 
der Tabelle 3.2 ist eine Übersicht der untersuchten Templatelängen und den 
daraus ermittelten mittleren Geschwindigkeitsabweichung Av des DIS zu den 
Referenzsensoren (GNSS und TS15) sowie die Standardabweichung o für die 
Fahrmanöver ID = 1 bis 33 zu sehen. Durch die Verwendung eines längeren 
Templates profitiert vor allem der Niedriggeschwindigkeitsbereich. Für eine 
Templatelänge von 4 Sekunden besitzen 93 % bzw. 97 % (vgl. Tabelle 3.3) 
aller gemessenen Geschwindigkeiten die über den DIS ermittelt wurden eine 
Geschwindigkeitsabweichung von weniger als 0,56 m/s und entspricht somit 
ebenfalls dem ETCS Standard [7]. 
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Tabelle 3.3: Abweichung zwischen der ermittelten Geschwindigkeit des DIS und der Referenzge- 
schwindigkeit für unterschiedliche Fahrmanöver (ID). 


ID=3 ID = 25 bis 33 
Templatelänge / (s) 4,0 4,0 
Av / (m/s) x0 0,05 
Re / (m/s) 0,12 0,23 


3.4 Zusammenfassung der Ergebnisse 


Um den Anforderungen moderner Zugkontroll- und Zugsicherungssysteme zu 
genügen und den Einschränkungen heutiger Geschwindigkeitssensoren ent- 
gegenzuwirken, wurde ein neuartiges Sensorkonzept entwickelt. Durch ein 
elektromagnetisches Messverfahren, welches für die Bestimmung der ferro- 
magnetischen Inhomogenitäten verwendet wird, kann die Geschwindigkeit des 
Schienenfahrzeugs berührungslos, d. h. ohne Verschleiß des Radumfangs oder 
Schlupf zwischen Schiene und Rad, ermittelt werden. Die Robustheit des Mess- 
verfahrens konnte anhand von Messungen auf dem PCW unter regnerischen 
Wetterbedingungen nachgewiesen werden. Das Messverfahren des DIS und die 
Geschwindigkeitsbestimmung wurden unter kontrollierten Laborbedingungen 
an einem Schienenprüfstand und auf einer Teststrecke des PCW unter rea- 
len Betriebsbedingungen evaluiert. Aus den Untersuchungen konnte eine ideale 
Montageposition in der Nähe des Drehgestells in einer Höhe von ca. 10 cm über 
dem Schienenkopf ermittelt werden. Mit einer Templatelänge von 4 Sekunden 
lagen die ermittelten Geschwindigkeiten in der Nähe der verwendeten Refe- 
renzsensorik. Die Referenzsensorik weißt dabei eine Genauigkeit auf, die weit 
über den Anforderungen heutiger Sensoren in diesem Anwendungsbereich lie- 
gen. Für alle untersuchten Szenarien ergab sich eine Geschwindigkeit, deren 
Abweichung zu der Referenzgeschwindigkeit die Anforderungen des ETCS 
im anspruchsvollen Niedriggeschiwndigkeitsbereich erfüllen. Die Ergebnisse 
zeigen deutlich das Potenzial der streckenseitigen ferromagnetischen Inhomo- 
genitäten in Verbindung mit dem DIS auf. 
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In dem vorherigen Kapitel wurde gezeigt, wie die ferromagnetischen Inhomo- 
genitäten der Schiene in Kombination mit dem DIS dazu verwendet werden, 
um die Geschwindigkeit eines Schienenfahrzeugs zu bestimmen. Wie bereits 
aus den vorangegangenen Untersuchungen bekannt, reicht die Kenntnis der 
Geschwindigkeit nicht aus, um eine gleisselektive Position des Schienenfahr- 
zeugs im Gleisnetz zu ermitteln. Hierfür sind zusätzliche Informationen, bspw. 
in Form von Merkmalen, die sich entlang der Gleisstrecke befinden, notwen- 
dig. Erst durch die Zuordnung der ferromagnetischen Inhomogenitäten an das 
jeweilige Schienenelement ergeben sich eindeutige Landmarken, die in einer 
Streckenkarte hinterlegt werden. Bei einer erneuten Befahrung der Gleisstrecke 
werden die gemessenen ferromagnetischen Inhomogenitäten mit den Signaturen 
in der Streckenkarte verglichen. Über die eindeutige Zuordnung der Landmar- 
ken an das jeweilige Schienenelement wird die absolute Position des Schienen- 
fahrzeugs ermittelt. Ein ähnliches Prinzip wird bereits im Bereich der Robotik 
und dem autonomen Fahrens angewandt. Hier werden überwiegend Bildmerk- 
male aus Kameras verwendet, um zusammen mit einer Merkmalskarte die Posi- 
tion des Fahrzeugs zu ermitteln [70]. Auch der ECS besitzt bereits ein grobes 
Lokalisierungsverfahren. Jedoch nutzt der ECS lediglich Gleisstrukturen, wie 
Weichen, zur Ermittlung der Position. Mehrdeutigkeiten in der Position kön- 
nen erst durch das Überfahren von mehreren Weichen reduziert werden [33]. 
Ein Vorteil des DIS ist die Möglichkeit, kontinuierlich anhand der ferroma- 
gnetischen Inhomogenitäten der Schiene, die Position des Schienenfahrzeugs 
zu bestimmen. Bereits nach wenigen Streckenmetern ist dadurch eine gleisse- 
lektive Position des Schienenfahrzeugs gegeben. Dadurch liefert der DIS eine 
schnellere und genauere Position des Schienenfahrzeugs. 

Im folgenden Kapitel wird das Verfahren zur Erstellung der ferromagneti- 
schen Streckenkarten erläutert. Anschließend wird mit den ferromagnetischen 
Inhomogenitäten die Position des Schienenfahrzeugs kontinuierlich entlang der 
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zurückgelegten Gleisstrecke bestimmt. Zur Untersuchung der Positionsgenau- 
igkeit werden die Positionsergebnisse aus dem DIS mit denen aus der Referenz- 
sensorik (vgl. Kapitel 3.3.1) verglichen. 


4.1 Erstellung der ferromagnetischen 
Streckenkarte 


Um die absolute Position des Schienenfahrzeugs auf der Gleisstrecke zu ermit- 
teln, werden zunächst die ferromagnetischen Inhomogenitäten der Schiene bei 
einer Kartierungsfahrt aufbereitet und in einer Streckenkarte hinterlegt. Nach- 
dem der Gleisabschnitt kartiert wurde, können die ferromagnetischen Inho- 
mogenitäten zur Lokalisierung verwendet werden. Durch die starke Bindung 
der ferromagnetischen Strukturen an die Schienenelemente wird für die Reprä- 
sentation der Karte ein Topogramm in Form eines Liniennetzplans gewählt. 
Die Generierung der Streckenkarte und die anschließende Lokalisierung des 
Schienenfahrzeugs kann unter alleiniger Verwendung des DIS durchgeführt 
werden und benötigt keine zusätzlichen Sensoren. Im Gegensatz zu etablier- 
ten Lokalisierungsverfahren im Schienen- oder Straßenverkehr werden keine 
geografischen Streckenkarten für eine gleis- bzw. spurgenaue Lokalisierung auf 
mehrspurigen Streckenabschnitten benötigt [71]. Dadurch entfallen aufwendige 
und teure Maßnahmen für die Erstellung von geeigneten georeferenzierten Stre- 
ckenkarten. Zur Erstellung der ferromagnetischen Streckenkarten werden die 
ferromagentischen Strukturen, wie bereits im Kapitel 3 gezeigt, zunächst von 
dem DIS in Abhängigkeit der Zeit erfasst. Das Messsignal besitzt dadurch eine 
geschwindigkeitsabhängige Komponente. Um den Einfluss der Geschwindig- 
keit zu kompensieren, werden die ferromagnetischen Strukturen in Abhängig- 
keit der zurückgelegten Wegstrecke aufbereitet und dem jeweiligen Strecken- 
abschnitt zugeordnet. Dies kann direkt über die ermittelte Geschwindigkeit des 
DIS (siehe Kapitel 3) oder unter Verwendung zusätzlicher Sensoren für die 
Geschwindigkeits- oder Distanzmessung erfolgen. Über die Geschwindigkeit 
des DIS wird die zurückgelegte Wegstrecke: 


Vta- + Vta 


> ) (4.1) 


N 
D:(N) = È (ln - tat) 


n=1 
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Abbildung 4.1: Darstellung der Transformation der ferromagnetischen Inhomogenitäten eines 
Detektors von der Zeit- in die Ortsdomäne. Dargestellt sind die ferromagnetischen Inhomogenitä- 
ten (unten), die Geschwindigkeit des Schienenfahrzeugs (Mitte) und die zurückgelegte Wegstre- 
cke (oben) für einen Zeitraum von 40 s. Links: Darstellung der ferromagnetischen Inhomogenitäten 
transformiert in die Ortsdomäne. Der rote Bereich zeigt exemplarisch einen Abgleich der ferroma- 
gnetischen Signatur (Länge von ca. 25 m) mit der hinterlegten ferromagnetischen Streckenkarte. 
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des Schienenfahrzeugs durch nummerische Integration (Heun-Verfahren) 
ermittelt (vgl. Abbildung 4.1 Ds). Sie beinhaltet die Einzelgeschwindigkei- 
ten ve, über alle Geschwindigkeiten N einer zusammenhängenden Wegstrecke. 
Für die Transformation der ferromagnetischen Inhomogenitäten in die Ortsdo- 
mäne werden zunächst Interpolationsfunktionen für die diskreten Messpunkte 
w und D, erstellt. Die Interpolationsfunktion Ju gibt die ferromagnetischen 
Inhomogenitäten eines Detektors in Abhängigkeit der Zeit an. Die Interpolati- 
onsfunktion fp, gibt die Zeit in Abhängigkeit der zurückgelegten Wegstrecke 
an. Beihohen Geschwindigkeiten des Schienenfahrzeugs ergeben sich mit kubi- 
schen Interpolationen im allgemeinen glattere Kurvenverläufe für die Ortsdar- 
stellung der ferromagnetischen Inhomogenitäten. Bei konstanten Signalstruk- 
turen neigt dieses Verfahren zu Schwingungen. Dies kann weitestgehend durch 
eine stückweise kubische Hermite-Interpolationsfunktion (PCHIP) verhindert 
werden. Die einzelnen Teilsegmente sind stetig differenzierbar und gehen an 
den Rändern nahtlos ineinander über. PCHIP verursacht aber in Bereichen 
mit oszillierenden Signalanteilen Artefakte. Aufgrund der hohen Messrate des 
DIS (30 kHz), verglichen mit der geringen Geschwindigkeit des Schienenfahr- 
zeugs, liegen die zeitdiskreten ferromagnetischen Signaturen der Schienen in 
einer ausreichend hohen Punktdichte vor (vgl. Kapitel 2). Dadurch kann auch 
eine lineare Interpolation verwendet werden, ohne die ferromagnetische Signa- 
tur bei der Transformation in die Ortsdomäne zu verfälschen. 

Über fiin(D;, ds) wird die zurückgelegte Wegstrecke D, in äquidistante Weg- 
punkte d, unterteilt. Mit der Gleichung 4.2 werden die zeitdiskreten ferroma- 
gnetischen Inhomogenitäten in die Ortsdomäne überführt (vgl. Abbildung 4.1). 


D' = füw(fpt(fin(d, ds))) (4.2) 


Die ferromagnetischen Inhomogenitäten des zweiten Detektors können eben- 
falls über die Gleichung 4.2 zu einer ferromagnetischen Streckenkarte verar- 
beitet werden. Für die Kartierungsfahrt ist die Befahrrichtung der Strecke nicht 
relevant. Einzig die Reihenfolge der ferromagnetischen Strukturen kehrt sich bei 
einer entgegengesetzten Fahrt auf derselben Strecke um. Im Liniennetzplan wird 
deshalb die jeweilige Befahrrichtung bei der Kartierungsfahrt mit angegeben. 
Dadurch können beide Befahrrichtungen mit demselben Liniennetzplan abge- 
deckt werden. In der Abbildung B.1 bis B.4 sind die ortsdiskreten ferromagneti- 
schen Inhomogenitäten entlang der zurückgelegten Wegstrecke für die Fahrma- 
növer aus der Abbildung 3.10 zu sehen. Ein Teilstück eines Liniennetzplans ist 
in der Abbildung 4.2 abgebildet. Die Linien in schwarz und cyan repräsentieren 
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Abbildung 4.2: Exemplarischer Aufbau eines Liniennetzplans mit einer ferromagnetischen Stre- 
ckenkarte des Streckenabschnittes [A, C, E, F, M] für die linke und rechte Schiene (schwarz, 
cyan). 


die Schienen der Gleisabschnitte (A bis M) und werden über Weichen (rote Drei- 
ecke) miteinander verbunden. Ein Gleiswechsel kann nur über eine der beiden 
angrenzenden Weichen erfolgen. Unter dem Liniennetzplan sind die ortsdiskre- 
ten ferromagnetischen Inhomogenitäten der Streckenabschnitte [A, C, E, F, M] 
für die linke und rechte Schiene (schwarz und cyan) abgebildet. Durch die 
Separation der Merkmalspunkte (ferromagnetischen Inhomogenitäten) von den 
eigentlichen Positionskoordinaten werden Mehrdeutigkeiten bei der Gleisbe- 
legung, die aufgrund von Messunsicherheiten (bspw. bei der GNSS-basierten 
Positionsschätzungen) in Bereichen mehrerer parallel zueinander verlaufender 
Gleisabschnitten zustandekommen (siehe [K], [L], [M] und [N]), verhindert. 
Die räumliche Nähe der einzelnen Schienen zueinander stellt somit mit dem 
DIS kein begrenzendes Element für die Positionsgenauigkeit des Schienenfahr- 
zeugs dar. Der Aufbau der ferromagnetischen Streckenkarte ist so konzipiert, 
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Abbildung 4.3: Darstellung der georeferenzierten ferromagnetischen Streckenkarte der rechten 
Schiene des T4. Offset der Position des Koordinatensystems; Z:32U; E:306076; N:5666477 [36]. 


dass diese mit zusätzlichen Informationen!, bspw. Bahnübergänge, Haltestel- 
len, Briicken, Balisen oder UTM-Koordinaten aus weiteren Sensoren, beliebig 
erweitert werden kann. In der Abbildung 4.3 ist die Zuordnung der ferro- 
magnetischen Inhomogenitäten der rechten Schiene des T4 zu einem UTM- 
Koordinatensystem zu sehen. In der gelben Ebene ist die zurückgelegte Weg- 
strecke mit der ferromagnetischen Streckenkarte der Schiene und der UTM- 
Koordinate der nördlichen Richtung zu sehen. Die rote Ebene zeigt die zurück- 
gelegte Wegstrecke und die UTM-Koordinate in östlicher Richtung. Die blaue 
Ebene zeigt die UTM-Koordinaten (N | E) der ferromagnetischen Streckenkarte 


an. 


! Die zusätzlichen Informationen sind optional und dienen nicht der Entscheidungsfindung für 
eine Schienenzuordnung. 
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4.2 Absolute Positionsbestimmung mit dem DIS 


Für die absolute Positionsbestimmung des Schienenfahrzeugs wird zunächst 
mit dem DIS die Geschwindigkeit anhand der zeitdiskreten ferromagnetischen 
Inhomogenitäten ermittelt. Anschließend werden die ferromagnetischen Inho- 
mogenitäten von der Zeit- in die Ortsdomäne (siehe Kapitel 4.1) überführt. 
Die ortsdiskreten Teilsegmente der ferromagnetischen Streckensegmente wer- 
den ähnlich der Geschwindigkeitsmessung aus Kapitel 3.1 mit dem Template- 
Matching-Verfahren in der ferromagnetischen Streckenkarte gesucht. Bei der 
Positionsbestimmung besteht das Template-Signal aus dem ortsdiskreten Teil- 
segment und das Image-Signal aus der ferromagnetischen Streckenkarte. Der 
Punkt, an dem die gemessenen ferromagnetischen Signaturen die höchste Über- 
einstimmung mit der ferromagnetischen Streckenkarte aufweisen, entspricht 
dem gesuchten Schienenelement und damit der Position des Schienenfahrzeugs 
im Schienennetz. In der Abbildung 4.1 ist für ein Teilsegment mit einer Länge 
von ca. 25 m (roter Bereich) die Ortsdiskretisierung der ferromagnetischen 
Signatur zusammen mit dem Abgleich über die ferromagnetische Strecken- 
karte (blau) zu sehen. Das Verfahren wird anschließend für die nachfolgen- 
den Streckensegmente wiederholt und so eine kontinuierliche Positionsbestim- 
mung des Schienenfahrzeugs erreicht. Mit zunehmender Länge des verwende- 
ten Templates wird die Unsicherheit bei der Positionsbestimmung auf Kosten 
der Rechenzeit verringert. Abhängig von der Größe der hinterlegten ferroma- 
gnetischen Streckenkarte genügt bereits ein kurzes Schienensegment, um die 
Position des Schienenfahrzeugs eindeutig zu bestimmen. Ein entscheidender 
Vorteil des Verfahrens besteht darin, dass im Gegensatz zu IMU oder Raden- 
codern keine initiale Position des Schienenfahrzeugs erforderlich ist. Unsicher- 
heiten, die bei der Transformation der ferromagnetischen Inhomogenitäten von 
der Zeit- in den Ortsbereich auftreten, wirken sich nur lokal auf das aktuelle 
Streckensegment aus. Im Vergleich zur Koppelnavigation akkumulieren sich 
Unsicherheiten durch den Abgleich mit der ferromagnetischen Streckenkarte 
nicht auf. 


61 


4 Kartierung und Positionsbestimmung 


4.3 Aktualisierung der ferromagnetischen 
Streckenkarte 


Für die Erstellung einer ferromagnetischen Streckenkarte ist eine Kartierungs- 
fahrt der Strecke erforderlich. Sind mehrere ferromagnetische Signaturen D. 
aus unterschiedlichen Befahrungen B derselben Strecke vorahnende, können 
diese über das Superpositionsprinzip zu einer gemeinsame ferromagnetischen 
Streckenkarte AD’ kombiniert werden (vgl. Gleichung 4.3). Uber o. können 
die ferromagnetischen Strukturen der einzelnen Befahrungen unterschiedlich 
stark gewichtet werden. Vorausgesetzt die ferromagnetischen Signaturen sind 
bereits zueinander ausgerichtet. 


B B 
AD! = Y anD, KE (4.3) 


n=0 n=0 


Die Ausrichtung kann auch direkt über den Positionsabgleich der ferroma- 
gnetischen Streckkarte geschehen. Die ferromagnetische Streckenkarte wird 
dadurch mit jeder Befahrung aktualisiert. Um eine Fragmentierung der ferro- 
magnetischen Streckenkarte zu vermeiden, ist ein zentrales Verfahren fiir den 
Abgleich erstrebenswert. In der Abbildung 4.4 ist hierfür ein mögliches Verfah- 
ren dargestellt. Abhängig von den Gewichtungen können mit dem Aktualisie- 
rungsverfahren Störungen aus der initialen Streckenkarte reduziert oder auch 
Veränderungen in den ferromagnetischen Strukturen in die ferromagnetische 
Streckenkarte übernommen werden. 


4.4 Versuchsaufbau zur Evaluation der 
Positionsbestimmung mit dem DIS 


Um die Positionsgenauigkeit des ferromagnetischen Lokalisierungsverfahrens 
zu untersuchen, ist eine ferromagnetische Streckenkarte des zu befahrenden 
Gleisabschnittes erforderlich. Im Falle des T4 (vgl. Kapitel 3.3.1) wird diese 
über eine Kartierungsfahrt der Strecke realisiert. Die zeitdiskreten ferromagne- 
tischen Inhomogenitäten werden über die Geschwindigkeit des Schienenfahr- 
zeugs (vgl. Kapitel 3.1 und Gleichung 4.2) in die ortsdomäne überführt und 
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Abbildung 4.4: Ein mögliches Verfahren zur Aktualisierung und Synchronisierung unterschiedli- 
cher ferromagnetischer Streckenkarten. 
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mit den optionalen UTM-Koordinaten aus der Referenzssensorik in eine Stre- 
ckenkarte aufgenommen. Für eine detaillierte Positionsauswertung werden die 
gemessenen Signalstrukturen beider Detektoren in die Streckenkarte (Karte A 
und Karte B) aufgenommen. Im Normalfall ist die Signatur eines Detektors 
ausreichend. Die ferromagnetische Streckenkarte der rechten Schiene des T4 
ist in der Abbildung 4.3 abgebildet. In der Abbildung 4.5 ist der nordöstliche 
Streckenabschnitt mit den dazugehörigen ferromagnetischen Inhomogenitäten 
eines Detektors für die linke und rechte Schiene im Detail zu sehen. 

Bei den anschließenden Fahrmanövern ID 1 bis 33 (siehe Abbildung 3.10) wird 
das ferromagnetische Lokalisierungsverfahren getestet. Dabei werden zunächst 
die zeitdiskreten ferromagnetischen Signaturen der Schiene kontinuierlich über 
die Geschwindigkeit des Schienenfahrzeugs in den Ortsbereich überführt. Die 
Positionsgenauigkeit des ferromagnetischen Lokalisierungsverfahrens wird mit 
Hilfe von zwei Konfigurationen untersucht: 


e Kombinierter Ansatz: 
Beim kombinierten Ansatz wird das zugrundeliegende Funktionsprinzip 
des ferromagnetischen Lokalisierungsverfahrens untersucht. Um etwaige 
systematische Fehler in der Messschaltung des DIS auszuschließen, wird 
die Ortsdiskretisierung der aktuellen ferromagnetischen Inhomogenitäten 
mit der Geschwindigkeit aus der Referenzssensorik durchgeführt. 


Stand-Alone-Ansatz: 

Beim Stand-Alone-Ansatz wird ausschließlich der DIS verwendet. Dazu 
wird zunächst die Geschwindigkeit aus dem zeitlichen Versatz der ferro- 
magnetischen Inhomogenitäten ermittelt, vgl. Kapitel 3. Die Geschwin- 
digkeit aus dem DIS wird anschließend verwendet, um die aktuellen 
ferromagnetischen Signaturen in die Ortsdomäne zu transformieren. 


Für die Positionsbestimmung werden bei beiden Ansätzen die aktuellen orts- 
diskreten ferromagnetischen Signaturen mit den ferromagnetischen Signaturen 
der Streckenkarte verglichen. Stimmen beide Signaturen überein, ist das Schie- 
nensegment und die Position des Schienenfahrzeugs eindeutig gegeben. Für 
die Ermittlung der Positionsgenauigkeit des DIS werden die optionalen UTM- 
Koordinaten der Streckenkarte mit denen der Referenzsensorik verglichen. 
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Abbildung 4.5: Schematische Darstellung der ortsdiskreten ferromagnetischen Inhomogenitäten 
der rechten und linken Schienen des T4. 


4.4.1 Kombinierter Ansatz 


Beim kombinierten Ansatz wird die Position des Schienenfahrzeugs durch 
den DIS, eine ferromagnetische Streckenkarte und die Geschwindigkeit aus 
der Referenzsensorik bestimmt. Letzteres wird verwendet, um die zeitdiskre- 
ten ferromagnetischen Inhomogenität in die Ortsdomäne zu überführen und 
etwaige systematische Fehler in der Messschaltung auszuschließen. Wird eine 
Strecke von einem Meter zurückgelegt, werden die letzten 15 Meter der orts- 
diskreten ferromagnetischen Signatur der Schiene mit der ferromagnetischen 
Streckenkarte abgeglichen. Beim Abgleich mit der Streckenkarte wird auf eine 
Initialisierung der Position oder eine Einschränkung des Suchraums durch eine 
zuvor getätigte Lokalisation verzichtet. Das Lokalisierungsverfahren wird für 
beide Detektorsignale A und B des DIS durchgeführt. In der Abbildung 4.6 ist 
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Abbildung 4.6: Oben: Häufigkeit der Positionsabweichung zwischen dem DIS (mit einer Templa- 
telänge von 15 m) und der Referenzsensorik ausgewertet für die Detektoren A und B. Unten: Häu- 
figkeit der Positionsabweichung zwischen den Detektoren A und B. 


die Häufigkeit der Positionsabweichung des DIS zu der Referenzposition für 
die Fahrmanöver ID = 1-33 (siehe Abbildung 3.10) zu sehen. Für den Detek- 
tor A/B (rot/blau) wird der Abgleich mit der jeweiligen ferromagnetischen 
Streckenkarte durchgeführt, die mit dem dazugehörigen Detektor A/B erstellt 
wurde. In cyan ist die Häufigkeit der Positionsabweichung zwischen den Detek- 
toren A und B zu sehen. Beim kombinierten Ansatz besitzen 99 % aller ermit- 
telten Positionen des DIS eine Positionsabweichung von weniger als 0,40 m zu 
der Referenzposition. Die Ergebnisse sind für die Detektoren A und B identisch. 
Ein ähnliches Verhalten zeigt sich bei der Positionsabweichung zwischen den 
Detektoren. In der Tabelle 4.1 ist eine detaillierte Übersicht der prozentualen 
Häufigkeit der Positionsabweichung des DIS zu der Referenzposition zu sehen. 


4.4.2 Stand-Alone-Ansatz 
Beim Stand-Alone-Ansatz wird ausschließlich der DIS, der Abstand zwi- 


schen den Detektoren und eine ferromagnetische Streckenkarte verwen- 
det, um die Position des Schienenfahrzeugs zu bestimmen. Es werden 
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Abbildung 4.7: Oben: Häufigkeit der Positionsabweichung zwischen dem DIS (beim Stand-Alone- 
Ansatz mit einer Templatelänge von 15 m) und der Referenzposition für die Detektoren A und B. 
Unten: Häufigkeit der Positionsabweichung zwischen den Detektoren A und B. 


keine zusätzlichen Sensorinformationen, initiale Positionen oder Positionen 
aus einem der vorherigen Zeitschritte einbezogen. Der Stand-Alone-Ansatz 
stellt unter diesen Bedingungen das technisch anspruchsvollste Szenario dar. 
Um die ferromagnetischen Inhomogenitäten in die Ortsdomäne zu überfüh- 
ren, wird zunächst aus den zeitdiskreten ferromagnetischen Inhomogenitä- 
ten die Geschwindigkeit (siehe Kapitel 3) des Schienenfahrzeugs bestimmt. 
Ausgehend von den Ergebnissen der Geschwindigkeitsbestimmung aus dem 
Kapitel 3.4 werden für das Template- und Imagesignal die Signallängen von 
4 s bzw. 6 s verwendet. Über die Geschwindigkeit des DIS wird kontinuier- 
lich die zurückgelegte Wegstrecke ermittelt und die zeitdiskreten ferromagne- 
tischen Inhomogenitäten in die Ortsdomäne überführt. Nachdem eine Strecke 
von einem Meter zurückgelegt wird, werden die letzten 15 Meter der orts- 
diskreten ferromagnetischen Signatur der Schiene mit der ferromagnetischen 
Streckenkarte abgeglichen. Beim Abgleich mit der Streckenkarte wird auf eine 
Initialisierung der Position oder eine Einschränkung des Suchraums durch eine 
zuvor getätigte Position verzichtet. Das Lokalisierungsverfahren wird für beide 
Detektorsignale A und B des DIS durchgeführt. In der Abbildung 4.7 ist die 
Häufigkeit der Positionsabweichung des DIS zu der Referenzposition für die 


67 


4 Kartierung und Positionsbestimmung 


10 Stand-Alone-Ansatz Kombinierter Ansatz 
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Abbildung 4.8: Häufigkeit des betragsmäßigen Positionsabstandes ermittelt aus den Absolutpo- 
sitionen der Detektoren (Det. A und B). Links: Die ortsdiskreten ferromagnetischen Streckenseg- 
mente wurden aus der Geschwindigkeit des DIS abgeleitet. Rechts: Die ortsdiskreten ferromagneti- 
sche Streckensegmente wurden aus der Referenzgeschwindigkeit abgeleitet. In beiden Fallen wurde 
der Abgleich für beide Streckenkarten (Kar. A und B) durchgeführt. Der reale Montageabstand der 
Detektoren beträgt 0,6 m (schwarze vertikale Linie). 


Fahrmanöver ID = 1-33 (siehe Abbildung 3.10) zu sehen. Für den Detek- 
tor A/B (rot/blau) wird der Abgleich mit der jeweiligen Streckenkarte durch- 
geführt, die mit dem dazugehörigen Detektor A/B erstellt wurde. In cyan ist 
die Häufigkeit der Positionsabweichung zwischen den Detektoren A und B zu 
sehen. Beim Stand-Alone-Ansatz, bei dem der gesamte Lokalisierungsprozess 
nur von dem DIS getätigt wird, besitzen 91 % aller ermittelten Positionen eine 
Positionsabweichung von unter 0,60 m zu der Referenzposition. Die Ergebnisse 
sind für die Detektoren A und B identisch. Ein ähnliches Verhalten zeigt sich bei 
der Positionsabweichung zwischen den Detektoren. In der Tabelle 4.1 ist eine 
detaillierte Übersicht der prozentualen Häufigkeit der Positionsabweichung des 
DIS zu der Referenzposition zu sehen. 
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4.4 Versuchsaufbau zur Evaluation der Positionsbestimmung mit dem DIS 


Tabelle 4.1: Prozentuale Häufigkeit der Positionsdifferenz zwischen den Detektoren A/B des DIS 
zu der Referenzposition. Stand-Alone-Ansatz: Beim Stand-Alone-Ansatz wurde die Position aus- 
schließlich über den DIS bestimmt. Kombinierter Ansatz: Bei dem kombinierten Ansatz wurde die 
Geschwindigkeit aus der Referenzsensorik verwendet. Die Länge des ferromagnetischen Strecken- 
elementes für den Abgleich mit der Streckenkarte betrug jeweils 15 m. 


Distanz/(m) < 1,0 < 0,6 < 0,5 < 0,4 < 0,3 
Stand-Alone- Detektor A 93 91 89 84 74 
Ansatz Detektor B 93 91 89 84 73 
Kombinierter Detektor A 99 99 99 99 95 
Ansatz Detektor B 99 99 99 99 95 


4.4.3 Relative Positionsgenauigkeit des DIS 


Für eine detailliertere Bewertung der Positionsgenauigkeit des DIS wird der 
Abstand zwischen den Beiden Detektoren A und B über die gemessene fer- 
romagnetische Signatur an den beiden Detektoren bestimmt. Aufgrund des 
räumlichen Abstandes der Detektoren am Schienenfahrzeug unterscheiden sich 
die gemessenen ferromagnetischen Signaturen an den Detektoren voneinander. 
Werden die Positionen über die ferromagnetische Streckenkarte korrekt ermit- 
telt, muss der Abstand der ermittelten Position dem räumlichen Abstand der 
Detektoren entsprechen. In der Abbildung 4.8 ist die Häufigkeit der Positions- 
differenz der Detektoren unter Verwendung der jeweiligen ferromagnetischen 
Streckenkarte A und B zu sehen. Für die ortsdiskreten ferromagnetischen Inho- 
mogenitäten wurde einmal die Geschwindigkeit über den DIS (links; Stand- 
Alone-Ansatz) und einmal die Referenzgeschwindigkeit (rechts; kombinierter 
Ansatz) verwendet. Die schwarze vertikale Linie kennzeichnet den reale Mon- 
tageabstand der beiden Detektoren (0,6 m). Aus den Einzelpositionen des DIS 
konnte für 85 % (Stand-Alone-Ansatz) bzw. 94 % (kombinierter Ansatz) aller 
Positionen der Abstand der Detektoren zueinander korrekt ermittelt werden. Im 
Detail beträgt der ermittelte mittlere Abstand der Detektoren Adref, = 0,63 m, 
O Ref. = 0,03 m bzw. Adpts = 0,63 m, ODIS = 0,05 m. 
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4.4.4 Ergebnisse der ferromagnetischen Lokalisierung 


Für den Positionsabgleich werden die ferromagnetischen Signaturen der letzten 
15 Streckenmeter in der hinterlegten Streckenkarte gesucht. Der Punkt mit der 
höchsten Übereinstimmung entspricht der Position des Schienenfahrzeugs. In 
der Tabelle 4.1 ist eine detaillierte Übersicht der prozentualen Häufigkeit der 
Positionsabweichung des DIS zu der Referenzposition für den Stand-Alone- 
Ansatz und den kombinierten Ansatz zu sehen. Für die Fahrmanöver ID 1 
bis 33 liegt die Positionsabweichung des DIS in derselben Größenordnung 
wie die Positionsunsicherheit der Referenzsensorik, vgl. Abbildung 4.9. Beim 
kombinierten Ansatz ist für 99 % aller ermittelten Positionen des DIS die Posi- 
tionsdifferenz zu der Referenzposition kleiner als 0,40 m. Beim Stand-Alone- 
Ansatz ist für 91 % aller ermittelten Positionen des DIS die Positionsdifferenz 
zu der Referenzposition kleiner als 0,60 m. Aufgrund der häufigen Anfahrvor- 
gänge (vgl. Abbildung 3.10) und der damit verbundenen geringen Signalvaria- 
bilität an den Detektoren ist die Positionsgenauigkeit beim Stand-Alone-Ansatz 
geringer. Mit der Verwendung eines längeren Templates oder einer Einschrän- 
kung der ferromagnetischen Streckenkarte kann dieses Verhalten kompensiert 
werden. Um detailliertere Aussagen über die Positionsgenauigkeit des DIS zu 
erhalten, wird die Positionsdifferenz der beiden Detektoren A und B zueinander 
betrachtet. Durch den Abstand der Detektoren zueinander unterscheiden sich die 
beiden ferromagnetischen Signalstrukturen geringfügig voneinander. Werden 
nun beide Signalstrukturen mit der Streckenkarte verglichen, muss die Positi- 
onsdifferenz der beiden Detektoren dem Abstand der Detektoren zueinander 
entsprechen. In der Abbildung 4.8 sind die Häufigkeiten der Positionsdifferen- 
zen für den Stand-Alone-Ansatz (links) und den kombinierten Ansatzes (rechts) 
zu sehen. Der reale Abstand der Detektoren? (0,60 m) liegt damit innerhalb der 
Standardabweichung von orrer. = 0,03 m bzw. oprys = 0,05 m des ferromagne- 
tischen Lokalisierungsverfahrens des DIS. 

Das ferromagnetische Lokalisierungsverfahren besitzt eine deutlich höhere 
Positionsgenauigkeit als Sensoren, die bisher für die fahrzeugseitige Positions- 
bestimmung im Schienenverkehr zum Einsatz kommen. Aus den Untersuchun- 
gen konnte gezeigt werden, dass der DIS eine Positionsgenauigkeit aufweist, 
die mit denen der verwendeten Referenzsesnorik (Galileo-GNSS und LASER- 
Tachymeter) vergleichbar ist. Dabei ist das hervorstechendste Merkmal die 


? Wird die räumliche Ausdehnung der Detektoren berücksichtigt, beträgt der Abstand der Detek- 
toren (0,60 + 0,03) m. 
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Abbildung 4.9: Histogramm der Häufigkeit der Positionsdifferenz zwischen dem Galileo GNSS- 
Receiver und dem LASER-Tachymeter. 


räumliche Trennung über die ferromagnetischen Signaturen der Schiene, wor- 
über sich eine gleisselektive Position des Schienenfahrzeugs ergibt. 


4.5 Zusammenfassung Lokalisierung 


Die einzigartigen ferromagnetischen Strukturen entlang der Schienen eignen 
sich für die Geschwindigkeits- und für die absolute Positionsbestimmung 
eines Schienenfahrzeugs. Über eine Streckenkarte werden die ferromagneti- 
schen Inhomogenitäten eines jeden Gleiselementes mit einer Ortsposition ver- 
bunden. Durch die Überführung der aktuell gemessenen ferromagnetischen 
Signalverläufe in die Ortsdomäne sowie einem anschließenden Vergleich der 
ferromagnetischen Signaturen mit denen der Streckenkarte, wird mit dem DIS 
die Grundlage für ein neuartiges Verfahren zur absoluten Positionsbestimmung 
für Schienenfahrzeuge geschaffen. Mit den ferromagnetischen Inhomogenitäten 
wird ein komplementäres Messverfahren bereitgestellt, welches die bisherigen 
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Verfahren für die absolute Positionsbestimmung ergänzt. Die ferromagneti- 
schen Inhomogenitäten der Schiene ermöglichen es, die Position eines Schie- 
nenfahrzeugs mit einer Genauigkeit anzugeben, die den Ansprüchen des ETCS? 
für eine Fahrzeugseitige Positionsbestimmung entspricht. Der DIS ist dabei so 
konzipiert, dass eine absolute Positionsbestimmung ohne zusätzliche Sensoren 
durchgeführt werden kann. Ebenfalls kann auf eine Initialisierung der Position 
verzichtet werden. Die Positionsabweichung zwischen LASER-Tachymeter und 
dem Galileo-GNSS-Receiver, vgl. Abbildung 4.9, die als Referenzsensoren ver- 
wendet wurden, liegt in der Größenordnung der Positionsabweichung des DIS 
zu der Referenzposition. Die Positionsgenauigkeit des DIS istdadurch mit denen 
der Referenzsensoren vergleichbar. Im Gegensatz zu GNSS werden räumlich 
nahe beieinander liegende Positionen über die ferromagnetische Signatur der 
Schiene eindeutig voneinander getrennt. Dadurch ist es möglich, die rechte und 
linke Schiene eines Gleissegmentes eindeutig voneinander zu trennen. Insbe- 
sondere für Bereiche mit mehreren parallel zueinander verlaufenden Gleisen 
ist eine eindeutige Positionsbestimmung des Schienenfahrzeugs gegeben. Jeder 
der Detektoren des DIS ist dabei in der Lage, die Position zu bestimmen. Durch 
einen erweiterten Positionsabgleich in Verbindung mit einer zuvor ermittelten 
Positionen lässt sich die Positionsgenauigkeit des Verfahrens weiter steigern 
und die Länge der ferromagnetischen Streckenelemente verringern. Durch die 
Verwendung von längeren ferromagnetischen Signalstrukturen oder die simul- 
tanen Erfassung der ferromagnetischen Inhomogenitäten der rechten und linken 
Schiene kann die Positionsgenauigkeit weiter verbessert werden. Darüberhin- 
aus können Abweichungen in den ferromagnetischen Strukturen aufgrund von 
baulichen Veränderungen an der Infrastruktur über die Streckenkarte erfasst 
werden. 


3 Die Distanzabweichung wird mit +(5 m + 5 % s) angegeben, wobei s der zurückgelegten Distanz 
des Schienenfahrzeugs entspricht [7]. 
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Streckensignatur 


In den vorangegangenen Kapiteln wurde gezeigt, wie mit Hilfe der ferroma- 
gnetischen Inhomogenitäten entlang der Schienen und des DIS die Geschwin- 
digkeit und die absoluten Position eines Schienenfahrzeugs ermittelt werden 
kann. Für eine schnelle und vor allem effiziente Positionsbestimmung über das 
Template-Matching-Verfahren (siehe Kapitel 3.1) soll die Länge des Templates 
möglichst klein gewählt werden. Gleichzeitig muss das Template groß genug 
sein, um Mehrdeutigkeiten aufgrund von ähnlichen Signalstrukturen in großen 
ferromagnetischen Streckenkarten auszuschließen. Unter Berücksichtigung der 
verwendeten ferromagnetischen Streckenkarte wird im folgenden Kapitel die 
maximal benötigte Länge der ferromagnetische Signalstruktur (Template), die 
für eine eindeutige Positionsbestimmung des Schienenfahrzeugs benötigt wird, 
abgeschätzt. Bei der Abschätzung werden Fehlmessungen, bspw. aufgrund von 
Inkonsistenzen in den ferromagnetischen Signalstrukturen, und das Messrau- 
schen an den Detektoren, welches einen eindeutigen Abgleich mit der fer- 
romagnetischen Streckenkarte erschwert, berücksichtigt. Zusätzlich wird ein 
erweitertes Suchverfahren vorgestellt, mit dem die benötigte Rechenzeit zum 
Auffinden korrespondierender Signalstrukturen in ausgedehnten ferromagneti- 
schen Streckenkarten deutlich verringert wird. 
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Abbildung 5.1: Darstellung der reduzierten ferromagnetischen Signalstruktur an einem Detektor. 
Komponenten mit D’ > 0 besitzen den Zustand H (rot) mit der eingeschlossenen Gesamtflä- 
che Ay, Komponenten mit D’ < 0 besitzen den Zustand L (cyan) mit der eingeschlossenen 
Gesamtfläche Ar. 


5.1 Vereinfachung der ferromagnetischen 
Signalstruktur 


Zur Bestimmung der maximalen Templatelänge, die erforderlich ist, um eine 
eindeutige Positionsbestimmung zusammen mit der hinterlegten ferromagneti- 
schen Streckenkarte zu ermöglichen, wird zunächst die Komplexität der ferro- 
magnetischen Signalstruktur reduziert. Dabei muss die Struktur der ferroma- 
gnetischen Inhomogenitäten erhalten bleiben, um über das Template-Matching- 
Verfahren weiterhin die Geschwindigkeit und Position des Schienenfahrzeugs 
bestimmen zu können. Eine Methode zur Reduktion der Komplexität ist in 
der Abbildung 5.1 zu sehen. Die kontinuierliche ferromagnetische Signalstruk- 
tur D’ wird über die Nulldurchgänge: 

= H. D. >0 
D,= E Dn (5.1) 


in einen diskreten Zustand L bzw. H überführt. Wobei H,, bzw. L,, die Inter- 
valllänge zwischen zwei Nulldurchgängen beschreibt. Über die Verkettung der 
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Abbildung 5.2: Oben: Längenhäufigkeitsverteilung der experimentell gewonnenen ferromagneti- 
schen Zustände. Unten: Gegenüberstellung der Längenhäufigkeitsverteilung aus den experimentell 
erhobenen Messdaten (grün) und den synthetisch generierten Zeichenfolgen (rot). Vergleich beider 
Längenverteilungen über den Kolmogorov-Smirnov-Test. 


einzelnen Zustände, bestehend aus L,„ bzw. Hy, bleibt die einzigartige ferroma- 
gnetische Signalstruktur erhalten. Die reduzierte ortsdiskrete Signalstruktur ist 
dann durch: 


D = {Ho, Lo, H1,Lı,...Hn-1,Ln-1,H„,Ln} (5.2) 


gegeben. Über ein Histogramm mit den Klassen {0, 1, 2, . . . } der Zustandslängen 


ergibt sich eine Häufigkeitsverteilung, die der Wahrscheinlichkeitsdichtefunk- 
tion: 


Rae- ep), x>0 (5.3) 


P E 
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(blaue Kurve in der Abbildung 5.2 oben) entspricht. Über einen Parameterfit 
wird der Versatz ly ~ O und der Skalierungsfaktor s = 1,77 ermittelt. Glei- 
chung 5.3 kann dann geschrieben werden als: 


f(x) = A-exp(-xA), x20 (5.4) 


wobei A = l/s die mittlere Anzahl der Nulldurchgänge pro Meter angibt. Die 
durchschnittliche Länge L (in Meter) der beiden Grundzustände H und L ist 
dann über die Gleichung 5.5 und 5.6 gegeben. 


L(H) =—=s (5.5) 


L(L)===s (5.6) 
Auffallig ist, dass die die Verteilungsdichte (Gleichung 5.4) einer Exponential- 
verteilung folgt. Wiederum letztere entspricht bekanntermaßen der Verteilung 
der Abstände zwischen zwei Ereignissen eines Poisson-Prozesses. Dies legt die 
Vermutung nahe, dass die Nulldurchgänge, gemäß eines Poisson-Prozesses, rein 
zufällig auftreten. Aus der exponentiell abfallenden Verteilung geht zudem her- 
vor, dass längere Streckenabschnitte mit einem konstanten ferromagnetischen 
Zustand seltener auftreten als kürzere Streckenabschnitte. Von der Variation 
der ferromagnetischen Signalstrukturen profitiert insbesondere das Template- 
Matching-Verfahren zur Bestimmung der Geschwindigkeit und Position des 
Schienenfahrzeugs. 


5.2 Abschätzung der ferromagnetischen 
Signallänge 


Für die Untersuchung der ferromagnetischen Signalstrukturen, die sich entlang 
der Schiene befinden, wird die Komplexität auf die beiden diskreten Zustände H 
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und L reduziert. Aus der durchschnittlichen Länge der Zustände H und L (siehe 
Gleichung 5.5 und 5.6) und den dazugehörigen Eintrittswahrscheinlichkeiten: 


Pipe 5.7 

Mu (5.7) 
es AL = = 

P(L) = ee 1- P(B), (5.8) 


die sich aus dem Verhältnis der Flächen Ay und Ar (vgl. Abbildung 5.1) 
zusammensetzen, lässt sich die reduzierte Signalstruktur als Bernoulli-Prozess 
beschreiben, wobei von einer festen Abtastrate von !/s ausgegangen wird. Für ein 
seltenes Ereignis, bspw. eine konstante Zeichenfolge (HHHHHH), vgl. Abbil- 
dung 5.2, auf der ferromagnetischen Streckenkarte (Grundgesamtheit) lässt sich 
die Eintrittswahrscheinlichkeit mit Hilfe der Poisson-Verteilung: 


P(X =x) = Si - exp(-A), mit: A=n-P(H)* (5.9) 
ermitteln. Aus der Gleichung 5.9 geht hervor, dass für seltene Ereignisse (Zei- 
chenfolgen mit k > 1) die Wahrscheinlichkeit gering ist, diese in den fer- 
romagnetischen Inhomogenitäten entlang der Gleisstrecke mehrmals (x > 2) 
vorzufinden. Die Eintrittswahrscheinlichkeit für eine beliebige Zeichenfolge, 
bspw. HLLHH. .., ist über: 


P(HLLHH... |Modell) 
= P(H) - P(L|H) - P(LIL) - PHIL) - PCHIH)... 
= P(H)"® . P(L)® (5.10) 


gegeben. Die Eintrittswahrscheinlichkeit hängt dabei nur von der Häufigkeit 
der Zustände (Hg und Lg) und nicht von der Abfolge der Zeichen innerhalb 
der Zeichenfolge oder einer vorherigen Zeichenfolge ab (vgl. Anhang A.5). 
Die synthetische Zeichenfolge (HHLH...LLHL) aus dem Modell besitzen 
damit, bis auf das Abtastintervall s, die gleiche Längenverteilung (vgl. Glei- 
chung 5.2). Die Ähnlichkeit der synthetischen Zeichenfolge aus dem Modell 
mit der Zeichenfolge aus den experimentell gewonnenen Signalstrukturen wird 
mit dem Kolmogorov-Smirnov-Test (KS-Anpassungstest) und einem Signifi- 
kanzniveau von a = 0,2! überprüft [72]. Das Testergebnis zeigt auf, dass keine 


' Die Nullhypothese soll bestätigt werden. 
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signifikanten Unterschiede zwischen der synthetischen und der experimentell 
ermittelten Zeichenfolge bestehen: Z(970k) = 0,02, p = 0,973 (Bestätigung 
der Nullhypothese). Die synthetische Längenverteilung aus dem Modell und 
die experimentell gewonnene Längenverteilung stammen demnach mit einer 
hohen Wahrscheinlichkeit von derselben Verteilungsform ab. Aus der Anzahl 
der enthaltenen Zeichen (Hg bzw. Lg) des H- und des L-Zustandes im Templa- 
tesignal und der Zeichenlänge der Streckenkarte (Kg) ergibt sich die erwartete 
Eintrittshäufigkeit für die Zeichenfolge des Templatesignals: 


N(F) = (Kg - (Hg + Lp - 1)) : (PH) ™ (PL) 


Term A Term B 
< Kp (Pri bis <1 (5.11) 


Für ausgedehnte Streckenkarten ist die Anzahl der Zeichen in der Streckenkarte 
um ein vielfaches Größer als die Anzahl der Zeichen des gesuchten Templa- 
tes (Term A). Für Kg > (Hp + Leg — 1) kann der Term A vernachlässigt werden. 
Die Eintrittshäufigkeit, eine bestimmte Zeichenfolge (Template) in einer sehr 
großen Zeichenfolge (Streckenkarte) zu finden, hängt dann nur noch von den 
Eintrittswahrscheinlichkeiten P(H) und P(L) und der Anzahl der Zeichen Hz 
und Le ab. Durch das Erfüllen der Ungleichung 5.11 ist die Zeichenfolge des 
Templates mit großer Wahrscheinlichkeit höchstens ein Mal in der Zeichen- 
folge der Streckenkarte enthalten, wodurch eine eindeutige Positionszuordnung 
des Schienenfahrzeugs gegeben ist. Mit dem Abtastintervall s ist die Länge des 
Templates (in Meter) durch: 


T, ~ (Lp + Hp) +s (5.12) 


gegeben. Für die Abschätzung der maximalen Templatelänge wird die Zeichen- 
folge mit der höchsten Eintrittswahrscheinlichkeit ermittelt (vgl. Gleichung 5.11 
und Anhang A.5). Unabhängig von der Länge des Templates enthält die Zei- 
chenfolge nur Zeichen des H-Zustandes (siehe Gleichung 5.11 und Abbil- 
dung A.7). Die Gleichung 5.12 vereinfacht sich dadurch und kann geschrieben 
werden als: 


Ts = TL = Hpg Er (5.13) 
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Für eine konservative Abschätzung der erforderlichen maximalen Template- 
länge wird der Skalierungsfaktor s auf den Wert s* = 1,77 — 2,00 m angehoben. 
Die maximale Templatelänge ist dann durch: 


Tt = Hg- s* > T? (5.14) 


max 


gegeben. 


5.2.1 Maximale Templatelänge 


Aus der Zeichenfolge mit der höchsten Eintrittswahrscheinlichkeit, dem Ska- 
lierungsfaktor und der Länge der Streckenkarte ergibt sich die maximale Tem- 
platelänge, die erforderlich ist, um ein Schienenfahrzeug auf einer gegebenen 
Streckenkarte eindeutig zu lokalisieren. In der Abbildung 5.3 ist auf der Ordi- 
natenachse die Länge des Templates- und auf der Abszissenachse die Länge 
der Streckenkarte in Metern aufgetragen. Zu sehen sind die maximal benötigte 
Templatelängen mit den dazugehörigen Restunsicherheiten, die sich über die 
Gleichung 5.11 für unterschiedlich große Streckenkarten ergeben. Die schwarz 
durchgezogene Linie kennzeichnet dabei die untere Grenze der maximalen 
Templatelänge. An dieser Grenze ist die Länge des Templates ausreichend, um 
eine eindeutige Position des Schienenfahrzeugs über die gesamte Streckenkarte 
zu ermitteln. Der gelbe Bereich kennzeichnet die Grenzen, an dem die Restun- 
sicherheit der Position 107° beträgt (vgl. Tabelle A.1). Die beiden vertikalen 
Linien entsprechen einer Streckenkarte, die der Größe des deutschen (grün; 
ca. 61000 km) bzw. des weltweiten Gleisnetzes (schwarz; ca. 1148186 km) ent- 
spricht [73], [74]. In dem grau hinterlegten Bereich kann es durch ein zu klein 
gewähltes Template bzw. eine zu große Streckenkarte zu Mehrdeutigkeiten bei 
der Positionsbestimmung kommen. Eine eindeutige Positionsbestimmung istin 
diesem Bereich nicht mehr gegeben. Für die maximale globale Gleislänge ist 
eine Templatelänge von 200 m erforderlich, um einen eindeutigen Positionsab- 
gleich mit einer Restunsicherheit von 10”? zu erhalten. 
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5.2.2 Speichergröße der ferromagnetischen Streckenkarte 


Für die Speicherung einer unkomprimierten ferromagnetischen Streckenkarte, 
bestehend aus den ferromagnetischen Inhomogenitäten der linken und rech- 
ten Schiene zusammen mit einem UTM32-Koordinatenpaar (E | N), müssen 
ca. 4,00 MB/km bei einer Streckenauflösung von d, = 0,01 m abgespeichert 
werden. In der Tabelle 5.1 sind für unterschiedlich große ferromagnetische 
Streckenkarten die Datenraten und Datenmengen für die Streckenauflösung von 
0,01 m und 0,10 m zu sehen. Das Speichervolumen des gesamten deutschen 
Streckennetzes in seiner höchsten Auflösung (0,01 m) beträgt dabei ca. 260 GB 
und ist durch handelsüblichen Speichermedien (Festplatte oder SD-Karte) rea- 
lisierbar. Die größte Herausforderung bei der absoluten Positionsbestimmung 
besteht darin, möglichst effizient die aktuellen ferromagnetischen Signaturen 
mit denen der ferromagnetischen Streckenkarte abzugleichen. 
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J 10720 
10-24 
10° 10728 
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Länge der Streckenkarte / (m) 
Abbildung 5.3: Templatelänge, die für eine eindeutige Positionsbestimmung erforderlich ist. Die 


vertikalen Linien entsprechen der gesamten deutschen (grün; ca. 61000 km) und weltweiten 
(schwarz; ca. 1148186 km) Gleisstreckenlänge [73], [74]. 
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Tabelle 5.1: Speicherbedarf der unkomprimierten ferromagnetischen Streckenkarten für unter- 
schiedliche Auflösungen und Inhalte. 


ferr. Streckenkarte (MB) Deutschland / (TB) Welt / (TB) 


Œ | N) / Det.A / 
Det.B / Distanz / 4,00 zx 0,26 zx 4,59 
Ad = 0,01 m 


(E| N)/ 
Det.A / Det.B / 3,20 z 0,22 x 3,67 
Ad = 0,01 m 


(E|N)/ 
Det A / Det.B / 0,32 = 0,02 = 0,37 
Ad = 0,10 m 


5.3 Erweitertes Lokalisierungsverfahren 


Abhängig von der Anzahl der Datenpunkte in der ferromagnetischen Strecken- 
karte steigt der Rechenaufwand und dadurch die Zeit, die für ein Positionsab- 
gleich benötigt wird. Unter regulären Betriebsbedingungen ist der Suchraum 
bereits durch eine zuvor ermittelte Position, bspw. über den DIS, Bluetooth-, 
WLAN beacon, GNSS Empfänger oder Funkmasten entlang der Strecke, auf 
ein Teilsegment der ferromagnetischen Streckenkarte beschränkt. Eine unge- 
fähre Position des Schienenfahrzeugs ist dabei ausreichend, um den Suchraum 
und damit die Berechnungszeit auf einen Bruchteil der ursprünglich benötigten 
Zeit zu reduzieren. Ist eine Einschränkung des Suchraumes, bspw. im Fall einer 
initialen Positionsbestimmung, nicht möglich, muss die aktuelle ferromagneti- 
sche Signatur der Schiene mit der gesamten ferromagnetischen Streckenkarte 
abgeglichen werden. Für eine ferromagnetische Streckenkarte mit der Länge von 
10000 km benötigt der Abgleich mit einem Intel Core i7-3740QM @ 2.70 GHz 
zwei Minuten?. 

Um das Auswerteverfahren zu beschleunigen, werden über eine DWT 


? Die Zeit, die für ein Abgleich benötigt wird, hängt entscheidend von der Leistungsfähigkeit 
der verwendeten Hardware ab. Durch eine Aufteilung der ferromagnetischen Streckenkarte in 
mehrere Teilsegmente lässt sich der Suchprozess parallel auf mehreren Prozessorkerne verteilen. 
Die Suchdauer reduziert sich dabei um den Faktor, der der Anzahl an zur Verfügung stehenden 
Prozessorkernen entspricht. 
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die Approximationskoeffizienten der ferromagnetischen Inhomogenitäten 
bestimmt. Dabei verringert sich mit jedem Level die Anzahl der zu verglei- 
chenden Datenpunkte. Werden anstelle der ferromagnetischen Inhomogenitäten 
die entsprechenden Approximationskoeflizienten der Signalstrukturen verwen- 
det, verringert dies die Rechenzeit. Durch den äquivalenten Informationsgehalt 
der beiden Darstellungen sind unterschiedliche Detailgrade in einer gemein- 
same ferromagnetischen Streckenkarte darstellbar, ohne dabei die Menge der 
zu speichernde Datenpunkte zu erhöhen, vgl. Kapitel 2.7.2. 


5.3.1 Wavelet-basiertes Lokalisierungsverfahren 


Das Wavelet-basierte Lokalisierungsverfahren fordert die Reduktion der Mess- 
daten. Anstelle der ferromagnetischen Inhomogenitäten werden für den Stre- 
ckenabgleich die Signalstrukturen der entsprechenden Approximationskoeffi- 
zienten aus dem jeweiligen Level / der DWT als Image- und Template für 
das Template-Matching-Verfahren verwendet (vgl. Kapitel 2.7.3 und 3.1). Der 
Index w;, an dem die Approximationskoeffizienten des Templates die höchste 
Übereinstimmung mit den Approximationskoeffizienten der ferromagnetischen 
Streckenkarte besitzen, wird mit der Gleichung: 


n=2!-w (1 = 0) (5.15) 
n=2!-w, +2) (1 >0) (5.16) 


auf den Positionsindex n der urspriinglichen ferromagnetischen Streckenkarte 
zurückgerechnet (Level /= 0 entspricht der höchsten Auflösung). Die bei der 
Rücktransformation entstehende Unsicherheit AD beträgt je nach Level: 


AD=0 (I = 0) (5.17) 
AD = +2") . d; (1 > 0). (5.18) 


Die maximale Positionsunsicherheit AD max ist dann durch: 
ADmax = [AD] = 2°") . dg (1 >0) (5.19) 
gegeben und symmetrisch um die tatsächliche Position in der ferromagnetischen 


Streckenkarte verteilt. Bei einer Kompressionsrate von bspw. 256 (Level / = 8) 
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entspricht dies einer maximalen Abweichung von ADmax = 1,28 m bei einer 
Auflösung der ferromagnetischen Streckenkarte von d, = 0,01 m. 


5.4 Simulationsaufbau 


Zur Untersuchung des Wavelet-basierten Lokalisierungsverfahrens wird mit 
dem Modell der ferromagnetischen Inhomogenitäten zunächst eine ausgedehnte 
Streckenkarte generiert. Ein Teilsegment aus der Streckenkarte dient dabei 
als Template. Durch eine gezielte Manipulation des Templates werden unter- 
schiedliche Störungen simuliert. Anschließend werden die Approximationsko- 
effizienten aus der Streckenkarte und des Templates ermittelt und miteinander 
verglichen. Aus der Berechnungszeit (Zeitspanne, die für den Abgleich mit der 
Streckenkarte benötigt wird) und der Abweichung zu der Referenzposition wird 
in Abhängigkeit des Level / und der Kompressionsrate k die jeweils maximal 
mögliche Kompressionsrate ermittelt, bei der ein sicherer Abgleich mit der 
Streckenkarte gegeben ist. 


5.4.1 Abgleich idealer Signaturen 


Mit dem Modell der ferromagentischen Inhomogenitäten wird eine synthetische 
Streckenkarte mit einer Länge von ca. 1342 km generiert. Aus der Streckenkarte 
wird ein zufälliges Template mit einer Länge von 200 m ausgewählt. Das Tem- 
plate wird anschließend mit dem Lokalisierungsverfahren der unkomprimierten 
ferromagnetischen Inhomogenitäten (Level 0) und mit dem Lokalisierungsver- 
fahren der Approximationskoeffizienten der Level 1 bis 12 in der Streckenkarte 
gesucht. Für jedes Level werden 100 Suchdurchläufe mit unterschiedlichen 
Templates an unterschiedlichen Positionen durchgeführt. In der Abbildung 5.4 
sind der Mittelwert und die Standardabweichung für die Berechnungszeit (blau) 
und die Positionsabweichung zu der Referenzposition? (rot) für die Approxi- 
mationskoeffizienten (Level 0 bis 12) und der dazugehörigen Kompressionsrate 
zu sehen. Ohne die Abweichung zu der Referenzposition zu erhöhen, reduziert 


3 Die Referenzposition ist aus der zuvor extrahierte Position des Templates bekannt. 
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sich die Berechnungszeit und die zu speichernde Datenmenge der Streckenkarte 
um den Faktor 256. 


5.4.2 Abgleich fehlerbehafteter Signaturen 


Das zum Abgleich mit der Streckenkarte verwendet Template-Matching- 
Verfahren ist robust gegenüber Signalrauschen oder Skalierung der Messergeb- 
nisse. Kommt es jedoch zu baulichen Veränderungen an der Strecke, bspw. zum 
Austausch eines defekten Schienensegments, treten Inkonsistenzen zwischen 
den gemessenen- und den in der Streckenkarte hinterlegten ferromagnetischen 
Strukturen auf. Abhängig von der Länge des veränderten Streckenabschnitts ist 
eine Positionsbestimmung nur noch bedingt gegeben.* Um die Auswirkungen 
von Messfehlern oder lokalen Veränderungen der ferromagnetischen Inhomo- 
genitäten abzuschätzen, wird die in Kapitel 5.4.1 beschriebene Versuchsreihe 
mit einem modifizierte ferromagnetische Signatur als Template wiederholt. 


Kompressionsrate k / (—) 
1 2 4 8 16 32 64 128 256 512 1024 2048 4096 


—— Berechnungszeit 
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Abbildung 5.4: Berechnungszeit und Abweichung zur Referenzposition in Abhangigkeit der dazu- 
gehörigen Wavelet-Koeffizienten für das jeweilige Level 1. Auf einer Streckenkarte von 1342 km 
wurden fiir jeweils 100 unterschiedliche Templates (mit einer Lange von jeweils 200 m) die Posi- 
tionen ermittelt. 


4 Gegebenenfalls muss die ferromagnetische Signatur des erneuerten Streckenabschnittes erneut 
ermittelt werden. 
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Kompressionsrate k / (—) 
1 2 4 8 16 32 64 128 256 512 10242048 4096 


Berechnungszeit t / (s) 
Abweichung AP / (km) 


t T t T t t t t t t t t 
0 1 2 3 4 5 6 7 8 9 10 11 12 
Level 1/(-) 


Abbildung 5.5: Berechnungszeit und Abweichung zur Referenzposition in Abhangigkeit der dazu- 
gehörigen Wavelet-Koeflizienten für das jeweilige Level 1. Auf einer Streckenkarte von 1342 km 
wurden für jeweils 100 unterschiedliche Templates (mit einer Länge von jeweils 200 m) die Posi- 
tionen ermittelt. Für jedes Template wurden an einer zufälligen Positionen im Signal 50 m (25 % 
der Daten) der ursprünglichen Signalsignatur modifiziert. 


Jedem Template wird dabei an einer zufälligen Position ein modifiziertes 50 m 
langes ferromagnetisches Signalstück (entspricht 25 % der Gesamtlänge des 
Templates) zugewiesen. In der Abbildung 5.5 sind der Mittelwert und die Stan- 
dardabweichung für die Berechnungszeit (blau) und die Positionsabweichung 
zu der Referenzposition (rot) für die Approximationskoeffizienten (Level 0 bis 
12) und der dazugehörigen Kompressionsrate zu sehen. Mit den modifizierten 
ferromagnetischen Signaturen des Templates ist weiterhin eine eindeutige Posi- 
tion gegeben. Die Berechnungszeit und die zu speichernde Datenmenge der 
Streckenkarte kann in diesem Fall bis zum Faktor 32 reduziert werden. 


5.4.3 Maximale Templatelänge 


Die maximale Templatelänge, die für eine ferromagnetische Streckenkarte mit 
einer Länge von ca. 1342 km für eine eindeutige Positionsermittlung benötigt 
wird, beträgt 29 m, vgl. Gleichung 5.11. Durch die reduzierte Länge des Tem- 
plates (vgl. Kapitel 5.4.1 und 5.4.2) verringert sich ebenfalls die maximal mög- 
liche Kompressionsrate. Um die Auswirkungen auf die Positionsbestimmung 
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Kompressionsrate E / (—) 
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Abbildung 5.6: Berechnungszeit und Abweichung zur Referenzposition in Abhangigkeit der dazu- 
gehörigen Wavelet-Koeffizienten für das jeweilige Level 1. Auf einer Streckenkarte von 1342 km 
wurden für jeweils 100 unterschiedliche Templates (mit einer Länge von jeweils 29 m) die Positio- 
nen ermittelt. 


zu untersuchen, wird die Versuchsreihe aus Kapitel 5.4.1 mit einer reduzierten 
Templatelänge von 29 m wiederholt. In der Abbildung 5.6 sind der Mittelwert 
und die Standardabweichung für die Berechnungszeit (blau) und die Positions- 
abweichung zu der Referenzposition (rot) für unterschiedliche Approximati- 
onskoeffizienten (Level 0 bis 10) und der dazugehörigen Kompressionsrate zu 
sehen. Ohne die Abweichung zu der Referenzposition zu erhöhen, kann die 
Berechnungszeit und die zu speichernde Datenmenge der Streckenkarte um bis 
zu 32 reduziert werden. 


5.4.4 Ergebnisse aus den Simulationen 


Der initiale Positionsabgleich über die ferromagnetische Signalstruktur der 
Schiene stellt — bedingt durch den rechenintensiven Vorgang — die größte tech- 
nische Herausforderung des Lokalisierungsverfahrens dar. Um den Abgleich 
mit der Streckenkarte zu beschleunigen, werden anstelle der ferromagnetischen 
Inhomogenitäten deren jeweiligen Approximationskoeffizienten aus der Wave- 
lettransformation verwendet. Die simulierten Testszenarien ergeben eine konse- 
quente Reduktion der Auswertezeit und Speicherkapazität. Selbst bei gezielter 
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Tabelle 5.2: Dauer in Sekunden, die für den Abgleich mit der Streckenkarte benötigt wird, um ein 
Template mit einer Länge von 200 m für verschieden lange Streckenkarten mit unterschiedlichen 
Kompressionsraten abzugleichen. Die Auflösung der Streckenkarte beträgt jeweils 0,01 m. 


Kompressionsrate k 1000 km Deutschland Weltweit 
1 x 14,6 zx 490,0 =~ 16766,8 

32 z 0,5 z 15,3 x 524,0 

128 x 0,1 z 3,8 x 131,0 

256 x 0,1 x 2,0 x 65,5 


Manipulation des Templates kann die Position über die Approximationskoef- 
fizienten der Streckenkarte eindeutig ermittelt werden. Die Zeit, die für den 
Abgleich erforderlich ist, kann dabei um das 256-fache reduziert werden, ohne 
dabei die Abweichung zu der Referenzposition zu erhöhen. Die längsseitige 
Unsicherheit, die bei der IDWT auftritt, beträgt dabei 2,56 m. Das Wavelet- 
basierte Kompressionsverfahren liefert dadurch eine effiziente Möglichkeit, den 
Suchprozess auf ausgedehnten Streckenkarten zu beschleunigen, ohne dabei die 
Gleisselektvität durch die ferromagnetische Signatur der Schiene zu verlieren. 
In der Tabelle 5.2 können die benötigten Zeiten und Kompressionsraten abge- 
lesen werden, die erforderlich sind, um auf den gegebenen Streckenkarten eine 
initiale Position des Schienenfahrzeugs ermitteln zu können. 


5.5 Zusammenfassung 


Die ferromagnetischen Inhomogenitäten entlang der Schiene ergeben einzig- 
artige Signalstrukturen, die als Landmarken für die Geschwindigkeits- und 
Positionsbestimmung von Schienenfahrzeugen verwendet werden können. Zur 
Abschätzung der maximalen Templatelänge, die für einen sicheren Abgleich 
mit der Streckenkarte erforderlich ist, wurde auf Basis der bisher gemesse- 
nen Signaturen ein Modell der ferromagnetischen Inhomogenitäten erstellt. 
Aus dem Modell wurde eine maximale Templatelänge von 200 m mit einer 
Restunsicherheit von 107!? ermittelt. Diese Templatelänge erlaubt selbst auf 
ausgedehnten Streckenkarten einen eindeutigen Positionsabgleich. Über das 
Modell können zudem die ferromagnetischen Signalstrukturen für beliebig 
lange Streckenabschnitte simuliert werden. Die synthetischen Daten zeigen 
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auf, wie mithilfe der DWT und den dazugehörigen Approximationskoeffizien- 
ten der Positionsabgleich der ferromagnetischen Signaturen mit der Strecken- 
karte beschleunigt werden kann. Dabei wird eine Reduktion der Rechnezeit 
und der zu speichernden Datenmenge um das 256-fache erreicht, ohne die 
Positionsgenauigkeit oder die Gleisselektivität des Lokalisierungsverfahrens 
zu beeinträchtigen. Somit kann selbst bei einer initialen Positionsbestimmung 
ein Abgleich mit der Streckenkarte in weniger als einer Sekunde erfolgen. 
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6 Zusammenfassung, Fazit und 
Ausblick 


In diesem Kapitel werden die Ergebnisse aus den Untersuchungen des Differenz- 
Induktivitäts-Sensors (DIS) zusammengefasst und interpretiert. Neben den 
Möglichkeiten, die Geschwindigkeit und eine gleisselektive Position über die 
ferromagnetischen Inhomogenitäten der Schiene zu bestimmen, wird aufge- 
zeigt, wie diese für bestehende und zukünftige Anwendungen genutzt werden 
können. Zusätzlich werden Vorteile des DIS gegenüber aktuellen Sensoren, die 
zur fahrzeugseitigen Geschwindigkeits- und Positionsbestimmung eingesetzt 
werden, aufgezeigt. Im Ausblick werden Erweiterungen und Verbesserungen 
an der Hardware des DIS vorgeschlagen. Im Anschluss wird auf die Notwen- 
digkeit weiterer Kartierungs- und Lokalisierungsfahrten eingegangen. 


6.1 Zusammenfassung 


In der vorliegenden Arbeit wurde ein neuartiges Verfahren zur bordautonomen 
Geschwindigkeits- und gleisselektiven Positionsbestimmung von Schienenfahr- 
zeugen vorgestellt. Der eigens für diese Aufgabe entwickelte DIS verwendet 
dazu ausschließlich die bereits vorhandenen ferromagnetischen Infrastruktur- 
elemente (Schienen, Radlenker, Herzstücke etc.) entlang der Gleisstrecke. In 
den ersten Kapiteln wurde der Aufbau, das elektromagnetische Messverfahren 
und die Funktionsweise des DIS erläutert und anschließend anhand von Simula- 
tionen und unterschiedlichen Experimenten validiert. Aus den zeitsynchronen 
Messwerten beider Detektoren konnte die DIS-Arbeitsfrequenz von 30 kHz 
bestätigt werden. Über ein Template-Matching-Verfahren wurde anschließend 
der zeitliche Versatz der zuvor über eine soft-thresholding (wavelet shrinkage) 
gefilterten ferromagnetischen Signaturen der beiden Detektoren bestimmt. Mit 
dem Abstand der Detektoren wurde die Geschwindigkeit des Schienenfahrzeugs 
ermittelt. Über das berührungslose Messverfahren ist die Geschwindigkeit 
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unempfindlich gegenüber Verschleiß und Schlupf zwischen Rad und Schiene. 
Gleichzeitig ergibt sich aus den einzigartigen ferromagnetischen Signaturen, die 
sich entlang der Gleisstrecke befinden, und der ferromagnetischen Streckenkarte 
eine Methode zur gleisselektiven Bestimmung der Position von Schienenfahr- 
zeugen. Die Leistungsfähigkeit des DIS zur Bestimmung der Geschwindigkeit 
und der gleisselektiven Position wurde auf der Teststrecke des PCW (T4) unter- 
sucht. Für den anspruchsvollen Niedriggeschwindigkeitsbereich (< 2 m/s) liegt 
die mittlere Geschwindigkeitsabweichung des DIS zur Referenzgeschwindig- 
keit bei Av = 0,02 m/s. 97 % aller Geschwindigkeiten, die mit dem DIS ermittelt 
wurden, besitzen eine Abweichung zur Referenzgeschwindigkeit von weniger 
als + 0,5 m/s. Der DIS erfüllt dadurch die Anforderungen des ETCS zur 
Geschwindigkeitsbestimmung im Niedriggeschwindigkeitsbereich. Zusätzlich 
zu der Gleisselektivität besitzen 99 % der ermittelten Positionen des DIS eine 
Positionsabweichung zu der Referenzposition von unter 0,60 m. Die relative 
Positionsabweichung beträgt weniger als 0,05 m. Aus den Messergebnissen des 
DIS wurde ein Modell entwickelt, welches die Struktur der ferromagnetischen 
Inhomogenitäten der Schienen wiedergibt. Basierend auf diesem Modell wurde 
eine maximale Templatelänge von 200 m abgeleitet, mit der eine eindeutige 
Position auf einer ferromagnetischen Streckenkarte mit einer Länge von mehr 
als einer Millionen km ermittelt werden kann. Um den Abgleich auf solch 
einer ausgedehnten ferromagnetischen Streckenkarte zu beschleunigen, werden 
anstelle der ferromagnetischen Inhomogenitäten die dazugehörigen Approxi- 
mationskoeffizienten für den Abgleich verwendet. Dadurch wird die Rechenzeit 
und die Datenmenge um das 256-fache reduziert, ohne dabei die Positionsge- 
nauigkeit oder die Gleisselektivität des Verfahrens zu verlieren. 


6.2 Fazit 


Anhand der bisherigen Untersuchungen wurde gezeigt, wie ausschließlich über 
die ferromagnetischen Signaturen der Schiene und mit dem DIS berührungs- 
los die Geschwindigkeit und eine gleisselektive Position des Schienenfahr- 
zeugs ermittelt werden kann. Im Vergleich zu den etablierten Verfahren bietet 
der DIS zusammen mit den ferromagnetischen Inhomogenitäten der Schiene 
diverse Vorteile. Das berührungslose Messverfahren liefert stets die Relativge- 
schwindigkeit zwischen dem Schienenfahrzeug und der Gleisstrecke. Dadurch 
ist das Messverfahren unempfindlich gegenüber Verschleiß des Radumfangs 
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oder Schlupf zwischen Rad und Schiene. Für die gleisselektive Positionsbe- 
stimmung des Schienenfahrzeugs sind weder zusätzliche Infrastrukturelemente 
entlang der Gleisstrecke, noch eine freie Sicht auf Satelliten, klare Sichtverhält- 
nisse oder eine Initialisierung der Position erforderlich. Der Aufbau des DIS 
und das zugrunde liegende Messverfahren sind dabei so konzipiert, dass aus- 
schließlich der Abstand der Detektoren für die Bestimmung der Geschwindig- 
keit und eine ferromagnetische Streckenkarte für die Bestimmung der Position 
erforderlich ist. Ein Abgleich oder eine Kalibrierung der Detektoren unterein- 
ander ist nicht notwendig. Der Aufbau der ferromagnetischen Streckenkarte 
profitiert dabei von der starken Bindung des Schienenfahrzeugs an die Gleis- 
strecke. Die ferromagnetische Streckenkarte ist so konzipiert, dass diese die 
Struktur der Schienen durch einen Liniennetzplan widerspiegelt. Hochauflö- 
sende und geografisch exakte Streckenkarten bzw. eine aufwendige Vermes- 
sung der Gleisstrecke entfallen. Ungeachtet dessen sind die linke und rechte 
Schiene eines Gleissegmentes eindeutig über die ferromagnetischen Inhomoge- 
nitäten voneinander unterscheidbar. Dadurch sind geografisch dicht beieinander 
liegende Streckenabschnitte eindeutig voneinander trennbar. Zur Validierung 
des DIS wurden ein LASER-Tachymeter und ein Galileo GNSS-Receiver als 
Referenzsensoren verwendet. Diese weisen eine wesentlich höhere Genauig- 
keit auf als Sensoren, die üblicherweise für die fahrzeugseitige Bestimmung 
der Geschwindigkeit und Position im Schienenverkehr zum Einsatz kommen. 
Die Positionsabweichung des DIS zu der Referenzsensorik liegt dabei in der 
Größenordnung der Positionsungenauigkeit der verwendeten Referenzsenso- 
ren. Der DIS erfüllt somit die Anforderungen des ETCS für die fahrzeugseitige 
Geschwindigkeits- und Positionsbestimmung. Verglichen mit dem ECS besitzt 
der DIS ein deutlich leistungsfähigeres Verfahren, um die absolute Position des 
Schienenfahrzeugs zu bestimmen. Im Gegensatz zu optischen Messverfahren 
ist das verwendete elektromagnetische Messverfahren des DIS unempfindlich 
gegenüber Verschmutzung der optischen Bauteile oder den vorherrschenden 
Lichtverhältnissen. Die Ergebnisse zeigen deutlich, welches Potenzial die fer- 
romagnetischen Inhomogenitäten der Schiene und der DIS besitzen. Es ergeben 
sich viele Vorzüge, wie eine autarke Arbeitsweise unter der alleinigen Verwen- 
dung des DIS (Stand-Alone-Ansatz), eine berührungslose Erfassung der ferro- 
magnetischen Inhomogenitäten, eine gleisselektive Positionsbestimmung sowie 
ein Messverfahren, dessen Messprinzip sich von den bisherigen Ansätzen und 
Sensoren abgrenzt (vgl. Tabelle 1.1). Hieraus ergeben sich Anwendungen, um 
bestehende Lokalisierungsverfahren zu ergänzen oder völlig neue Verfahren im 
Bahnbereich zu etablieren. 
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6 Zusammenfassung, Fazit und Ausblick 


6.3 Ausblick 


Wie bereits aus den Untersuchungen hervorgeht, besitzen sowohl die linke als 
auch die rechte Schiene eines Gleissegmentes eine einzigartige ferromagne- 
tische Signatur, mit der die Geschwindigkeit und gleisselektive Position des 
Schienenfahrzeugs ermittelt werden kann. Um das bisherige Messverfahren 
weiter zu verbessern, ist eine Erweiterung des DIS von zwei auf vier oder mehr 
Detektorspulen, welche die ferromagnetischen Signaturen beider Schienen zur 
selben Zeit erfassen, erstrebenswert. Durch die zusätzlichen Messungen und die 
Möglichkeit, unterschiedliche Abstände zwischen den Detektorspulen zu rea- 
lisieren, steigt die Signalvariation mit jedem zusätzlichen Detektor an. Davon 
profitiert die Genauigkeit und Robustheit zur Bestimmung der Geschwindigkeit 
und Position des Schienenfahrzeugs. 

Neben der hardwareseitigen Erweiterung des DIS sind zusätzliche Kartierungs- 
fahrten erforderlich, um die bestehenden ferromagnetischen Streckenkarten 
weiter auszubauen. Im Rahmen des Localization Reference Train! (LRT) [8], 
dem Deutschen Zentrum für Luft- und Raumfahrt e. V.? (DLR) und des 
Living Lab? [75] bestehen bereits konkrete Vorhaben, bei denen bereits Stre- 
ckenabschnitte zwischen Annaberg-Buchholz und Schwarzenberg (im Erzge- 
birge) kartiert wurden. Zukünftig ist eine permanente Installation des DIS auf 
einem Schienenfahrzeug vorgesehen, um kontinuierlich die Position des Schie- 
nenfahrzeugs zu ermitteln und die Langzeitstabilität der ferromagnetischen 
Signaturen entlang der Gleisstrecke zu untersuchen. 


I TU Braunschweig, RWTH Aachen, Thales Deutschland GmbH, DLR Institut für Verkehrstech- 
nik, Hörmann Vehicle Engineering Chemnitz, TU Darmstadt. 

2 Institut für Verkehrstechnik. 

3 In Kooperation mit Thales Deutschland GmbH. 
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A Induktivitätsmessung 


A.1 Phasenverschiebung 


Induktivitäten Z und Kapazitäten C an Wechselspannungen verursachen in elek- 
trischen Schaltungen eine Phasenverschiebung y zwischen der Spannung und 
der Stromstärke. Im allgemeinen Fall kann deshalb in der Wechselstromtechnik 
nicht mehr davon ausgegangenen werden, dass die Spannung und Stromstärke 
zu jedem Zeitpunkt dieselbe Phase besitzen. Für den harmonischen Verlauf 
wird der Momentanwert der Spannung U und des Stroms / durch 


U(t) = Uge (A.1) 
I(t) = heite, w=2nfo (A.2) 


mit den Extremwerten Up bzw. Io und der Kreisfrequenz w, dargestellt. Die 
Verschiebung der Phase hat zur Folge, dass die Extremwerte und Nulldurch- 
gänge zu unterschiedlichen Zeitpunkten durchlaufen werden. Ist die Frequenz 
und der Versatz in der elektrischen Schaltung bekannt, lässt sich die Phase 
präzise ermitteln. Für den Fall einer idealen Induktivität (elektrische Spule), 
die an einer idealen Wechselspannungsquelle angeschlossen ist, wird durch die 
zeitliche Änderung der Stromstärke eine Spannung 


dI 
Ur = -L— A.3 
b Er (A.3) 
in der Spule induzierte, die der anliegenden äußeren elektrischen Spannung 
entgegengerichtet ist. Das 2. Kirchhoffsches Gesetz besagt, dass die Summe 
aller Teilspannungen gleich Null sein muss: 


d/ 
a E (A.4) 
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A Induktivitätsmessung 


Im 
wL 


R Re 
Abbildung A.1: Links: Zeigerdiagramm einer realen Spule. Rechts: Ersatzschaltbild einer realen 
Spule. 


Die Gleichung A A geht zusammen mit A.2 über in 


U=iwLI (A.5) 
= wLie@t-¥*?), (A.6) 


Wird Gleichung A.5 durch den Strom J dividiert, ergibt sich der induktive 
Widerstand der Spule 


U 
Z, = 7 =ioL. (AT) 


Aus der imaginären Einheit i (siehe Gleichung A.5 bzw. A.6) folgt, dass für 
eine reine induktive Last die Spannung gegenüber dem Strom um 7/2 = 90° 
vorauseilt. Für eine reale Spule, bestehend aus einem langen Draht der um einen 
Spulenkörper gewickelt wird (siehe Gleichung 2.1), ergibt sich auch immer ein 
reellwertiger ohmscher Widerstand. Die Spule kann dadurch nicht mehr als 
reine Induktivität aufgefasst werden. Das Ersatzschaltbild einer realen Spule ist 
in der Abbildung A.1 zu sehen und besteht aus einem idealen ohmschen Wider- 
stand der mit einer idealen Induktivität in Serie geschaltet ist. Daraus ergibt sich 
ein Gesamtwiderstand bestehend aus einem Realteil (ohmscher Widerstand) R 
und einem komplexen (induktiven Widerstand) Zz (siehe Gleichung A.8). 


Z=R+iwLl. (A.8) 
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A.T Phasenverschiebung 


Aus dem Zeigerdiagramm (siehe Abbildung A 1 links) ergibt sich durch die 
vektorielle Addition von Realteil und Imaginärteil der Scheinwiderstand: 


|Z, | = 4R? + (wL)?. (A.9) 


Der von dem Scheinwiderstand und der reellen Achse eingeschlossene Winkel: 


L 
Inte) = (A.10) 
R 
gibt den Phasenwinkel zwischen der Spannung und der Stromstärke an. Für 
einen festen Widerstand R und eine konstante Kreisfrequenz w lässt sich die 
Phase nur über die Induktivität L bzw. Permeabilitätszahl u, verändern. 


dy Rwo 
aL "e ES 
In der Abbildung A 23 ist die Änderung der Phase © Ar £ über die Gleichung A 11. 
in Abhängigkeit von der Messfrequenz der Spule, zu sehen. Die rote Linie kenn- 
zeichnet den Bereich mit einer Messfrequenz von fo = 30 kHz und einem Seri- 
enwiderstand von R = 10 Q. Unabhängig von der Detektorfrequenz (f >0Hz) 
und dem Serienwiderstand (R > 0 Q) ergibt sich eine maximale Änderung der 
Phase für eine Induktivität der Spule L + 0 H. Durch das Einbringen einer 
zusätzliche Kapazität C (in Serie zu der Spule) lässt sich der Arbeitspunkt der 
maximalen Empfindlichkeit einstellen und dadurch die vorhandene Grundin- 
duktivität Lg kompensieren. Die effektive Induktivität der Spule wird durch die 
Gleichung A.12 beschrieben. 


L=Lgt+AL (A.12) 


Dadurch ist gewährleistet, dass eine kleine Induktivitätsänderung AL zu einer 
großen Phasenverschiebung zwischen der Spannung und der Stromstärke führt. 
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A Induktivitätsmessung 


250 -10° 
200 -10° 
150 -10° 
100 -10° 
50 -10° 
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Ké 10 15 N 
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Abbildung A.2: Änderung der Phase we für verschiedene Induktivitäten L. Die rote Linie markiert 
die Messfrequenz bei 30 kHz. 


A.2 Erzwungene elektrische Schwingung 


Durch das Einbringen einer zusätzlichen Kapazität, die in Serie zu dem Wider- 
stand und der Spule geschaltet wird, ergibt sich ein Serienschwingkreis (siehe 
Abbildung A.3 rechts). Aus den Teilspannungen: 


Ur =I-R (A.13) 
1 
Uc = = fra (A.14) 
dI 
=z=L— Al 
UL J (A.15) 
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A.2 Erzwungene elektrische Schwingung 


Im UL =iwLl 


4 


Abbildung A.3: Links: Zeigerdiagramm fiir den Serienschwingkreis bestehend aus einer realen 
Spule und einem Kondensator. Rechts: Ersatzschaltbild des Serienschwingkreises. 


der harmonischen Eingangsspannung (siehe Gleichung A.1) und dem 2. Kirch- 
hoffschen Gesetz: 


Ur+Uc =U +Uz (A.16) 


ergibt sich durch Differentiation nach der Zeit die Schwingungsgleichung: 


I Rd 1 Uo ; 
—+—— bf =iw—e™ A.17 
ao La te mee 

R 1 Uo ; 
=w +iw = + — =iw—e'” Al 
w +tWwr +77 (A.18) 

. Dn uk 

= iw— cos(p) +isin(p) (A.19) 

ToL 


für den Serienschwingkreis. Nach dem Einschwingvorgang stellt Gleichung A.2 
eine Lösung der seriellen Schwingungsgleichung dar. Werden die zeitlichen 
Ableitungen: 


dI ; 
ES = iwe o) (A.20) 
di ën A 
FEB Tout) 
dei =—-W Toe ? (A.21) 
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A Induktivitätsmessung 


5.0 | r -AL = - 100 nH} 
| |-AL= —50pH 
| -AL = 0 pH 

4.0 — b-AL= 50 pH 


|-AL= 1004H 


Io/(A) 


fo / (kHz) 
Abbildung A.4: Stromamplitude eines Serienschwingkreises (Grundinduktivität LG = 670 uH, 


Serienwiderstand R = 1 Q und Kapazität C = 42 nF) für unterschiedliche Induktivitätsänderun- 
gen AL. 


in die Gleichung A.17 eingesetzt, müssen Realteil und Imaginärteil: 


> 1 U 

Re: ven Se sin(p) (A.22) 
R Upw 

Im: gen A.23 

m Lo = FE O (A.23) 


erfüllt sein. Durch Addition der Quadrate (Gleichungen A.22 und A.23) ist die 
partikuläre Lösung: 


ae (A.24) 


= (A.25) 
VR? + (Xi - Xc)? 


fiir die Stromamplitude des eingeschwungenen Zustandes gegeben. In der 
Abbildung A A sind die Stromamplituden für verschiedene Induktivitäten des 
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A.2 Erzwungene elektrische Schwingung 


Serienschwingkreises zu sehen. Bei einer festen Messfrequenz führt bereits 
eine geringe Änderung der Spuleninduktivität zu einer großen Anderung der 
Stromamplitude. Die Impedanz des Serienschwingkreises ergibt sich aus: 


IZrrel = VR? + (XL - Xoy (A.26) 


und besteht aus dem Blindwiderstand der Induktivität Xz = wL und der Kapa- 
zität Xc = a. Die Gleichung A.25 besitzt die Form des ohmschen Gesetzes. 
Fiir die Resonanzfrequenz: 


(A.27) 


fo (A.28) 


u 2nVLC 


bei der sich X, und Xc gegenseitig aufheben, ist die maximal aufnehmbare 
Stromstärke im Serienschwingkreis durch den Serienwiderstand R begrenzt. 
Die Phase zwischen Spannung und Stromstärke im Serienschwingkreis ist aus 
dem Zeigerdiagramm (siehe Abbildung A.3 links) gegeben: 


wL- wilo+AL)- e XL-Xc 
IECH SE EE (A.29) 
Die Grundinduktivität! Lg der Spule wird durch die zusätzliche Kapazität 
kompensiert. Kommt es nun zu einer Änderung der ferromagnetischen Zusam- 
mensetzung entlang der Schiene, verändert sich die Permeabilitätszahl und 
dadurch die Induktivität der Detektorspule (Induktivitätsänderungen AZ). Mit 
der Arbeitsfrequenz lässt sich die Phase auch als Laufzeitdifferenz (zeitlicher 
Versatz) zwischen Spannung und Stromstärke ausdrücken. Die Gleichung A.29 


geht dann über in: 


D Pu,i _ Pu,i 
wo 27 fo 


z (5) (A.31) 


At 


(A.30) 


R 


! Die Grundinduktivität der Detektorspule sollte so gewählt werden, dass ein möglichst großer 
Arbeitsbereich entsteht. Denkbar wäre die mittlere Induktivität L bzw. mittlere Permeabilität 0. 
entlang der Schiene. 
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A Induktivitätsmessung 


A.3 Bestückungsplan der Messschaltung 


Der Bestückungsplan aus der Abbildung A.5 zeigt die jeweiligen Positionen 
der einzelnen Baugruppen für eine der beiden Messschaltungen (vgl. Abbil- 
dung 2.3). 


Abbildung A.5: Bestückungsplan für einen der beiden Messschaltungen des DIS. 
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A.4 Montageposition des DIS 


AA Montageposition des DIS 


In der technischen Zeichnung der Abbildung A.6 ist die Montageposition des 
DIS und der Detektoren A/B über dem Schienenkopf zu sehen (Zeichnung mit 
freundlicher Genehmigung von Bernd Kröper). 
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Abbildung A.6: Technische Zeichnung des DIS (mit freundlicher freundlicher Genehmigung von 
Bernd Kröper). 
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A Induktivitätsmessung 


A.5 Modell der ferromagnetischen 
Inhomogenitaten 


Die Häufigkeit der Signalstrukturen HHHHHH, HHHHHL, HHHHLH, ..., 
LLLLLL wird über eine Simulation (bestehend aus 10° Zeichen) ermittelt und 
anschließend mit der theoretischen Häufigkeit, die aus dem ferromagnetischen 
Modell ermittelt wurde (siehe Gleichung 5.11), normiert. In der Abbildung A.7 
ist die normierte Häufigkeit aus der Simulation zu sehen. Die Zeichenfolge mit 
der höchsten Eintrittswahrscheinlichkeit besteht demnach nur aus den Zustän- 
den H. 


1.0 ai 


0.8 


Normierte Häufigkeit / (—) 


Abbildung A.7: Normierte Häufigkeit der Signalstrukturen HHHHHH, HHHHHL, HHHHLH, 
... LLLLLL, die mit Hilfe einer Simulation, bestehend aus 10% Zeichen, ermittelt wurde. 
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A.6 Safety Integrity Level SIL 


A.6 Safety Integrity Level SIL 


Für den sicheren Betrieb von Maschinen und Anlagen sind sogenannte Risi- 
koanalysen erforderlich. Zur Quantifizierung der Risikoreduzierung wird mit 
Risikobewertungsmethoden die jeweilige SIL-Anforderung erstellt. In der 
Tabelle A.1 ist die Tolerierbare Gefährdungsrate (THR) und die dazugehörige 
Sicherheitsanforderungsstufe zu sehen [76]. 


Tabelle A.1: Tolerierbare Gefährdungsrate (THR) und Sicherheitsanforderungsstufe (SIL) 
(EN 50129, zitiert nach [77]). 


Tolerierbare Gefährdungsrate (THR) Sicherheitsanforderungsstufe (SIL) 
pro Stunde und pro Funktion 


107° < THR < 1078 4 
1078 < THR < 1077 3 
10-7 < THR < 10° 2 
10°° < THR < 10° 1 
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B Ferromagnetische Signaturen 


Darstellung der ortsdiskreten ferromagnetischen Inhomogenitäten für die Fahr- 
manöver ID = 1 bis 33 (vgl. Abbildung 3.10) zusammen mit den superpositio- 
nierten ferromagnetischen Inhomogenitäten der rechten Schiene (ID = 1-7), der 
linken Schien (ID = 12-33) und der inneren Schiene (ID = 8-11). 
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B Ferromagnetische Signaturen 
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Abbildung B.1: Darstellung der ortsdiskreten ferromagnetischen Inhomogenitäten der Fahrmanö- 
ver ID = 1 bis 7 (vgl. Abbildung 3.10) und der Superposition aus ID = 1-7. 
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Abbildung B.2: Darstellung der ortsdiskreten ferromagnetischen Inhomogenitäten der Fahrmanö- 
ver ID = 8 bis 11 (vgl. Abbildung 3.10) und der Superposition aus ID = 8-11. 
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Abbildung B.3: Darstellung der ortsdiskreten ferromagnetischen Inhomogenitäten der Fahrmanö- 
ver ID = 12 bis 22 (vgl. Abbildung 3.10). 
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Abbildung B.4: Darstellung der ortsdiskreten ferromagnetischen Inhomogenitäten der Fahrmanö- 
ver ID = 23 bis 33 (vgl. Abbildung 3.10) und der Superposition aus ID = 12-33. 
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